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Diversity  of  the  Marine  Cyanobacterium  Trichodesmium: 
Characterization  of  the  Woods  Hole  Culture  Collection  and 
Quantification  of  Field  Populations 

by 

Annette  Michelle  Hynes 

Submitted  to  the  Department  of  Biology 
on  August  28,  2009,  in  partial  fulfillment  of  the 
requirements  for  the  degree  of 
Doctor  of  Philosophy  in  Biological  Oceanography 

Abstract.  Trichodesmium  is  a  colonial,  INb-fixing  cyanobacterium  found  in  tropical 
oceans.  Species  of  Trichodesmium  are  genetically  similar  but  several  species  exist  to¬ 
gether  in  the  same  waters.  In  order  to  coexist,  Trichodesmium  spp.  may  occupy  differ¬ 
ent  niche  spaces  through  differential  utilization  of  resources  such  as  nutrients  and  light, 
and  through  responses  to  physical  characteristics  such  as  temperature  and  turbulence.  To 
investigate  niche  differentiation  in  Trichodesmium ,  I  characterized  cultured  strains  of  Tri¬ 
chodesmium ,  identified  and  enumerated  Trichodesmium  clades  in  the  field,  and  investi¬ 
gated  P  stress  and  N2  fixation  in  field  populations.  Species  of  Trichodesmium  grouped 
into  two  clades  based  on  sequences  from  16S  rDNA,  the  internal  transcribed  spacer  (ITS), 
and  the  heterocyst  differentiation  gene  hetR.  Clade  I  contained  Trichodesmium  erythraeum 
and  Trichodesmium  contortum,  and  clade  II  contained  Trichodesmium  thiebautii,  Tri¬ 
chodesmium  tenue,  Trichodesmium  hildebrandtii,  and  Trichodesmium  pelagicum.  Each 
clade  was  morphologically  diverse,  but  species  within  each  clade  had  similar  pigmen¬ 
tation.  I  developed  a  quantitative  polymerase  chain  reaction  (qPCR)  method  to  distin¬ 
guish  between  these  two  clades.  In  field  populations  of  the  Atlantic  and  Pacific  Oceans, 
the  qPCR  method  revealed  that  clade  II  Trichodesmium  spp.  were  more  prominent  than 
clade  I  in  the  open  ocean.  Concentrations  of  Trichodesmium  did  not  correlate  with  nutri¬ 
ent  concentrations,  but  clade  I  had  wider  temperature  and  depth  distributions  than  clade 
II.  Temperature  and  light  are  physical  characteristics  that  may  define  niche  spaces  for 
species  of  Trichodesmium.  Clade  I  and  II  concentrations  correlated  with  each  other  in  the 
Pacific  but  not  in  the  Atlantic,  indicating  that  the  two  clades  were  limited  by  the  same  fac¬ 
tors  in  the  Pacific  while  different  factors  were  limiting  the  abundance  of  the  two  clades  in 
the  Atlantic.  Trichodesmium  populations  in  the  North  Atlantic  were  more  P  stressed  and 
had  higher  N2  fixation  rates  than  populations  in  the  western  Pacific.  While  nutrient  con¬ 
centrations  didn’t  directly  correlate  with  Trichodesmium  concentrations,  the  contrasting 
nutrient  regimes  found  in  the  Atlantic  and  Pacific  Oceans  might  influence  distributions  of 
the  two  clades  differently.  Unraveling  the  differences  among  species  of  Trichodesmium 
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begins  to  explain  their  coexistence  and  enables  us  to  understand  factors  controlling  global 
N?  fixation. 


Thesis  Supervisor:  Scott  C.  Doney 
Title:  Senior  Scientist 

Thesis  Supervisor:  John  B.  Waterbury 
Title:  Scientist  Emeritus 
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CHAPTER  1 


Introduction 


Statement  of  Purpose 

Trichodesmium  is  a  colonial,  N2-fixing  cyanobacterium  and  a  significant  source  of  fixed 
nitrogen  (N)  in  oligotrophic  oceans.  Its  roles  in  the  biogeochemical  cycling  of  N,  phos¬ 
phorus  (P),  iron  (Fe),  and  carbon  (C)  have  made  it  the  target  of  a  wide  variety  of  physio¬ 
logical,  ecological,  optical,  and  modeling  studies.  Trichodesmium  spp.  form  colonies  visi¬ 
ble  to  the  naked  eye,  making  field  collections  and  observations  easier  than  for  many  other 
phytoplankton.  There  is  currently  a  genome  available  for  Trichodesmium  erythraeum 
IMS  101,  aiding  genetic  studies  of  this  genus.  Teasing  out  the  factors  that  may  lead  to 
niche  differentiation  in  Trichodesmium  addresses  Hutchinson’s  “paradox  of  the  plank¬ 
ton,”  whereby  many  more  species  of  plankton  coexist  than  would  be  expected  based  on 
the  number  of  limiting  resources  available  (Hutchinson,  1961).  Lessons  learned  from 
Trichodesmium  can  be  applied  to  other  species  of  phytoplankton  which  are  not  as  well- 
studied. 

There  are  several  morphologically  distinct  but  genetically  similar  species  of  Tri¬ 
chodesmium  that  coexist  in  tropical  and  subtropical  waters.  I  hypothesized  that  the  co¬ 
existence  of  several  Trichodesmium  species  is  due  in  part  to  niche  differentiation.  Tri¬ 
chodesmium  spp.  occupy  different  niche  spaces  to  utilize  resources  such  as  nutrients  and 
light  and  to  respond  to  physical  characteristics  such  as  temperature  and  turbulence.  To 
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explore  the  niche  spaces  occupied  by  Trichodesmium  spp.,  I  investigated  cultured  strains 
available  in  the  Woods  Hole  culture  collection  and  natural  communities  in  the  Atlantic 
and  Pacific  Oceans.  I  characterized  cultured  Trichodesmium  strains  by  sequencing  the  het¬ 
erocyst  differentiation  gene  hetR ,  16S  rDNA,  and  the  internal  transcribed  spacer  (ITS)  re¬ 
gions;  photographing  and  measuring  trichomes  (filaments);  and  analyzing  the  absorption 
spectra  of  phycobiliproteins  (PBP).  I  used  the  hetR  sequences  to  develop  a  real-time  quan¬ 
titative  polymerase  chain  reaction  (qPCR)  method  to  identify  and  quantify  the  two  major 
clades  of  Trichodesmium  in  cultures  and  field  samples.  I  applied  this  qPCR  method  to  re¬ 
solve  spatial  distributions  of  Trichodesmium  from  samples  collected  with  Niskin  bottles 
across  the  Equatorial  Atlantic,  the  west  Pacific  warm  pool,  and  the  South  Pacific.  Phos¬ 
phorus  stress  was  assayed  using  enzyme-labeled  fluorescence  (ELF)  of  alkaline  phos¬ 
phatase  (AP)  and  N2  fixation  was  measured  using  acetylene  reduction  in  Trichodesmium 
populations  across  the  west  Pacific  warm  pool  and  the  western  North  Atlantic. 

Ecological  roles  of  Trichodesmium 

Cyanobacteria  evolved  3500  Ma  ago  and  dramatically  changed  the  redox  state  of  the 
Earth’s  atmosphere  through  oxygenic  photosynthesis  (Schopf,  2000).  Also  known  as 
blue-green  algae,  cyanobacteria  are  found  in  both  fresh  and  marine  waters  and  a  wide 
variety  of  terrestrial  habitats.  They  are  ubiquitous,  prone  to  blooms,  and  can  produce 
toxins  (Cohen  and  Gurevitz,  2006).  Trichodesmium  spp.  are  filamentous,  colonial,  non- 
heterocystous,  N2-fixing  cyanobacteria  belonging  to  the  order  Oscillatoriales  (Waterbury, 
2006).  The  genus  name  comes  from  the  Greek  root  “trikh-”  meaning  “hair”,  referring  to 
the  filaments  known  as  trichomes.  Trichodesmium  is  found  in  tropical  and  subtropical 
oligotrophic  waters  with  shallow  mixed  layer  depths  and  warm  temperatures. 

Some  cyanobacteria  such  as  Trichodesmium  are  diazotrophs  and  contribute  to  new 
N  inputs  in  oligotrophic  systems  (Capone  et  al.,  1997).  Fixed  N  is  considered  to  be  the 
proximal  limiting  nutrient  in  ocean  ecosystems.  For  a  system  in  steady-state,  new  produc¬ 
tion  is  stimulated  by  inputs  of  new  N,  resulting  in  excess  production  available  for  export 
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Figure  1-1:  Trichodesmium  and  other  diazotrophs  fuel  new  production  through  N2  fixa¬ 
tion.  Photo  shows  a  surface  slick  of  Trichodesmium  in  the  western  South  Pacific  Ocean. 
Inputs  of  new  nitrogen  stimulate  the  growth  of  phytoplankton.  New  production  can  be 
exported  out  of  the  euphoric  zone,  drawing  down  surface  concentrations  of  carbon.  Photo 
courtesy  of  Daniel  Ohnemus.  Illustration  by  E.  Paul  Oberlander,  WHOI. 


(Fig.  1-1)  (Eppley  and  Peterson,  1979).  New  N  comes  from  eddy  diffusion  processes, 
seasonal  deep  mixing,  atmospheric  deposition,  lateral  advection,  and  N2  fixation  (Karl, 
2002).  In  the  North  Atlantic,  N2  fixation  rates  of  Trichodesmium  can  equal  or  exceed  the 
vertical  flux  of  NO3  into  surface  waters  (Capone  et  al.,  2005),  while  in  the  North  Pacific, 
N2  fixation  is  the  source  of  up  to  half  of  the  new  N  (Karl  et  al.,  1997). 

Trichodesmium  is  also  a  significant  contributor  of  fixed  C,  through  its  own  pho¬ 
tosynthesis  and  through  stimulation  of  primary  production  of  other  organisms  (Fig.  1-1). 
Trichodesmium  can  account  for  8-47%  of  total  primary  production  in  the  tropical  North 
Atlantic  (Carpenter  et  al.,  2004).  In  a  global  marine  ecosystem  model,  the  presence  of 
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N2-fixers  increased  primary  production  by  diatoms  and  small  phytoplankton  through  the 
excretion  of  dissolved  organic  nitrogen  (DON)  in  N-limited  areas  (Moore  et  al.,  2002; 
Moore  and  Doney,  2007).  Model  results  by  Coles  et  al.  (2004)  showed  that  N2  fixation 
increases  dissolved  inorganic  nirogen  (DIN),  phytoplankton  biomass,  primary  production 
and  export  flux  in  the  subtropics. 

Factors  affecting  the  growth  of  Trichodesmium  include  nutrients  such  as  P  and 
Fe,  temperature,  light  intensity,  mixed  layer  depth  (MLD),  and  turbulence  (Capone  et  al., 
1997;  Post,  2005).  The  amelioration  of  N  limitation  through  N2  fixation  can  cause  a  draw¬ 
down  of  other  nutrients  such  as  P  and  Fe  (Karl,  2002).  The  North  Atlantic  receives  more 
dust  input  than  the  North  Pacific,  and  the  Fe  available  in  terrestrial  dust  is  necessary  for 
the  enzyme  nitrogenase  used  in  N2  fixation.  Low  excess  N  relative  to  excess  P  in  dis¬ 
solved  nutrients  is  a  necessary  but  not  sufficient  condition  for  selection  for  N2-fixing  or¬ 
ganisms  in  oligotrophic  regions;  in  N-limited  scenarios,  dominant  competitors  such  as 
diatoms  cannot  grow  well  and  diazotrophs  or  diatoms  with  diazotrophic  symbionts  can 
succeed  (Karl,  2002).  When  nutrient  N:P  <  16,  the  classical  Redfield  ratio,  diazotrophs 
are  able  to  compete,  drawing  down  P  while  increasing  N  and  making  the  system  more 
P-limited  (Karl,  2002).  Trichodesmium  spp.  have  been  shown  to  be  able  to  use  organic 
sources  of  P  (Dyhrman  et  al.,  2006;  Sohm  and  Capone,  2006)  and  are  capable  of  lux¬ 
ury  uptake  of  Fe  and  P  (Kustka  et  al.,  2003;  White  et  al.,  2006a).  In  a  numerical  model 
simulating  vertical  migration  of  Trichodesmium ,  colonies  >  1000  /urn  were  shown  to  be 
capable  of  P-mining:  sinking  by  carbohydrate  ballasting  below  the  phosphocline,  tak¬ 
ing  up  P,  and  then  rising  back  up  to  the  euphoric  zone  (White  et  al.,  2006b).  Differential 
adaptations  to  nutrient  limitation  including  alternative  nutrient  sources,  uptake  rates,  half¬ 
saturation  constants,  luxury  uptake,  and  P-mining  can  lead  to  niche  differentiation  or  co¬ 
existence  of  species  limited  by  different  resources  (Titman,  1976). 

In  addition  to  its  roles  in  N2  and  C  fixation,  Trichodesmium  provides  substrate  and 
shelter.  A  variety  of  organisms  ranging  from  bacteria  to  crustacean  larvae  are  enriched 
by  2-5  orders  of  magnitude  in  association  with  Trichodesmium  colonies  (Sheridan  et  al., 
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2002).  The  copepod  Macrosetella  gracilis  uses  Trichodesmium  as  a  substrate  for  juve¬ 
nile  development  as  well  as  a  food  source  (O’Neil,  1998).  Trichodesmium  blooms  can  be 
undesirable  in  regions  like  the  Gulf  of  Mexico,  however,  because  they  can  produce  tox¬ 
ins  and  contribute  little  to  the  marine  food  web,  being  grazed  upon  by  only  a  few  pelagic 
harpacticoid  copepods:  M.  gracilis,  Miracia  efferata,  and  Oculosetella  gracilis  (O’Neil 
and  Roman,  1994;  Post,  2005). 

Species  identifications  and  distributions 

There  are  six  well-described  species  of  Trichodesmium'.  T.  erythraeum,  T.  thiebautii,  Tri¬ 
chodesmium  contortum,  Trichodesmium  tenue,  Trichodesmium  hildebrandtii  and  Tri¬ 
chodesmium  pelagicum  (formerly  Katagnymene  pelagica  and  Katagnymene  spiralis). 
Traditionally,  species  of  filamentous,  non-heterocystous  cyanobacteria  such  as  Trichodesmium 
were  identified  using  morphological  characteristics  including  cell  width,  cell  length,  fila¬ 
ment  spirality,  color,  gas  vacuoles,  presence  or  absence  of  a  sheath,  and  colony  morphol¬ 
ogy  (Table  1.1)  (Janson  et  al.,  1995;  Orcutt  et  al.,  2002;  Waterbury,  2006).  However,  these 
characteristics  can  often  prove  to  be  problematic  and  many  studies  refer  to  field  popu¬ 
lations  merely  as  Trichodesmium  spp.  or  identify  them  by  colony  morphology:  spheri¬ 
cal  “puffs,”  fusiform  “rafts”  or  “tufts,”  and  “bowties”  (Janson  et  al.,  1995).  Morpholog¬ 
ical  traits  used  to  identify  Trichodesmium  spp.  “correlate  poorly”  with  genetic  data,  as 
shown  by  the  clustering  of  the  species  formerly  known  as  K.  pelagica  and  K.  spiralis 
among  species  of  Trichodesmium  (Lundgren  et  al.,  2005).  The  different  species  of  Tri¬ 
chodesmium  are  closely  related.  Sequences  of  the  nitrogenase  gene  nifH,  the  heterocyst 
differentiation  gene  hetR,  16S  rRNA,  and  the  ITS  region  as  well  as  fingerprinting  of  HIP  1 
show  low  genetic  diversity  among  Trichodesmium  (Ben-Porath  et  al.,  1993;  Orcutt  et  al., 

2002;  Lundgren  et  al.,  2005). 
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Table  1.1:  Physical  characteristics  and  colony  morphologies  of  Trichodesmium 
spp.(Janson  et  al.,  1995;  Orcutt  et  al.,  2002;  Lundgren  et  al.,  2005) 


Species 

puff 

tuft 

bowtie 

trichomes 

Cell  width 
Cum) 

Cell  length 
(urn) 

Other 

T.  aureum 

/ 

20-25 

gold 

T.  contortion 

/ 

30-40 

T.  erythraeum 

/ 

/ 

6-9 

4-10 

red-brown 

T.  hildebrandtii 

/ 

/ 

/ 

14-17 

5-12 

T.  pelagicum 

K.  pelagica 

/ 

15-35 

straight  trichome 

K.  spiralis 

/ 

12-23 

4-7 

coiled  trichome 

T.  tenue 

/ 

/ 

/ 

5-7 

10-19 

T.  thiebautii 

/ 

/ 

/ 

7-9 

6-12 

The  distribution  of  Trichodesmium  is  influenced  by  temperature.  Temperature  sets 
an  overall  upper  limit  on  growth  rate  (Eppley,  1972).  Increased  temperature  within  a  vi¬ 
able  range  will  increase  enzyme  activity,  hence  increasing  cellular  functions  such  as  res¬ 
piration;  photosynthesis,  particularly  the  Calvin  Cycle;  and  N?  fixation.  In  the  eastern  At¬ 
lantic,  Trichodesmium  is  found  at  highest  densities  between  0  —  15°N,  with  a  complete 
absence  south  of  30°S  (Tyrrell  et  al.,  2003).  Blooms  are  found  in  surface  waters  with 
temperatures  >  25°C,  but  trichomes  can  be  found  in  deeper  waters  with  temperatures  of 
21  —  23°C  (Karl,  2002). 

While  often  found  concurrently,  the  species  of  Trichodesmium  have  different  dis¬ 
tributions.  Table  1.2  gives  an  overview  of  the  dominant  species  found  in  different  ocean 
basins.  While  abundance  varies  in  time  and  space,  T.  thiebautii  is  generally  most  common 
in  the  Northern  Hemisphere  while  T.  erythraeum  is  most  common  in  the  Great  Barrier 
Reef  (Table  1.2)  (Orcutt  and  Gundersen,  2003;  Bell  et  al.,  2005).  At  the  Bermuda  Atlantic 
Time-series  Station  (BATS),  T.  thiebautii  puffs  are  present  year-round,  but  during  the 
summer  maximum  of  Trichodesmium  abundance,  T.  thiebautii  tufts  are  the  predominant 
morphology,  followed  by  T.  thiebautii  puffs  and  T.  erythraeum  tufts  (Orcutt  and  Gunder¬ 
sen,  2003).  In  addition  to  temporal  and  horizontal  spatial  variation,  colony  morphologies 
show  different  vertical  distributions.  In  the  Gulf  of  Aqaba,  Red  Sea,  tuft  colonies  were 
found  near  the  surface  while  puff  colonies  were  found  in  the  lower  half  of  the  photic  zone 
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(Post  et  al.,  2002).  During  a  trans-Atlantic  survey  with  a  video  plankton  recorder  (VPR), 
tufts  were  more  abundant  in  the  upper  50  m  while  puffs  were  abundant  throughout  the 
euphoric  zone  (Davis  and  McGillicuddy,  2006). 


Table  1.2:  Distribution  of  Trichodesmium  spp. 


Region 

Dominant 

Other  Common 

Rare  Species 

Source 

Species 

Species 

N.  Atlantic 

T.  thiebautii 

T.  erythraeum 

1,2,3 

N.  Pacific 

T.  thiebautii , 

T.  erythraeum 

T.  contortum 

4 

S.  Pacific  (Great  Barrier 
Reef) 

T.  erythraeum 

5 

S.  Pacific  (New  Caledo¬ 

T.  thiebautii , 

T.  erythraeum 

6,7 

nia) 

T.  tenue 

S.  Pacific  (open  ocean) 

T.  thiebautii 

T.  erythraeum, 

T.  pelagicum 

7 

Red  Sea 

T.  thiebautii, 

T.  erythraeum 

T.  tenue 

T.  hildebrandtii 

8 

Western  Indian  Ocean 

T.  erythraeum 

T.  thiebautii,  T.  tenue, 
Trichodesmium  sp. 

T.  contortum 

9 

1  =  O’Neil  (1998),  2  =  Carpenter  et  al.  (2004),  3  =  Capone  et  al.  (2005),  4  =  Letelier  and  Karl  (1996),  5 
=  Bell  and  Fu  (2005),  6  =  Neveux  et  al.  (2006),  7  =  Chapt.  4,  8  =  Post  et  al.  (2002),  9  =  Lugomela  et  al. 
(2002). 


Within  an  ocean  basin,  there  is  high  spatial  and  temporal  variation  in  Trichodesmium 
densities  (Carpenter  et  al.,  2004;  Capone  et  al.,  2005).  Mesoscale  eddies  can  influence 
phytoplankton  communities  by  upwelling  or  downwelling  nutrients  and  affecting  temper¬ 
ature.  In  the  North  Atlantic,  diatoms  and  dinoflagellates  bloom  in  upwelling  mode-water 
eddies  while  cyanobacteria  such  as  Synechococcus  are  prevalent  in  cyclones  (Sweeney 
et  al.,  2003),  and  Trichodesmium  are  associated  with  warm,  downwelling  anti-cyclonic 
eddies  (Davis  and  McGillicuddy,  2006).  Trichodesmium  are  also  associated  with  smaller 
physical  features  such  as  windrows  and  Langmuir  cells  (Carpenter  and  Price,  1977;  Car¬ 
penter  and  Capone,  1992).  Buoyant  Trichodesmium  colonies  and  trichomes  accumulate  at 
the  surface  in  convergence  zones  of  these  features. 
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Competition  theory 


According  to  resource  competition  theory,  the  species  which  is  best  able  to  utilize  and 
deplete  the  limiting  resource  will  exclude  competing  species  in  a  steady- state  environ¬ 
ment  (Titman,  1976).  Species  with  different  limiting  resources  could  coexist  (Tilman, 
1977).  An  environment  at  steady-state  is  predicted  to  be  able  to  hold  only  as  many  com¬ 
peting  species  as  it  has  limiting  resources.  However,  there  are  many  more  species  of  phy¬ 
toplankton  than  there  are  limiting  resources,  resulting  in  Hutchinson’s  “plankton  para¬ 
dox”  (Hutchinson,  1961).  Processes  which  support  coexistence  of  species,  and  hence  ad¬ 
dress  the  plankton  paradox,  include  spatial  and  temporal  heterogeneity,  differing  mortal¬ 
ity  rates,  and  niche  differentiation  along  physical  gradients  as  well  as  resource  gradients 
(Levins,  1979;  Tilman  et  al.,  1981). 

Competition  between  two  species  can  be  illustrated  with  the  Lotka-Volterra  com¬ 
petition  model: 

(n  —  otiiM  —  CC12N2)  (1-1) 

(/2  —  0C2iiVi  —  OC22N2),  (1-2) 

where  Nj  is  the  population  size  of  species  i,  r,-  is  the  maximal  growth  rate  for  species  i,  a„ 
is  the  intraspecific  competition  coefficient  (a„  =  ^  where  K,  is  the  carrying  capacity), 
and  CLij  is  the  competition  coefficient  parameterizing  the  strength  of  effect  of  species  j  on 
species  i  (Kot,  2001).  This  model  incorporates  both  inter-  and  intraspecific  competition 
for  resources.  The  success  of  one  species  over  the  other  depends  on  the  relative  strengths 
of  the  competition  coefficients  a ,y.  Species  1  wins  if  ^  and  ^  >  2L.  In  this 
case,  intraspecific  competition  of  species  1  is  relatively  greater  than  the  competitive  ef¬ 
fect  of  species  2  and  the  effect  of  species  1  on  2  is  greater  than  intraspecific  competition 
of  species  2.  Likewise,  species  2  wins  if  ^  ^  and  7^-  <  7^-.  Coexistence  occurs  if 

7^  <  7^-  and  if  intraspecific  competition  is  greater  than  interspecific  competi- 


1  dN  1 
N 1  dt 
1  dN2 
Ni  dt 
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tion.  Either  species  dominates,  depending  on  initial  conditions  if  ^  ^  and  ^  ^ 

(Shigesada  and  Kawasaki,  1997).  The  Lotka-Volterra  model  is  descriptive  and  assumes 
an  environment  in  equilibrium  where  r,  and  a,/  are  fixed  (Tilman  et  al.,  1981).  However, 
in  dynamic  environments  such  as  the  open  ocean,  coexistence  may  occur.  Trichodesmium 
spp.  r,  and  a,7  may  vary  depending  on  the  physical  and  nutrient  environments. 

A  more  mechanistic  model  of  competition  uses  the  Monod  equation,  similar  to 
Michaelis-Menton  enzyme  kinetics.  Monod’s  equation  relates  the  growth  rate  of  a  popu¬ 
lation  to  the  amount  of  resource  available  in  a  continuous  flow  environment: 


1  dNi 
Ni  dt 

dSi 

dt 


mm 


nSj 


j  \  kij  +  Sj 


D(Sjo-Sj) 


f  NjnSj 

i= t  (kij  + 


(1.3) 

(1.4) 


where  Sj  is  the  external  concentration  of  resource  j,  kij  is  the  half  saturation  constant  for 
species  i  limited  by  resource  j,  D  is  the  dilution  rate  for  a  continuous  flow  system  and  the 
steady-state  growth  rate,  Sjo  is  the  influent  concentration  of  resource  j,  n  is  the  number 
of  species,  and  Yij  is  the  yield  of  species  i  limited  by  resource  j  measured  in  cells  per  re¬ 
source  (Tilman,  1977).  The  growth  of  species  i  is  limited  by  the  resource  which  gives  the 
smallest  per  capita  growth  rate.  The  steady-state  solution  to  the  Monod  model  is  analo¬ 
gous  to  the  steady-state  form  of  the  Lotka-Volterra  model,  but  the  Monod  model  allows 
explicit  examination  of  an  environment  which  is  not  in  steady-state  (Tilman,  1977).  Co¬ 
existence  is  possible  if  each  species  is  limited  by  a  different  resource  and  two  species  lim¬ 
ited  by  the  same  resource  can  coexist  if  (Tilman,  1977).  For  example,  if  one 

species  of  Trichodesmium  is  more  limited  by  P  while  the  other  is  more  limited  by  light  or 
Fe,  they  may  be  able  to  coexist. 

A  third  model  of  resource  competition  is  a  variable  internal  stores  model  (Droop, 
1974).  This  model  also  follows  Michaelis-Menten  kinetics  but  relates  growth  to  internal 
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stores  of  resources,  or  cell  quotas: 


1  dNi 

Nj  dt 

dQij 

dt 

^ I 

dt 


min( - 

]  V  Qij 


V, 


kjj  +  Sj 


ri{Qij  8ij)i 


D(sJo-Sj)-'EmJ- 


i=  i  kij  +  Sj 


(1.5) 

(1.6) 

(1.7) 


where  qij  is  the  minimum  internal  store  of  resource  j  by  species  i,  Qij  is  the  cell  quotient 
of  resource  j  in  species  i,  and  is  the  maximal  uptake  rate  of  resource  j  by  species  i 
(Tilman,  1977).  The  internal  stores  model  is  more  complicated  than  the  Monod  model 
because  it  allows  for  luxury  uptake  of  nutrients  by  phytoplankton.  At  equilibrium,  the  two 
models  are  equivalent,  and  although  they  showed  different  behavior  under  non-equilibrium 
conditions,  the  superior  competitor  in  the  Monod  model  was  generally  the  superior  com¬ 
petitor  in  the  variable  internal  stores  model  (Grover,  1990).  Variable  stores  predicts  three 
outcomes  depending  on  the  flow  rate  D  and  ratio  of  the  two  nutrients:  both  species  limited 
by  nutrient  1,  both  species  limited  by  nutrient  2,  or  coexistence. 

The  above  models  examine  up  to  two  limiting  resources.  When  models  include 
three  or  more  limiting  resources,  oscillations  and  chaos  can  result  in  a  monod-type  model 
(Huisman  and  Weissing,  1999).  The  internal  dynamics  are  complex  enough  to  generate 
non-equilibrium  conditions,  and  several  competing  species  can  persist  on  just  a  few  re¬ 
sources  (Huisman  and  Weissing,  1999).  A  variable  internal  stores  model  also  had  oscilla¬ 
tions  and  chaos  that  allowed  more  species  than  resources  to  coexist,  and  the  oscillations 
were  damped  when  compared  to  those  of  a  Monod  model  (Revilla  and  Weissing,  2008). 
These  chaotic  oscillations  have  also  been  observed  in  long-term  culture  experiments  with 
multiple  trophic  levels  and  constant  physical  conditions  (Beninca  et  al.,  2008).  Physiol¬ 
ogy  and  life  history  also  play  a  role  in  competition  dynamics;  incorporation  of  a  trade-off 
between  competitive  abilities  and  a  cyclic  relation  between  competitive  ability  and  re¬ 
source  content  into  a  competiton  model  with  three  resources  generated  oscillations  which 
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enabled  more  than  100  species  to  coexist  (Huisman  et  al.,  2001).  When  resource  limita¬ 
tion  is  combined  with  physical  factors  such  as  temperature,  pH,  or  salinity,  competitive 
advantages  can  shift  from  one  species  to  another  (Tilman  et  al.,  1981).  Trichodesmium 
spp.  compete  with  each  other  and  with  other  phytoplankton  for  nutrients  such  as  P  and  Fe 
and  for  light.  The  effects  of  pigmentation,  temperature,  and  salinity  can  affect  the  kinet¬ 
ics  of  nutrient  uptake  as  well  as  rates  of  N2  fixation,  potentially  differentiating  niches  for 
Trichodesmium. 

In  addition  to  the  internal  nonequilibrium  dynamics  of  a  competitive  system,  the 
ocean  itself  is  rarely,  if  ever,  in  steady-state  when  considering  the  space  and  time  scales 
of  a  phytoplankter.  Phytoplankton  are  distributed  discretely,  so  mass  conservation  models 
such  as  resource-based  competition  models  may  not  always  be  valid  (Siegel,  1998).  Rates 
of  competition  increase  with  cell  size  and  abundance,  so  in  the  oligotrophic  ocean  where 
abundance  is  low,  competitive  exclusion  may  take  so  many  generations  that  other  factors 
such  as  nutrient  pulses  or  changes  in  physical  properties  become  more  important  in  pop¬ 
ulation  dynamics  than  competition  (Siegel,  1998).  Stochastic  perturbations  (noise)  from 
the  environment  can  cause  shifts  in  model  dynamics  and  decrease  the  predictability  of  a 
system  (Bailey  et  al.,  2004).  Episodic  events  such  as  dust  deposition,  mesoscale  eddies, 
and  turbulent  mixing  can  affect  community  dynamics  of  Trichodesmium.  The  spatial  and 
temporal  heterogeneity  of  the  ocean  as  well  as  relevant  scales  must  be  considered  when 
modeling  competition. 

Niche  differentiation  in  cyanobacteria 

A  niche  is  the  combination  of  the  actions  and  relationships  of  an  organism  in  its  commu¬ 
nity  (Hutchinson,  1957).  In  classical  niche  theory,  overlapping  niches  compete  strongly 
with  each  other  and  differentiation  enables  coexistence.  Genomic  and  culture  studies  of 
unicellular  cyanobacteria  have  shown  evidence  of  niche  differentiation.  Two  closely- 
related  Synechococcus  isolates  from  the  Baltic  Sea  have  different  pigmentation,  red  and 
blue-green,  which  allow  them  to  coexist  by  partitioning  the  light  spectrum  (Stomp  et  al., 
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2004).  Prochlorococcus  has  a  high  light  adapted  strain  (MED4)  that  uses  NH4  as  a  N 
source  and  a  low  light  adapted  strain  (MIT9313)  that  uses  NH4  and  N02  as  N  sources 
(Rocap  et  al.,  2003).  Uncultured  strains  of  Prochlorococcus  in  the  North  Pacific  have 
been  found  to  have  assimilation  genes  for  both  N02  and  N02  in  both  high  and  low  light 
adapted  ecotypes,  so  N  sources  do  not  directly  correspond  to  light  adaptation  (Martiny 
et  al.,  2009).  Differences  in  the  vertical  distribution  of  Prochlorococcus  ecotypes  were 
found  between  the  western  and  eastern  North  Atlantic  using  qPCR,  probe  hybridization, 
and  flow  cytometry  (Zinser  et  al.,  2006).  A  north-south  transect  of  the  Altantic  Ocean 
showed  distributions  of  ecotypes  of  Prochlorococcus  following  patterns  of  temperature, 
light,  and  nutrients  (Johnson  et  al.,  2006).  These  biogeographical  patterns  have  emerged 
in  an  ecosystem  model  of  multiple  functional  groups  with  randomly  generated  parame¬ 
ters,  including  Prochlorococcus- like  functional  groups  (Follows  et  al.,  2007). 

The  ecology  of  field  populations  of  Trichodesmium  varies  according  to  colony 
morphology  and  species.  Puffs  and  tufts  have  different  depth  distributions  (Post  et  al., 
2002;  Davis  and  McGillicuddy,  2006).  Trichodesmium  colony  morphologies  also  respond 
differentially  to  physical  forcing;  at  BATS,  T.  thiebautii  puffs  rather  than  tufts  become  the 
predominant  form  in  surface  waters  when  the  sea  surface  temperature  is  greater  than  26°C 
and  winds  maintain  a  turbulent  mixed  layer  (Orcutt  and  Gundersen,  2003).  Puffs,  tufts 
and  bowties  also  show  varying  activities  in  enzymes  such  as  alkaline  phosphatase,  an  en¬ 
zyme  which  can  cleave  phosphate  from  organic  phosphate  and  is  expressed  under  P  stress 
(Stihl  et  al.,  2001;  Webb  et  al.,  2007;  Hynes  et  al.,  2009).  N2  fixation  rates  can  differ  sig¬ 
nificantly  between  colony  types  locally  (Webb  et  al.,  2007).  Grazing  rates  by  copepods  on 
Trichodesmium  vary  with  species  and  colony  morphology,  presumably  due  to  presence  of 
a  neurotoxin;  M.  gracilis  and  M.  efferata  had  higher  ingestion  rates  on  T.  erythraeum  than 
on  T.  thiebautii  and  higher  ingestion  rates  on  puffs  than  tufts  of  T.  thiebautii  (O’Neil  and 
Roman,  1994). 

Cyanobacteria  have  chlorophyll  a  (Chlr/)  as  the  primary  photo  synthetic  pigment 
with  phycobiliproteins  (PBP)  including  phycoerythrin  (PE),  phycocyanin  (PC),  and  al- 
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lophycocyanon  (AP)  as  accessory  light  harvesting  pigments  (Waterbury  et  al.,  1986). 
Species  of  Trichodesmium  have  different  pigmentation  and  their  colors  range  from  dark 
red-brown  ( T.  erythraeum  IMS  101)  to  light  salmon  ( T.  tenue  tenue)  to  light  green  ( Tri¬ 
chodesmium  sp.  H9-4).  T.  erythraeum  has  more  PE  than  T.  thiebautii  (Carpenter  et  al., 
1993).  These  observations  of  variable  characteristics  across  colony  morphologies  and 
species  of  Trichodesmium  point  to  possible  niche  differentiation. 

Overview  of  thesis  chapters 

Trichodesmium  contributes  a  significant  amount  of  fixed  N  to  marine  ecosystems  (Capone 
et  al.,  1997),  and  measurements  of  Trichodesmium  have  been  used  to  parameterize  N2  fix¬ 
ation  in  a  variety  of  models  (Moore  et  al.,  2004;  Coles  et  al.,  2004).  However,  much  is 
unknown  about  its  ecology  and  physiology.  Mortality  rates  and  causes  of  mortality  are  not 
well  characterized,  distributions  are  patchy,  and  potential  differences  among  species  of 
Trichodesmium  have  not  been  thoroughly  studied  and  are  often  ignored.  A  vast  majority 
of  the  culture  experiments  have  used  T.  erythraeum  IMS  101,  which  belongs  to  a  clade  that 
branches  separately  from  the  numerically  predominant  species  of  T.  thiebautii  (Chaps. 

2,  3).  This  species  is  also  used  to  parameterize  N2  fixation  in  global  ecosystem  models 
(Moore  et  al.,  2004).  To  understand  the  roles  Trichodesmium  plays  in  the  ecosystem  and 
the  factors  that  control  its  distribution,  it  is  important  to  understand  the  differences  among 
the  species  of  Trichodesmium  and  the  potential  for  niche  differentiation.  I  aimed  to  an¬ 
swer  the  questions: 

•  What  species  of  Trichodesmium  do  we  have  in  the  Woods  Hole  culture  collection? 

•  How  are  these  cultured  strains  related  to  each  other? 

•  Do  the  cultured  strains  have  characteristics  which  might  lead  to  niche  differentia¬ 
tion? 

•  How  are  held  populations  of  Trichodesmium  distributed? 

•  What  factors  control  the  distributions  of  Trichodesmium  spp.? 
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In  my  PhD  thesis,  I  characterized  cultured  strains  of  Trichodesmium,  developed  a 
molecular  method  to  distinguish  between  the  two  clades  of  Trichodesmium,  identified  and 
enumerated  Trichodesmium  clades  in  the  field,  and  investigated  P  stress  and  N2  fixation  in 
field  populations. 

Chapter  2:  Characterization.  The  Woods  Hole  Oceanographic  Institution  has 
an  extensive  culture  collection  of  Trichodesmium  spp.,  established  and  maintained  over 
fifteen  years  through  the  work  of  Dr.  John  Waterbury.  This  collection  gave  me  the  unique 
opportunity  to  take  measurements  from  a  variety  of  isolates.  The  genome  of  T.  erythraeum 
IMS  101  is  currently  available,  allowing  for  the  development  of  primers  to  target  genetic 
markers  for  identification  such  as  hetR,  16S  rDNA  and  ITS.  Species  of  Trichodesmium 
vary  by  cell  length,  cell  width,  and  colony  morphology.  Using  phase  contrast  light  mi¬ 
croscopy  and  differential  interference  contrast  (DIC)  microscopy,  I  measured  cell  width 
and  length  and  observed  gas  vesicles  of  cultured  strains.  I  analyzed  the  PBP  absorption 
spectra  for  the  cultured  strains  of  Trichodesmium  and  measured  PUB:PEB  and  PE:PC 
ratios  by  the  relative  peak  heights.  I  found  that  the  cultured  species  of  Trichodesmium 
clustered  into  two  clades:  one  containing  T.  erythraeum  and  T.  contortum,  the  other  con¬ 
taining  T.  thiebautii,  T.  tenue,  T.  pelagicum ,  and  T.  hildebrandtii.  Each  clade  has  diverse 
cell  morphology,  but  distinct  PBP  absorption  spectra. 

Chapter  3:  qPCR.  The  hetR  gene  is  putatively  used  in  cyanobacterial  hetero¬ 
cyst  differentiation  and  may  be  used  in  diazocyte  differentiation  in  Trichodesmium  (El- 
Shehawy  et  al.,  2003).  Using  hetR  sequences  from  Lundgren  et  al.  (2005)  and  cultured 
strains  (Chapt.  2),  I  developed  clade- specific  primers  for  approximately  140  bp  regions 
targeted  for  qPCR.  The  qPCR  assay  was  optimized  using  known  concentrations  of  cul¬ 
tured  strains  and  standardized  with  cloned  PCR  products  from  each  clade.  I  also  tested  a 
variety  of  DNA  extraction  methods  to  minimize  PCR  inhibition. 

Chapter  4:  Quantification  of  field  populations.  I  quantified  field  samples  of 
Trichodesmium  from  depth  profiles  along  horizontal  gradients  of  the  equatorial  Atlantic 
Ocean,  the  west  Pacific  warm  pool,  and  the  South  Pacific  Ocean.  I  collected  water  at  var- 
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ious  depths  including  the  surface  and  deep  chlorophyll  maximum  (DCM)  using  a  rosette 
of  Niskin  bottles.  Trichodesmium  clade  I  and  II  abundances  and  distributions  were  ana¬ 
lyzed  with  ancillary  data  such  as  temperature,  salinity,  depth,  and  nutrient  concentrations. 
Clade  II  was  more  numerous  than  clade  I  in  both  the  Pacific  and  Atlantic  Oceans.  Clade 
I  had  deeper  and  colder  distributions  than  clade  II.  There  was  a  correlation  between  clade 
I  and  clade  II  concentrations  in  the  Pacific  but  not  the  Atlantic.  There  were  no  patterns  or 
Trichodesmium  concentrations  with  respect  to  P  or  Fe. 

Chapter  5:  Phosphorus  stress  and  N2  fixation.  Phosphorus  (P)  can  limit  growth 
and  N2  fixation  rates  of  Tricodesmium.  Under  P  stress,  some  phytoplankton  use  alkaline 
phosphatase  (AP)  to  cleave  phosphate  from  organic  sources  (Stihl  et  al.,  2001).  I  analyzed 
P  stress  and  N2  fixation  rates  of  Trichodesmium  populations  from  the  Atlantic  and  Pa¬ 
cific  Oceans.  I  collected  colonies  of  Trichodesmium  and  analyzed  them  for  endogenous 
AP  activity  using  enzyme-labeled  fluorescence  (ELF)  and  for  nitrogenase  activity  using 
acetylene  reduction.  Both  P  stress  and  N2  fixation  rates  were  higher  in  the  western  North 
Atlantic  than  in  the  Pacific,  indicating  that  factors  other  than  P  were  constraining  N2  fixa¬ 
tion  rates  of  Trichodesmium  in  the  Pacific. 
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CHAPTER  2 


Diversity  of  the  N2-fixing  Cyanobacterium 
Trichodesmium:  Characterization  of  the  Woods 
Hole  Culture  Collection 

ABSTRACT: 

The  filamentous,  colonial  cyanobacterium  Trichodesmium  has  six  well-described 
species,  but  many  more  names.  Traditional  classification  was  based  on  morphological 
characteristics  such  as  cell  width  and  length,  gas  vesicle  distribution,  and  colony  mor¬ 
phology,  which  can  lead  to  inconsistent  results.  I  used  the  Woods  Hole  culture  collec¬ 
tion  of  Trichodesmium  to  identify  cultured  strains  by  species  using  cell  morphology; 
phycobiliprotein  absorption  spectra;  and  sequences  of  16S  rDNA,  the  16S-23S  inter¬ 
nal  transcribed  spacer  (ITS),  and  the  heterocyst  differentiation  gene  hetR.  There  were 
two  major  clades  of  Trichodesmium'.  clade  I  consisting  of  Trichodesmium  erythraeum 
and  Trichodesmium  contortum,  and  clade  II  made  up  of  Trichodesmium  thiebautii,  Tri¬ 
chodesmium  tenue,  Trichodesmium  pelagicum,  and  Trichodesmium  hildebrandtii.  These 
clades  were  genetically  coherent  with  similar  phycobiliprotein  composition,  but  morpho¬ 
logically  diverse.  In  the  continual  revision  of  cyanobacterial  taxonomy,  genetic  and  bio¬ 
chemical  information  are  useful  and  informative  complements  to  morphology  for  inform¬ 
ing  classification  schemes. 
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Introduction 


Identification  of  Trichodesmium.  Trichodesmium  spp.  (Cyanobacteriales,  order  Oscilla- 
toriales)  live  in  tropical  and  subtropical  oceans.  This  filamentous  phytoplankter  exists  as 
both  single  trichomes  and  colonies  that  are  visible  to  the  naked  eye,  earning  the  nickname 
“sea  sawdust”  (Cook,  1842).  They  occasionally  form  brown  or  red  blooms  that  are  visible 
from  space  and  gave  the  Red  Sea  its  name  (Ehrenberg,  1830).  In  addition  to  contributing 
fixed  carbon  (C)  to  the  ecosystem,  Trichodesmium  is  a  significant  contributor  of  new  ni¬ 
trogen  (N)  through  N2  fixation  in  oligotrophic  regions  (Capone  et  al.,  1997).  While  only 
a  few  grazers  feed  on  Trichodesmium  (O’Neil  and  Roman,  1994),  colonies  offer  habitat 
by  providing  an  oasis  in  the  open  ocean,  harboring  a  community  of  heterotrophic  bacteria 
and  invertebrates  (Sheridan  et  al.,  2002;  Hewson  et  al.,  2009). 

Classical  identification  of  Trichodesmium  was  based  on  cell  width  and  length, 
sheath  characteristics,  distribution  of  gas  vesicles,  and  colony  morphology.  There  are  six 
well-described  species  of  Trichodesmium :  Trichodesmium  contortum,  Trichodesmium 
erythraeum,  Trichodesmium  hildebrandtii,  Trichodesmium  pelagicum  (formerly  Katag- 
nymene  spiralis  and  Katagnymene  pelagica ),  Trichodesmium  tenue,  and  Trichodesmium 
thiebautii.  These  species  of  Trichodesmium  are  genetically  similar,  but  morphologically 
distinct.  Classification  by  cell  and  colony  morphologies  can  be  confusing  and  mislead¬ 
ing  due  to  variable  and  overlapping  characteristics  (Anagnostidis  and  Komarek,  1988; 
Janson  et  al.,  1995).  Many  genera  in  Oscillatoriaceae  were  identified  by  characteristics 
of  the  sheath,  which  can  vary  with  environmental  conditions  and  proved  to  be  an  unreli¬ 
able  criterion  (Rippka  et  al.,  1979;  Hoffmann,  1988;  Anagnostidis  and  Komarek,  1988). 
Trichodesmium  colonies  come  in  a  variety  of  morphologies  including  spherical  puffs, 
fusiform  tufts,  and  bowties.  Some  species  can  have  more  than  one  morphology,  for  ex¬ 
ample,  T.  thiebautii  forms  both  puff  and  tuft  colonies,  and  different  species  can  have 
the  same  colony  morphology:  T.  thiebautii ,  T.  hildebrandtii ,  and  T.  erythraeum  all  form 
tufts.  Single  colonies  may  also  contain  a  variety  of  filaments,  indicating  that  they  are  not 
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clonal  and  suggesting  that  colonies  may  coalesce  from  single  trichomes  of  several  species 
(Hynes  et  al.,  2009). 

Even  more  perplexing  is  the  list  of  names  for  species  currently  or  formerly  known 
as  “ Trichodesmium .”  Many  species  originally  identified  as  Trichodesmium  were  brought 
under  the  umbrella  of  the  genus  Oscillatoria  and  then  later  were  separated  back  to  Tri¬ 
chodesmium  (Geitler,  1932;  Rippkaet  al.,  1979;  Anagnostidis  and  Komarek,  1988).  Species 
of  the  genus  Katagnymene  were  found  to  be  genetically  similar  to  T.  thiebautii  with  re¬ 
spect  to  the  nitrogenase  gene  nifH,  the  16S-23S  internal  transcribed  spacer  (ITS),  and 
the  heterocyst  differentiation  gene  hetR,  so  Katagnymene  spp.  have  been  included  in  the 
genus  Trichodesmium  (Lundgren  et  al.,  2001;  Orcutt  et  al.,  2002;  Lundgren  et  al.,  2005). 
Due  to  the  inconsistent  nature  of  cell  morphology,  a  large  number  of  synonyms  have  been 
in  use  (Guiry  and  Guiry,  2008).  A  selection  of  synonyms  is  summarized  in  Table  2.1. 

Previous  studies  have  grouped  Trichodesmium  spp.  into  a  variety  of  taxonomic 
schemes.  The  classification  of  species  based  on  morphometric  characteristics  and  cell 
structure  resulted  in  two  groups  (Fig.  2-1  A):  (1)  species  with  peripheral  gas  vesicles  (T. 
tenue,  a  spherical-shaped  colony  from  an  undetermined  Trichodesmium  sp.,  and  T.  ery- 
thraeum);  and  (2)  species  with  randomly  dispersed  gas  vesicles  (T.  thiebautii,  T.  hilde- 
brandtii,  and  T.  contortum)  (Janson  et  al.,  1995).  In  phylogenetic  trees  of  Trichodesmium , 

T.  erythraeum  branched  off  from  all  the  other  groups  (Orcutt  et  al.,  2002).  Trees  con¬ 
structed  from  16S  rRNA  and  hetR  sequences  of  field  and  culture  samples  have  shown 
four  clades:  (1)  T.  hildebrandii,  T.  thiebautii,  and  T.  pelagicum;  (2)  a  novel  morphotype 
named  Trichodesmium  aureum;  (3)  T.  erythraeum  and  a  large,  dark-pigmented  trichome; 
and  (4)  T.  contortum  and  T.  tenue  (Lundgren  et  al.,  2005).  This  grouping  is  represented  in 
the  dendrogram  in  Fig.  2- IB. 

Photosynthetic  pigments.  Members  of  the  genus  Trichodesmium  have  a  photosyn¬ 
thetic  system  typical  of  other  cyanobacteria.  Their  primary  photo  synthetic  pigment  is 
chlorophyll  a  (Chla)  with  a  light-harvesting  phycobilisome  made  up  of  phycobiliproteins 
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Table  2.1:  List  of  species  names  for  Trichodesmium,  taken  partly  from  Guiry  and  Guiry 
(2008). 


Species 

Synonyms 

Ref. 

Trichodesmium  aureum 

1 

Trichodesmium  clevei 

Pelagothrix  clevei  J.  Schmidt 

2,3 

* Trichodesmium  contortion 

Xanthotrichum  contortion 

4 

Skujaella  contorta 

5 

* Trichodesmium  erythraeum 

Oscillaria  erythraea  (Ehrenberg) 

6 

Oscillatoria  erythraea  (Ehrenberg) 

7 

Trichodesmium  ehrenbergii  Montagne 

5 

Trichodesmium  erythraeum  Ehrenberg 

8 

Trichodesmium  flosaquae 

5 

Trichodesmium  havanum 

1 

* Trichodesmium  hildebrandtii 

Oscillatoria  hildenbrandtii 

7 

Skujaella  hildebrandtii 

5 

Trichodesmium  ehrenbergii  f  indicum 

2 

Trichodesmium  hildenbrandtii 

9 

Trichodesm  i urn  indicum 

5 

Trichodesmium  hindsii 

7 

Trichodesmium  iwanoffianum  Nygaard 

Oscillatoria  lacustris 

2 

Trichodesmium  lacustre 

Oscillatoria  lacustris 

2 

Skujaella  lacustris 

5 

Trichodesmium  lenticulare 

Haliarachne  lenticularis 

7 

Trichodesmium  maccii 

10 

* Trichodesmium  pelagicum 

Katagnymene  pelagica 

1,7 

Katagnymene  spiralis 

7 

Tri chodesm  i urn  sp i rali s 

11 

Trichodesmium  scoboideum 

Skujaella  scoboidea 

5 

* Trichodesmium  tenue 

Trichodesmium  radians 

2 

Heliotrichum  radians 

12 

* Trichodesmium  thiebautii 

Oscillatoria  thiebautii 

7 

Skujaella  thiebautii 

5 

Heliotrichum  radians 

13 

*Species  represented  in  the  Woods  Hole  culture  collection  (past  and  present). 

1  =  Lundgren  et  al.  (2005),  2  =  Anagnostidis  and  Komarek  (1988),  3  =  Silva  et  al.  (1996),  4  = 
Schiitt  (1893),  5  =  Drouet  (1968),  6  =  Kutzing  (1845),  7  =  Geitler  (1932),  8  =  Hallegraeff  et  al. 
(2003),  9  =  Gomont  (1892),  10  =  Li  et  al.  (1984),  11  =  Dyhrman  et  al.  (2006),  12  =  Golubic 
(1977),  13  =  Wille  (1904) 
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T.  tenue 

T.  erythraeum 

T.  thiebautii 
T.  hildebrandtii 
T.  contortum 
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T.  tenue 
T.  contortum 
T.  erythraeum 
T.  hildebrandtii 


T.  thiebautii 

T.  pelagicum 
T.  aureum 


Figure  2-1:  Dendrograms  of  classification  groupings  of  Trichodesmium  in  previous 
studies.  (A)  Grouping  based  on  morphological  characteristics  (Janson  et  al.,  1995).  (B) 
Grouping  based  on  hetR  sequences  (Lundgren  et  al.,  2005). 


Figure  2-2:  Typical  absorption  peak  values  for  photosynthetic  pigments  in  cyanobacteria: 
chlorophyll  a  (Chla),  phycourobilin  (PUB),  phycoerythrobilin  (PEB),  phycocyanin  (PC), 
allophycocyanin  (AP)  (Rippka  et  al.,  1974). 


(Glazer,  1987).  Trichodesmium  phycobilisomes  are  dominated  by  the  phycobiliprotein 
phycoerythrin  (PE),  which  has  absorbance  peaks  at  490  -  500  nm  for  its  brown  phycouro¬ 
bilin  (PUB)  component  and  at  545  -  565  nm  for  its  red  phycoerythrobilin  (PEB)  compo¬ 
nent.  The  phycobilisome  also  contains  the  blue-green  phycocyanin  (PC),  which  absorbs 
around  620  nm,  and  the  blue  allophycocyanin  (AP),  which  absorbs  at  650  nm  (Fig.  2-2). 
In  contrast,  Chkz  has  absorption  peaks  at  about  440  nm  and  680  nm  (Rippka  et  al.,  1974). 

Chromatic  adaptation  has  been  described  in  field  populations  and  cultured  Tri¬ 
chodesmium  (Subramaniam  et  al.,  1999;  Bell  and  Fu,  2005).  PUB:PEB  ratios  in  Tri¬ 
chodesmium  have  been  shown  to  increase  at  midday  in  the  Caribbean  Sea  (Subrama- 
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niam  et  al.,  1999).  Great  Barrier  Reef  strain  T.  erythraeum  GBRTRLI101  showed  chro¬ 
matic  adaptation  by  changing  its  PE:Chki  ratio  as  well  as  its  PE:PC  ratio;  PE  increased 
under  green  light,  PC  increased  under  red  light,  and  Chla  increased  under  blue  and  red 
light  (Bell  and  Fu,  2005).  However,  chromatic  adaptation  has  not  been  found  in  other 
field  studies  (McCarthy  and  Carpenter,  1979;  Neveux  et  al.,  2006).  Pigmentation  of  Tri- 
chodesmium  in  the  field  and  in  culture  varies  considerably  with  colors  ranging  from  red  to 
brown  to  green.  This  color  variation  is  partly  due  to  phycobiliprotein  composition;  Car¬ 
penter  et  al.  (1993)  noted  that  T.  erythraeum  had  more  phycoerythrin  per  colony  and  a 
higher  PE:Chk/  than  T.  thiebautii.  Color  may  also  be  influenced  by  growth  state;  senes¬ 
cent  Trichodesmium  cells  appear  more  brown  or  green  than  healthy  cells  due  to  chlorosis. 

Characterization  of  the  Woods  Hole  culture  collection.  The  Woods  Hole  Oceano¬ 
graphic  Institution  (WHOI)  has  an  extensive  collection  of  Trichodesmium  cultures  from 
all  over  the  world  (Table  2.2).  To  classify  these  strains  and  understand  the  relationships 
among  the  various  species,  I  characterized  the  culture  collection  using  cell  morphology, 
phycobiliprotein  spectra,  and  sequences  of  16S  rDNA,  ITS,  and  hetR.  The  results  of  this 
study  showed  that  there  are  two  major  clades  of  Trichodesmium  based  on  sequence  analy¬ 
sis.  Surprisingly,  the  clades  were  not  morphologically  homogeneous,  but  were  consistent 
with  respect  to  phycobiliprotein  content. 
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Table  2.2:  Collection  details  of  the  Woods  Hole  culture  collection  of  Trichodesmium 

Species  Strain  Origin  Latitude  Longitude  Isolated  by:  Year  isolated 

T.  erythraeum  IMS  101  N.  Carolina  coast.  Gulf  Stream  1  1992 

St.  5-1A  Northern  Sargasso  Sea  2 
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Methods 


Culture  conditions.  Batch  cultures  were  grown  at  26°C,  50  //E  m~2s_1,  under  soft 
white  fluorescent  bulbs  with  a  14:10  light-dark  cycle  in  RMP  medium  made  from  75% 
sterile  Sargasso  Sea  water  amended  with  trace  metals  and  vitamin  B 12,  as  described  in 
Webb  et  al.  (2001),  and  with  2  x  10~6  M  EDTA,  15  x  10~6  M  phosphoric  acid,  1  x  10~7 
M  Fe  (ferric  citrate).  Cultures  were  maintained  in  Citranox-  and  acid-washed  250  mL 
polycarbonate  baffled  flasks  closed  with  cotton  stoppers  or  75  mL  square  polycarbonate 
culture  bottles  (Nalgene). 

Micrographs.  Wet  mounts  of  culture  were  examined  using  a  Zeiss  Axioplan  2  micro¬ 
scope  using  phase  contrast  with  a  Plan-Neofluar  lOOx  objective  lens  and  differential  in- 
tereference  contrast  (DIC)  using  an  a  Plan-Fluar  lOOx  objective  lens.  Micrographs  were 
taken  using  a  Zeiss  Axiocam  HRC  digital  camera  and  Axiovision  4.6.3  software.  Images 
were  measured,  and  cell  metrics  were  determined  using  the  MATLAB  image  processing 
toolbox. 

Phycobiliprotein  extraction  and  absorption  spectra.  Trichodesmium  phycobilipro- 
tein  extracts  were  prepared  as  described  in  Rippka  et  al.  (1974)  with  minor  modifications. 
Cultures  were  filtered  on  25  mm,  5  /urn  polycarbonate  filters  and  stored  at  -20°C.  Two 
mL  of  pH  7.0,  5  mmol  L^1  phosphate  buffer  containing  0.1  mol  L^1  NaCl  were  added  to 
each  sample  and  vortexed  to  remove  cells  from  the  filter.  Cells  were  lysed  at  4°C  using  a 
digital  sonifier  (Branson;  3x  for  1  min  at  65%  with  on/off  cycles  for  0.5  s  each),  a  tissue 
homogenizer,  or  a  small  French  pressure  cell  (Aminco).  The  extracts  were  centrifuged  in 
a  Beckman  ultracentrifuge  for  45  minutes  at  1 12,000  ref  and  4  °C.  The  supernatant  was 
pipetted  into  a  1.5  mL  microcentrifuge  tube  and  centrifuged  in  an  Eppendorf  581  OR  cen¬ 
trifuge  for  15  minutes  at  20,000  ref  and  4°C.  Seventy-five  qL  of  supernatant  was  pipetted 
into  a  quartz  microcuvette  and  the  absorption  spectra  from  400  -  700  nm  (0.5  nm  step) 
was  measured  using  a  Shimadzu  UV-1601  spectrophotometer.  Blank  spectra  of  the  phos- 
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phate  buffer  alone  were  subtracted  from  the  absorption  spectra  of  the  samples.  Replicate 
spectra  were  averaged  and  normalized  to  the  maximum  absorption. 

Sequences  and  phytogeny.  Cultured  strains  of  Trichodesmium  were  filtered  on  25  mm, 

5  n m  polycarbonate  filters  and  stored  in  liquid  N2  until  extraction.  DNA  was  extracted  us¬ 
ing  the  DNeasy  Tissue  Kit  (Qiagen),  following  the  manufacturer’s  instructions  for  animal 
tissue.  Extracted  DNA  was  stored  at  -20°C.  All  polymerase  chain  reactions  (PCR)  were 
performed  on  a  Biorad  iCycler  with  25  /uL  reactions.  Primers  are  summarized  in  Table 
2.3.  Amplification  of  16S  rDNA  was  performed  using  the  high-fidelity  iProof  kit  (Biorad) 
to  get  a  662  bp  fragment.  Each  reaction  contained  1  x  HF  Buffer,  200  qM  dNTPs,  0.5  qM 
each  of  primers  CYA-106F  and  CYA-781R(a)  (Niibel  et  al.,  1997),  0.02  U  iProof  DNA 
polymerase,  and  4  ng  of  genomic  DNA.  Reactions  were  cycled  using  the  temperature  pro¬ 
file  of  98°C  for  30  s;  30  cycles  of  98°C  for  10  s,  65.5°C  for  20  s,  and  72°C  for  20s;  and 
one  cycle  at  72°C  for  7  min.  The  ITS  region  was  amplified  using  the  MasterThg  Kit  (Ep- 
pendorf,  5  Prime)  to  get  a  1026  bp  fragment.  Reactions  consisted  of  1  x  Th^Master  PCR 
Enhancer;  1  x  Taq  Buffer;  3.5  mM  Mg2+  (including  Mg2+  in  the  Tag  Buffer);  200  //M 
dNTPs;  0.5  /jM  each  of  primer  tril6S-1247F  (Orcutt  et  al.,  2002)  and  primer  tri-23SR, 
modified  from  23S-241R  (Rocap  et  al.,  2002)  to  be  specific  to  Trichodesmium  based  on 
the  T.  erythraeum  IMS  101  genome;  1.25  U  Taq  DNA  polymerase;  and  4  ng  genomic 
DNA.  Thermal  cycling  conditions  for  ITS  were  95°  for  2  min;  30  cycles  of  95°  for  1  min, 
55°  for  1  min,  and  72°for  1  min;  and  a  final  cycle  of  72°  for  10  min.  Amplification  of 
hetR  was  conducted  using  the  iProof  kit  (Biorad)  to  get  a  448  bp  fragment.  The  reac¬ 
tion  mixture  was  made  up  of  1  x  HF  Buffer,  200  /rmol  L”1  dNTP,  0.5  /jmol  L_1  each  of 
primers  PHI  and  PH2  (Lundgren  et  al.,  2005),  and  0.02  U  g L  1  iProof  DNA  polymerase, 
and  4ng  of  genomic  DNA.  Cycling  conditions  for  hetR  were  98°  for  30  s;  30  cycles  of  98° 
for  10  s,  57.7°  for  20  s,  and  72°  for  20  s;  and  one  cycle  of  72°  for  7  min. 
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Table  2.3:  Summary  of  PCR  primers. 


Gene 

Primers 

Sequence  (5’  to  3’) 

Source 

16S 

For 

CYA-106F 

CGGACGGGTGAGTAACGCGTGA 

(Niibel  et  al.,  1997) 

Rev 

CYA-781R(a) 

GACTACTGGGGTATCTAATCCCATT 

(Niibel  et  al.,  1997) 

ITS  (external) 

For 

tril6S-1247F 

CGTACTACAATGGTTGGG 

(Orcutt  et  al.,  2002) 

Rev 

tri-23SR 

TTCGCTCACCGCTACA 

This  study 

ITS  (internal) 

For 

AlaF 

TWTAGCTCAGTTGGTAGAG 

(Rocap  et  al.,  2002) 

Rev 

AlaR 

CTCTACCAACTGAGCTAWA 

(Rocap  et  al.,  2002) 

hetR 

For 

PHI 

TG  Y  GCKATTTAYATGACCTA 

(Lundgren  et  al.,  2005) 

Rev 

PH2 

ATG  A  AN  GGTATKCCCC  A  A  GGA 

(Lundgren  et  al.,  2005) 

Five  /rL  of  each  PCR  product  was  separated  on  a  1%  agarose  gel,  stained  with 
ethidium  bromide,  and  photographed  using  a  Gel  Logic  200  Imaging  System  and  Molec¬ 
ular  Imaging  Software  (Kodak).  16S,  hetR,  and  ITS  samples  with  only  one  band  were 
purified  using  a  QIAquick  PCR  Purification  Kit  (Qiagen).  Samples  of  ITS  that  showed 
multiple  bands  were  separated  on  1%  agarose  gels  stained  with  ethidium  bromide  and 
purified  using  the  QIAquick  Gel  Extraction  Kit  (Qiagen).  Sequencing  of  purified  PCR 
products  in  both  directions  was  performed  in  duplicate  or  triplicate  for  each  primer  on  a 
3730XL  capillary  sequencer  (Applied  Biosystems)  by  the  Keck  Facility  at  the  Josephine 
Bay  Paul  Center,  Marine  Biological  Laboratory,  Woods  Hole.  BigDye  Terminator  v3.1 

i  lh 

chemistry  was  used  in  ^  reactions:  6  n L  total  volume  with  lx  Sequencing  Buffer,  5 
pmol  primer  (Table  2.3),  0.5  /uL  BigDyeTerminator  v3.1,  0.1  /uL  DMSO,  and  10  -  30  ng 
purified  PCR  product.  Thermal  cycling  conditions  were  60  cycles  at  96°C  for  10  s,  50°C 
for  5  s,  and  60°  for  4  min.  Sequencing  reactions  were  precipitated  with  isopropanol  and 
resuspended  in  5  /uL  Hi-Di  formamide  (Applied  Biosystems).  To  get  complete  coverage 
of  the  ITS  region,  internal  primers  AlaF  and  AlaR,  were  also  used  (Rocap  et  al.,  2002).  A 
few  of  the  16S  and  hetR  purified  PCR  products  were  sent  to  Northwoods  DNA,  Inc.  for 
sequencing. 

Sequences  were  edited  using  Sequencher  4.6  and  aligned  using  either  ClustalX 
or  the  MATLAB  Bioinformatics  Toolbox.  Alignments  were  edited  in  either  Mac  Vec¬ 
tor  8.0  or  the  MATLAB  multiple  alignment  viewer  and  similarity  matrices  were  calcu¬ 
lated,  rounding  to  the  nearest  integer.  Sequences  acquired  from  GenBank  are  listed  be- 
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low  with  accession  numbers  in  parenthesis  (Benson  et  al.,  2005).  The  16S  alignment  was 
trimmed  to  564  bp  and  included  field  samples  T.  hildebrandtii  (AF091322),  T.  contor¬ 
tion  (AF013028),  and  T.  tenue  (AF013029)  from  Janson  et  al.  (1999);  the  Woods  Hole 
cultured  strain  T.  pelagicum  ZK  (AF5 18769);  field  sample  Trichodesmium  havanum 
(AF5 18770)  from  Lundgren  et  al.  (2005);  and  cultured  strain  Oscillatoria  sancta  PCC 
7515  (AY768400).  The  ITS  alignment  was  trimmed  to  exclude  the  16S  and  23S  sequences 
(497  bp).  Trichodesmium  ITS  includes  the  alanine  and  isoleucine  tRNA  gene  sequences, 
which  were  identical  among  the  strains  and  included  in  the  phylogeny.  The  alignment  in¬ 
cluded  Woods  Hole  cultured  strains  T.  pelagicum  ZK  (AF399652)  and  T.  hildebrandtii  II- 
4  (AF399650)  from  Orcutt  et  al.  (2002),  and  cultured  strain  O.  sancta  PCC  7515  (EF178272). 
The  hetR  alignment  was  trimmed  to  447  bp  and  included  field  samples  T.  pelagicum 
(AF490696),  T.  hildebrandtii  (AF490679),  Trichodesmium  aureum  (AF490680),  and  T. 
aureum  (AF490682)  from  Lundgren  et  al.  (2005);  field  samples  T.  contortum  (AF013031), 

T.  contortum  (AF490685),  and  T.  tenue  (AF013033)  from  Janson  et  al.  (1999);  and  Arthrospira 
platensis  HZ01  (EU427544).  16S,  ITS,  and  hetR  sequences  were  concatenated  to  form 
a  1507  bp  sequence  (Gadagkar  et  al.,  2005).  The  concatenated  alignment  included  Gen- 
Bank  sequences  for  field  samples  of  T.  hildebrandtii  (AF490679,  AF399650,  and  AF490679) 
and  the  16S  sequence  for  cultured  strain  T.  pelagicum  ZK  (AF5 18769).  Phylogenetic  trees 
were  created  for  each  genetic  marker  and  the  concatenated  alignment  in  Paup  4.0b  10 
using  parsimony,  neighbor-joining,  and  maximum  likelihood  methods  and  rooted  us¬ 
ing  the  sequences  from  the  closely-related  O.  sancta  or  A.  platensis  (order  Oscillatori- 
ales).  Each  tree  was  bootstrapped  2000  times  except  for  the  ITS  maximum  likelihood 
tree,  which  proved  to  be  computationally  intensive.  Trees  were  edited  using  Mesquite  2.6 
(http://mesquiteproject.org). 

Results 

Genetic  characterization  of  Trichodesmium.  Two  distinct  clades  of  Trichodesmium 
could  be  distinguished  using  each  of  the  three  genetic  markers:  clade  I  containing  T.  ery- 
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thraeum  and  T.  contortum,  and  clade  II  containing  T.  thiebautii,  T.  tenue,  T.  pelagicum, 
and  T.  hildebrandtii  (Figs.  2-3  -  2-14). 

The  partial  16S  sequences  were  highly  similar  for  all  strains  of  Trichodesmium 
(>  97%).  Similarity  to  the  outgroup  O.  sancta  ranged  from  95  -  97%.  For  clade  I,  intra- 
clade  similarities  were  between  99  and  100  %.  Subgroups  of  clade  I  found  consistently 
in  all  three  phylogenetic  analyses  were  (1)  T.  contortum  strains  20-70  and  21-74  and  T. 
havanum ;  and  (2)  T.  erythraeum  strains  XIII-13,  St.  5-1A,  St.  6-1,  St.  6-2,  and  St.  6-5 
(Figs.  2-3  -  2-5).  The  16S  sequence  for  T.  pelagicum  20-71  was  not  successful.  Within 
clade  II,  sequence  similarities  were  100%  with  the  exception  of  T.  hildebrandtii ,  which 
was  99%  similar  to  all  the  other  sequences  in  the  clade.  Trichodesmium  sp.  H9-4,  T.  tenue 
tenue,  and  T.  pelagicum  ZK  formed  a  consistent  subgroup  of  clade  II  and  had  one  base- 
pair  difference  from  the  other  sequences  in  the  clade  (Figs.  2-3  -  2-5).  Field  samples  of  T. 
contortum  and  T.  tenue  from  Janson  et  al.  (1999)  formed  their  own  clade  (Figs.  2-3  -  2-5). 

The  ITS  region  is  >  89%  identical  among  all  strains  of  Trichodesmium  and  simi¬ 
larity  to  the  outgroup  O.  sancta  was  80  -  81%.  Similarities  between  the  two  major  clades 
of  Trichodesmium  were  89  -  90%.  The  three  phylogenetic  analyses  showed  slightly  dif¬ 
ferent  topologies.  The  clade  I  sequences  were  >  98%  identical.  The  neighbor-joining  and 
parsimony  analyses  had  a  subgroup  of  clade  I  containing  T.  erythraeum  strains  St.  5-1  A, 
St.  6-1,  St.  6-2,  and  St.  6-5  (Figs.  2-6,  2-7),  but  the  maximum  likelihood  analysis  did  not 
(Fig.  2-8).  The  maximum  likelihood  analysis  yielded  a  clade  I  subgroup  containing  T. 
erythraeum  strains  IMS  101,  21-75,  K-02  #2,  GBRTRLI101,  and  GBRTRLIN201  (Fig. 
2-8).  Sequences  within  clade  II  were  98  -  100%  identical.  All  three  trees  showed  a  sub¬ 
group  of  clade  II  containing  T.  thiebautii  strains  II-3,  VI- 1  and  St.  3-4  (Figs.  2-6  -  2-8). 
Only  the  neighbor-joining  analysis  contained  a  clade  II  subgroup  of  Trichodesmium  sp. 
H9-4,  T.  tenue  tenue,  and  T.  pelagicum  ZK  (Fig.  2-6). 

Trichodesmium  hetR  sequences  were  >  90%  identical  and  were  74  -  75%  simi¬ 
lar  to  the  outgroup  Arthrospira  platensis.  Inter-clade  similarity  between  clades  I  and  II 
was  90  -  91%.  Members  of  clade  I  were  100%  identical  except  for  T.  erythraeum  21-75, 
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which  was  99%  similar  to  the  other  strains.  All  three  phylogenetic  analyses  clustered  T. 
erythraeum  strains  K- 11  #131,  St.  5-1  A,  St.  6-1,  St.  6-2,  and  St.  6-5  together  as  a  sub¬ 
group  of  clade  I  (Figs.  2-9-2-11).  Intra-clade  similarity  in  clade  II  was  97  -  100%.  All 
three  trees  had  subgroups  of  clade  II  consisting  of  (1)  T.  thiebautii  strains  St.  3-4  and  VI- 
1,  and  (2)  T.  thiebautii  strains  St.  1-2  and  St.  5-1JW  and  T.  hildebrandtii  (Figs.  2-9-2-11). 
The  neighbor-joining  tree  also  showed  a  clade  II  subgroup  containing  T.  thiebautii  II-3, 

T.  tenue  tenue,  and  Trichodesmium  H9-4  (Fig.  2-9).  The  Janson  et  al.  (1999)  field  sam¬ 
ples  of  T.  contortum  and  T.  tenue  formed  their  own  clade  as  did  the  Lundgren  et  al.  (2005) 
field  samples  of  T.  aureum. 

The  concatenated  Trichodesmium  \6S-YTS-hetR  sequences  were  >  93%  identical. 
Inter-clade  similarities  were  93%  and  intra-clade  similarities  were  99  -  100%.  In  clade  I, 
all  three  phylogenetic  analyses  showed  subgroups  of  (1)  T.  erythraeum  strains  St.  5-1A, 
St.  6-1,  St.  6-2  and  St.  6-5;  and  (2)  T.  contortum  strains  20-70  and  21-74  (Figs.  2-12— 
2-14).  Within  clade  II,  the  subgroups  (1)  T.  tenue  tenue  and  Trichodesmium  sp.  H9-4, 
and  (2)  Trichodesmium  strains  St.  1-2  and  St.  5-1JW  and  T.  hildebrandtii  were  consistent 
across  all  three  trees  (Figs.  2-12-2-14).  The  neighbor-joining  and  parsimony  analyses 
also  yielded  a  clade  II  subgroup  containing  T.  thiebautii  strains  VI- 1  and  St.  3-4  (Figs. 
2-12-2-13). 
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96 


Trichodesmium  sp.  H9-4 
Trichodesmium  tenue  tenue 
Trichodesmium  pelagicum  ZK  (AF518769) 
r  Trichodesmium  thiebautii  VI-1 
Trichodesmium  thiebautii  St.  5-1 J W 
Trichodesmium  thiebautii  St.  1-2 


75 


99 


II 


55  (-  Trichodesmium  thiebautii  11-3 
Trichodesmium  thiebautii  St.  3-4 

-  Trichodesmium  hildebrandtii  (AF091322) 

I —  Trichodesmium  erythraeum  21-75 
Q-  Trichodesmium  contortum  20-70 
££  L-  Trichodesmium  contortum  21-74 

-  Trichodesmium  havanum  (AF518770) 
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Figure  2-3:  16S  rDNA  tree  analyzed  by  neighbor-joining.  The  tree  was  bootstrapped 
2000  times  and  Oscillatoria  sancta  PCC  7515  was  used  as  an  outgroup  to  root  the  tree. 
Numbers  at  nodes  reflect  percent  bootstrap  consensus  >  50%.  Boldfaced  text  is  used  to 
indicate  sequences  from  the  Woods  Hole  culture  collection.  Sequences  from  GenBank  are 
indicated  by  accession  numbers  in  parentheses. 
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Figure  2-4:  16S  rDNA  tree  based  on  parsimony.  The  tree  was  bootstrapped  2000  times 
and  Oscillatoria  sancta  PCC  7515  was  used  as  an  outgroup  to  root  the  tree.  Numbers  at 
nodes  reflect  percent  bootstrap  consensus  >  50%.  Boldfaced  text  is  used  to  indicate  se¬ 
quences  from  the  Woods  Hole  culture  collection.  Sequences  from  GenBank  are  indicated 
by  accession  numbers  in  parentheses. 
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Figure  2-5:  16S  rDNA  tree  based  on  maximum  likelihood.  The  tree  was  bootstrapped 
2000  times  and  Oscillatoria  sancta  PCC  7515  was  used  as  an  outgroup  to  root  the  tree. 
Numbers  at  nodes  reflect  percent  bootstrap  consensus  >  50%.  Boldfaced  text  is  used  to 
indicate  sequences  from  the  Woods  Hole  culture  collection.  Sequences  from  GenBank  are 
indicated  by  accession  numbers  in  parentheses. 
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Figure  2-6:  ITS  tree  analyzed  by  neighbor-joining.  The  tree  was  bootstrapped  2000x 
and  Oscillatoria  sancta  PCC  7515  was  used  as  an  outgroup  to  root  the  tree.  Bootstrap 
values  >  50%  are  indicated  at  the  nodes  as  space  allowed.  The  bootstrap  value  for  the 
T.  tenue  tenue  and  Trichodesmium  sp.  H9-4  node  is  91;  the  value  for  the  T.  thiebautii 
VI- 1,  T.  thiebautii  St.  3-4,  and  T.  thiebautii  II-3  node  is  64;  the  value  for  the  T.  thiebautii 
VI-1,  T.  thiebautii  St.  3-4  node  is  67;  the  value  for  the  T.  contortum  21-74  and  T.  con¬ 
tortum  20-70  node  is  77;  and  the  value  for  the  T.  erythraeum  St.  6-1,  T.  erythraeum  St. 

6-2,  T.  erythraeum  St.  6-5,  and  T.  erythraeum  St.  5-1 A  is  99.  Boldfaced  text  is  used  to 
indicate  sequences  from  the  Woods  Hole  culture  collection.  Sequences  from  GenBank  are 
indicated  by  accession  numbers  in  parentheses. 
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Figure  2-7:  ITS  tree  analyzed  by  parsimony  tree.  The  tree  was  bootstrapped  2000  times 
and  rooted  using  Oscillatoria  sancta  PCC  7515.  Bootstrap  values  >  50%  are  indicated 
at  the  nodes  as  space  allowed.  The  bootstrap  value  for  the  T.  thiebautii  VI- 1,  T.  thiebautii 
St.  3-4,  and  T.  thiebautii  II-3  node  is  62  and  the  value  for  the  T.  erythraeum  St.  6-1,  32 
erythraeum  St.  6-2,  T.  erythraeum  St.  6-5,  and  T.  erythraeum  St.  5-1A  is  73.  Boldfaced 
text  is  used  to  indicate  sequences  from  the  Woods  Hole  culture  collection.  Sequences 
from  GenBank  are  indicated  by  accession  numbers  in  parentheses. 
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Figure  2-8:  ITS  tree  analyzed  by  maximum  likelihood.  Oscillatoria  sancta  PCC  7515  was 
used  as  an  outgroup  to  root  the  tree.  Boldfaced  text  is  used  to  indicate  sequences  from 
the  Woods  Hole  culture  collection.  Sequences  from  GenBank  are  indicated  by  accession 
numbers  in  parentheses.  This  tree  was  not  bootstrapped. 
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Trichodesmium  tenue  (AF013033) 
■  Arthrospira  platensis  FIZ01  (EU427544) 


Figure  2-9:  hetR  tree  based  on  neighbor-joining.  The  tree  was  bootstrapped  2000  times 
and  rooted  using  Arthrospira  platensis  HZ01.  Bootstrap  values  >  50%  are  indicated  at 
the  nodes  as  space  allowed.  The  bootstrap  value  for  the  Trichodesmium  tenue  tenue,  Tri¬ 
chodesmium  sp.  H9-4,  and  I.  thiebautii  II-3  node  is  85;  the  value  for  the  Trichodesmium 
sp.  H9-4  and  T.  thiebautii  II-3  node  is  68;  and  the  value  for  the  T.  thiebautii  St.  3-4  and 
VI- 1  node  is  88.  Boldfaced  text  is  used  to  indicate  sequences  from  the  Woods  Hole 
culture  collection.  Sequences  from  GenBank  are  indicated  by  accession  numbers  in 
parentheses. 
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■  Arthrospira  platensis  HZ01  (EU427544) 


Figure  2-10:  hetR  tree  based  on  parsimony.  The  tree  was  bootstrapped  2000  times  and 
rooted  using  Arthrospira  platensis  HZ01.  Bootstrap  values  >  50%  are  indicated  at  the 
nodes  as  space  allowed.  The  bootstrap  value  for  the  T.  thiebautii  St.  3-4  and  VI- 1  node  is 
63  and  the  value  for  the  T.  erythraeum  K-11  #131,  St.  6-2,  St.  6-5,  St.  6-1,  and  St.  5-1A 
node  is  78.  Boldfaced  text  is  used  to  indicate  sequences  from  the  Woods  Hole  culture 
collection.  Sequences  from  GenBank  are  indicated  by  accession  numbers  in  parentheses. 
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Figure  2-11:  hetR  tree  based  on  maximum  likelihood.  The  tree  was  bootstrapped  2000 
times  and  rooted  using  Arthrospira  platensis  HZ01.  Bootstrap  values  >  50%  are  indi¬ 
cated  at  the  nodes  as  space  allowed.  The  bootstrap  value  for  the  T.  thiebautii  St.  3-4  and 
VI-1  node  is  63  and  the  value  for  the  T.  erythraeum  K-11  #131,  St.  6-2,  St.  6-5,  St.  6-1, 
and  St.  5-1 A  node  is  63.  Boldfaced  text  is  used  to  indicate  sequences  from  the  Woods 
Hole  culture  collection.  Sequences  from  GenBank  are  indicated  by  accession  numbers  in 
parentheses. 
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Figure  2-12:  Tree  of  concatenated  sequences  of  16S,  ITS  and  hetR  analyzed  by  neighbor¬ 
joining.  The  tree  was  bootstrapped  2000  times  and  percent  bootstrap  consensus  values  > 
50%  are  indicated  at  the  nodes  as  space  allowed.  The  bootstrap  value  for  the  T.  thiebautii 
VI- 1,  T.  thiebautii  St.  3-4,  and  T.  thiebautii  II-3  node  is  65;  the  value  for  the  T.  thiebautii 
VI- 1  and  T.  thiebautii  St.  3-4  node  is  96;  the  value  for  the  T.  thiebautii  St.  5-1JW  and 
T.  thiebautii  St.  1-2  is  76;  and  the  value  for  the  T.  erythraeum  St.  6-5,  T.  erythraeum  St. 
6-1,  T.  erythraeum  St.  5-1A,  and  T.  erythraeum  St.  6-2  is  96.  Sequences  which  originated 
completely  from  this  study  are  in  boldface.  Sequences  which  contain  sequences  from 
GenBank  are  indicated  by  normal  text. 
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Figure  2-13:  Tree  of  concatenated  sequences  of  16S,  ITS  and  hetR  analyzed  by  parsimony 
tree.  The  tree  was  bootstrapped  2000  times  and  percent  bootstrap  consensus  values  > 

50%  are  indicated  at  the  nodes.  Sequences  which  originated  completely  from  this  study 
are  in  boldface.  Sequences  which  contain  sequences  from  GenBank  are  indicated  by 
normal  text. 
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Trichodesmium  erythraeum  St.  6-5 
|-  Trichodesmium  erythraeum  St.  6-1 
Trichodesmium  erythraeum  St.  5-1 A 
Trichodesmium  erythraeum  St.  6-2 
L  Trichodesmium  erythraeum  K  11-131 
Trichodesmium  erythraeum  XIII-13 
Trichodesmium  erythraeum  K-04  #20 
|  Trichodesmium  contortum  20-70 
'  Trichodesmium  contortum  21-74 
Trichodesmium  erythraeum  IMS101 
Trichodesmium  erythraeum  GBRTRLIN201 
Trichodesmium  erythraeum  K-02  #2 
I-  Trichodesmium  erythraeum  GBRTRLI101 
' —  Trichodesmium  erythraeum  21-75 

y 

0.01  substitutions/site 
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gg  |  Trichodesmium  sp.  H9-4 

Trichodesmium  tenue  tenue 
h  Trichodesmium  thiebautii  St.  3-4 
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-  Trichodesmium  thiebautii  11-3 

-  Trichodesmium  thiebautii  VI-1 
—  Trichodesmium  pelagicum  (20-71  and  ZK) 

Trichodesmium  thiebautii  St.  5-1 J W 
H"  Trichodesmium  thiebautii  St.  1  -2 

Trichodesmium  hildebrandtii  (11-4  and  field) 


II 


I 


Figure  2-14:  Tree  of  concatenated  sequences  of  16S,  ITS  and  hetR  analyzed  by  maximum 
likelihood.  The  tree  was  bootstrapped  2000  times  and  percent  bootstrap  consensus  values 
>  50%  are  indicated  at  the  nodes  as  space  allowed.  The  bootstrap  value  for  the  T.  ery¬ 
thraeum  St.  6-5,  T.  erythraeum  St.  6-1,  T.  erythraeum  St.  5-1A,  and  T.  erythraeum  St.  6-2 
node  is  70  and  the  value  for  the  T.  contortum  20-70  and  T.  contortum  21-74  node  is  76. 
Sequences  which  originated  completely  from  this  study  are  in  boldface.  Sequences  which 
contain  sequences  from  GenBank  are  indicated  by  normal  text. 
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Cell  morphology.  Both  clades  of  Trichodesmium  exhibited  wide  diversity  in  cell  mor¬ 
phology.  In  clade  I,  the  cell  shape  of  all  strains  of  T.  erythraeum  ranged  from  square  to 
wider  than  long  (Table  2.4,  Fig.  2-15).  Widths  ranged  from  5.3  -  13.6  pm  while  lengths 
ranged  from  3.9  -  12.2  p m.  Clade  I  T.  contortum  cells  were  all  disc-shaped  with  widths 
ranging  from  19.5  -  34.5  p m  and  lengths  ranging  from  2.0  -  4.3  pm  (Table  2.4,  Fig.  2-15). 
Clade  II  cell  shapes  ranged  from  being  slightly  wider  than  long  to  being  longer  than  wide 
(Table  2.4,  Fig.  2-16).  Cell  widths  in  T.  thiebautii  strains  ranged  from  6.0  -  11.8  pm  and 
lengths  ranged  from  5.6  -  22.4  p m.  T.  tenue  tenue  was  longer  than  wide  (7.9  -  10.3  p m 
compared  to  5.0  -  7.9  p m,  respectively),  while  Trichodesmium  sp.  H9-4  was  shorter  than 
wide  (4.4  -  8.6  p m  compared  to  5.7  -  8.9  pm,  respectively).  T.  pelagicum  20-71  was  the 
largest  of  the  clade  at  12.1  -  18.4  pm  wide  and  5.2  -  9.6  pm  long  (Table  2.4,  Fig.  2-16). 
Gas  vesicles  can  be  seen  as  bright  regions  under  differential  interference  contrast  (DIC) 
(Figs.  2-15,  2-16).  Gas  vesicle  distribution  within  the  clade  I  varied  from  random  dis¬ 
tribution  in  strains  IMS101,  20-70,  and  21-74  to  peripheral  distribution  in  all  the  other 
strains  in  the  clade.  Clade  II  gas  vesicle  distribution  ranged  from  random  distribution  in 
strains  St.  1-2,  St.  5-1JW,  tenue,  and  20-71  to  peripheral  distribution  in  strains  II-3,  VI-1, 
St.  3-4,  and  H9-4  (Fig.  2-16). 
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Table  2.4:  Summary  of  the  physical  characteristics  of  Trichodesmium  cultures.  Numbers  in  parentheses  are  averages  over  > 
100  cells.  PUB  =  phycourobilin,  PEB  =  phycoerythrobilin,  PE  =  phycoerythrin,  and  PC  =  phycocyanin. 

Species  Strain  Cell  Size  (/jm)  Absorption  peaks  (nm)  PUB:PEB  PE:PC 
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H9-4  5.7  -  8.9  (7.4)  4.4  -  8.6  (6.3)  493.0  555.0  615.0  1.54  3.68 

T.  pelagicum  20-71  12.1  -  18.4  (14.6)  5.2  -  9.6  (7.3) 

(. Katagnymene ) 
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Figure  2-15:  Trichodesmium  clade  I  micrographs  taken  using  phase  contrast  and  differ¬ 
ential  interference  contrast  (DIC):  T.  erythraeum  (A  -  V)  and  T.  contortum  (W  -  Z).  The 
white  bar  in  (A)  represents  20  /urn.  (A)  IMS  101,  phase  contrast;  (B)  IMS  101,  DIC;  (C) 
XIII-13,  phase  contrast;  (D)  XIII-13,  DIC;  (E)  St.  5-1A,  phase  contrast;  (F)  St.  5-1A, 
DIC;  (G)  St.  6-1,  phase  contrast;  (H)  St.  6-1,  DIC;  (I)  St.  6-5,  phase  contrast;  (J)  St.  6-5, 
phase  contrast;  (K)  K-02  #2,  phase  contrast;  (L)  K-02  #2,  DIC;  (M)  K-04  #20,  phase 
contrast;  (N)  K-04  #20,  DIC;  (O)  K-ll  #131,  phase  contrast;  (P)  K-ll  #131,  DIC;  (Q) 
GBRTRLI101,  phase  contrast;  (R)  GBRTRLI101,  DIC;  (S)  GBRTRLIN201,  phase  con¬ 
trast;  (T)  GBRTRLIN201,  DIC;  (U)  21-75,  phase  contrast;  (V)  21-75,  DIC;  (W)  20-70, 
phase  contrast;  (X)  20-70,  DIC;  (Y)  21-74,  phase  contrast;  (Z)  21-74,  DIC. 
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Figure  2-16:  Trichodesmium  clade  II  micrographs  taken  using  phase  contrast  and  differ¬ 
ential  interference  contrast  (DIC):  T.  thiebautii  (A  -  J),  T.  tenue  (K  -  N),  and  T.  pelagicum 
(O  -  P).  The  white  bar  in  (A)  represents  20  pm.  (A)  II-3,  phase  contrast;  (B)  II-3,  DIC; 

(C)  VI-1,  phase  contrast;  (D)  VI-1,  DIC;  (E)  St.  1-2,  phase  contrast;  (F)  St.  1-2,  DIC; 

(G)  St.  3-4,  phase  contrast;  (H)  St.  3-4,  DIC;  (I)  St.  5-1JW,  phase  contrast;  (J)  St.  5-1JW, 
DIC;  (K)  tenue,  phase  contrast;  (L)  tenue,  DIC;  (M)  H9-4,  phase  contrast;  (N)  H9-4,  DIC; 
(O)  20-71,  phase  contrast;  (P)  20-71,  DIC. 
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Absorption  spectra  of  phy cobilip roteins.  The  major  phycobiliprotein  in  all  the  mea¬ 
sured  strains  of  Trichodesmium  was  PE  with  the  ratio  of  PUB:PEB  varying  between  species 
clusters.  All  strains  had  a  PE:PC  ratio  >  1  and  the  AP  peak  was  not  detectable  (Table  2.4, 
Fig.  2-17). 

In  clade  I,  all  strains  of  T.  erythraeum  and  T.  contortum  had  a  PUB:PEB  ratio  < 

1  (0.79  -  0.92).  PE:PC  ratios  ranged  from  1.75  -  3.64  for  T.  erythraeum  and  5.11  -  5.19 
for  T.  contortum  (Table  2.4).  In  this  clade,  the  wavelength  of  maximum  absortion  (kmax) 
for  PUB  was  497.0  -  498.0  nm,  Xmax  for  PEB  was  543.0  -  546.5  nm,  and  Xmcx  for  PC  was 
613.5  -  620.0  nm  (Table  2.4,  Fig.  2-17  A).  Some  strains  in  this  clade  showed  a  shoulder 
on  the  PEB  peak  around  555  nm:  T.  erythraeum  strains  IMS  101,  K-02  #2,  K-04  #20, 

XIII- 13,  St.  5-1A,  GBRTRLI101,  GBRTRLIN201,  and  21-75  and  T.  contortum  20-70 
(Fig  2-17  A).  T.  erythraeum  strains  GBRTRLI101  and  GBRTRLIN201  had  absorption 
peaks  at  436.0  and  671.0  nm,  indicating  small  amounts  of  chlorophyll  contamination  (Fig. 
2-17  A). 

In  clade  II,  all  strains  of  T.  thiebautii  and  T.  tenue  had  a  PUB:PEB  ratio  >  1  (1.15 
-  1.54),  PE:PC  ratios  ranged  from  2.25  -  5.54,  and  T.  tenue  tenue  did  not  have  a  discernible 
PC  peak  (Table  2.4,  Fig.  2-17  B).  In  this  clade,  rkmax  for  PUB  ranged  from  492.5  -  496.0 
nm,  A .max  for  PEB  ranged  from  540.0  -  555.0  nm,  and  \max  for  PC  ranged  from  615.0  - 
620.0  nm  (Table  2.4,  Fig.  2-17  B). 

Discussion 

The  phylogenetic  trees  show  two  distinct  clades  of  cultured  Trichodesmium ,  and  within 
each  clade  are  consistent  subgroups.  In  clade  I,  T.  erythraeum  strains  St.  5-1A,  St.  6-1,  St. 
6-2,  and  St.  6-5  group  together  and  T.  contortum  strains  20-70  and  21-74  group  together. 

In  clade  II,  strains  tenue  and  H9-4  cluster  together  and  T.  thiebautii  strains  St.  1-2  and  St. 
5-1JW  cluster  with  T.  hildebrandtii.  In  general,  the  subgroups  are  morphologically  similar 
to  each  other  in  size  and  gas  vesicle  distribution.  The  exception  to  this  is  the  tenue/H9-4 
subgroup. 
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Wavelength  (nm) 


Figure  2-17:  Absorption  spectra  of  phycobiliproteins  from  cultured  strains  of  Tri- 
chodesmium.  Absorbance  is  normalized  to  the  maximum  absorbance.  Graphs  are  offset 
vertically  to  show  details.  (A)  Clade  I  spectra  with  PUB:PEB  <  1.  Solid  lines  indicate  T. 
erythraeum  and  dashed  lines  indicate  T.  contortum.  (B)  Clade  II  spectra  with  PUB:PEB 
>  1.  Solid  lines  indicate  T.  thiebautii  and  dashed  lines  indicate  T.  tenue.  Vertical  lines 
denote  typical  wavelengths  of  absorption  peaks  for  PUB  (498  nm),  PEB  (545  nm),  and 
PC  (618  nm). 


Genotypes  of  Trichodesmium  are  circumtropical.  The  North  Atlantic  strain  T.  ery¬ 
thraeum  IMS  101  is  more  similar  to  strains  from  the  South  Pacific  such  as  GBRTRLIN201 
and  K-02  #2  than  it  is  to  other  North  Atlantic  strains  such  as  St.  6-1.  Likewise,  the  Indian 
Ocean  strain  T.  tenue  tenue  is  most  similar  to  the  North  Pacific  strain  H9-4. 

The  genetic  grouping  of  species  differs  greatly  from  classical  groupings  based  on 
morphology,  in  which  T.  tenue  was  clustered  with  T.  erythraeum,  and  T.  contortum  was 
clustered  with  T.  thiebautii  (Janson  et  al.,  1995).  Our  culture  collection  does  not  encom- 
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pass  the  full  range  of  Trichodesmium  diversity.  Genotypes  identified  as  T.  contortion  and 
T.  tenue  in  Janson  et  al.  (1999)  and  T.  aureum  and  T.  havanum  (Lundgren  et  al.,  2005)  are 
not  represented  among  our  cultures.  The  Janson  T.  contortum  and  T.  tenue  formed  a  clade 
separate  from  our  cultured  strains.  The  contradiction  in  genotypes  of  T.  contortum  and  T. 
tenue  found  between  our  study  and  Janson  et  al.  (1999)  may  be  due  to  different  species 
being  identified  with  the  same  name.  T.  contortum  specimens  in  Janson  et  al.  (1999)  oc¬ 
curred  in  colonies  and  were  “a  pale,  straw-like  color.”  Ours  was  bright  red,  indicating 
high  PEB  content,  and  occurred  as  single  filaments.  The  T.  tenue  colonies  in  Janson  et  al. 
(1999)  were  larger  and  darker  red  than  usually  seen  in  T.  tenue  and  had  peripheral  gas 
vesicles.  Our  T.  tenue  tenue  has  dispersed  gas  vesicles.  Culturing  efforts  must  continue  to 
incorporate  this  diversity,  and  classification  systems  must  be  updated  and  revised  as  new 
data  is  available  (Anagnostidis  and  Komarek,  1988). 

Morphological  characteristics  can  aid  in  the  identification  of  Trichodesmium  spp., 
but  they  are  not  infallible.  Gas  vesicles  are  variable  and  may  not  be  a  good  characteris¬ 
tic  for  identification.  Cells  of  the  same  species  may  display  both  random  distribution  and 
peripheral  arrangement  of  gas  vesicles.  Cell  morphometries  have  a  broad  range  within 
a  clade  and  overlap  between  clades.  Even  IMS  101,  the  most  widely-used  cultural  repre¬ 
sentative  of  T.  erythraeum,  was  initially  identified  as  T.  thiebautii  (Prufert-Bebout  et  al., 
1993).  Although  T.  contortum  is  genetically  very  similar  to  T.  erythraeum,  they  look 
vastly  different.  Similarly,  T.  pelagicum  is  so  morphologically  different  from  T.  thiebautii 
that  until  genetic  work  was  introduced,  it  was  considered  a  separate  genus,  Katagnymene 
(Lundgren  et  al.,  2001;  Orcutt  et  al.,  2002;  Lundgren  et  al.,  2005). 

Phycobiliprotein  composition  was  consistent  with  the  phylogeny  of  the  two  ma¬ 
jor  clades  of  Trichodesmium.  This  result  is  unexpected  since  phycobiliproteins  have  not 
been  reliable  taxonomic  markers.  For  example,  PUB:PEB  ratios  in  closely-related  strains 
of  Synechococcus  ranged  from  0.4  -  2.0  while  unrelated  strains  had  similar  PUB:PEB 
ratios  (Toledo  et  al.,  1999).  Trichodesmium  strains  in  clade  I  were  all  brick  red  in  color 
due  to  the  prevalence  of  PEB  (PUB:PEB  =  0.79  -  0.92).  The  spectra  in  this  clade  were 
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highly  similar.  Some  of  the  strains  displayed  a  shoulder  on  the  red  side  of  the  PEB  peak. 

The  two  T.  contortum  strains  were  similar  in  that  they  had  much  less  PC  than  the  T.  ery- 
thraeum  strains,  shown  by  the  PE:PC  ratio:  5.11  -  5.19  compared  to  1.75  -  3.64  in  T. 
erythraeum.  The  two  T.  contortum  strains  differed  in  that  20-70  had  the  PEB  shoulder 
while  21-74  did  not.  Clade  II  varied  in  color  and  Xmax  of  the  PE  subunits,  but  they  all  had 
PUB  as  the  dominant  subunit  of  PE.  With  increasing  PUB:PEB,  cells  become  greener: 

T.  tenue  tenue  is  salmon  pink  (PUB:PEB  =  1.15),  T.  thiebautii  strains  tend  to  be  brown 
(PUB:PEB  =  1.15  -  1.39),  and  Trichodesmium  sp.  H9-4  is  sage  green  (PUB:PEB  =  1.54). 
The  Xmax  for  PUB  (493.0  -  496.0  nm)  is  blue-shifted  for  this  clade  compared  to  that  of 
clade  I  (497.0  -  498.0  nm).  In  clade  II,  the  'kmax  for  PEB  was  red-shifted  compared  to  that 
of  clade  I  (548.0  -  549.0  nm  compared  to  543.0  -  546.5  nm,  respectively).  Trichodesmium 
sp.  H9-4  had  the  most  red-shifted  'kmax  for  PEB  at  555  nm  and  T.  tenue  tenue  had  the 
most  blue-shifted  Tmax  for  PEB  at  540.0  nm,  the  shortest  wavelength  for  all  strains  re¬ 
gardless  of  species. 

The  consistency  of  absorption  spectra  within  monophyletic  groups  of  Trichodesmium 
indicated  that  pigment  composition  may  be  used  to  aid  in  species  identification  where  cell 
morphology  may  be  ambiguous.  In  South  Pacific  field  populations  of  Trichodesmium , 
the  spectra  of  field  T.  erythraeum  matched  the  absorption  spectra  presented  here  for  clade 
I,  and  the  spectrum  of  a  T.  lhiehautii-\\kc  green  colony  resembled  the  spectra  for  clade 
II  (Neveux  et  al.,  2006).  However,  South  Pacific  T.  thiebautii  and  T.  pelagicum  ( Katag - 
nymene)  samples  had  lower  PUB:PEB  ratios  than  clade  II  cultured  samples:  0.98  -  1.00 
and  0.68,  respectively,  compared  to  1.15  -  1.54  (Neveux  et  al.,  2006).  A  Caribbean  Sea 
population  identified  as  a  mixture  of  T.  erythraeum  and  T.  thiebautii  had  absorption  spec¬ 
tra  similar  to  that  of  clade  II  (Subramaniam  et  al.,  1999).  This  study  indicates  that  Tri¬ 
chodesmium  with  low  PUB:PEB  ratios  belong  to  clade  I  while  those  with  high  PUB:PEB 
ratios  belong  to  clade  II,  but  further  studies  are  required  to  understand  the  phenotypic 
plasticity  within  each  clade. 

Diverse  PUB:PEB  ratios  point  to  the  possibility  of  niche  differentiation  in  Tri- 
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chodesmium.  In  the  South  Pacific,  green  colonies  of  Trichodesmium  with  high  PUB  were 
found  at  50  -  120m,  deeper  than  the  high  PEB  Trichodesmium  collected  (Neveux  et  al., 
2006).  PUB  absorbs  shorter  wavelengths  of  light  than  PEB  and  hence  is  adapted  for  greater 
depths  (Neveux  et  al.,  2006).  Conversely,  Subramaniam  et  al.  (1999)  found  that  PUB:PEB 
increased  from  1.7  at  night  to  2.1  at  midmoming  for  Caribbean  Trichodesmium  specimens 
at  15  m.  The  increase  in  PUB  may  downregulate  the  energy  supply  to  photosystem  II  to 
prevent  damage  to  the  reaction  center  at  high  light  (Subramaniam  et  al.,  1999).  Interpreta¬ 
tion  of  the  role  of  PUB  in  Trichodesmium  ecology  requires  environmental  and  physiologi¬ 
cal  context. 

The  green  field  colonies  described  in  Neveux  et  al.  (2006)  may  be  closely  related 
to  our  cultured  strain  H9-4.  Orcutt  et  al.  (2008)  reported  that  the  DNA  fingerprint  of  the 
green  colonies  was  identical  to  that  of  T.  thiebautii.  While  the  green  color  has  been  sug¬ 
gested  to  be  an  indication  of  senescence  (Orcutt  et  al.,  2008),  the  green  field  colonies 
were  viable  (Neveux  and  Tenorio,  2008),  and  green  appears  to  be  the  native  color  of 
H9-4. 

H9-4  is  a  taxonomic  enigma.  Its  cell  morphology  is  most  similar  to  that  of  T. 
thiebautii;  its  genetic  markers  most  resemble  T.  tenue;  it  forms  single  trichomes,  puffs, 
and  tufts  all  in  the  same  flask;  and  its  color  is  like  nothing  else:  green  with  a  long  Xmax 
for  PEB  (555.0  nm).  The  morphological  characteristics  of  H9-4  and  T.  tenue  tenue  are  di¬ 
vergent.  H9-4  usually  has  peripheral  gas  vesicles  while  tenue’s  are  randomly  distributed, 
H9-4  cells  are  mostly  square-shaped  while  tenue  cells  are  much  longer  than  they  are  wide, 
H9-4  cells  are  green  while  tenue  cells  are  pink,  and  H9-4  has  different  PUB  and  PEB 
Xmax  s  than  tenue:  493.0  nm  compared  to  496.0  nm  and  555.0  nm  compared  to  540.0  nm, 
respectively.  In  spite  of  the  morphological  and  pigmentation  differences,  I  have  classified 
H9-4  as  T.  tenue  based  on  its  genetic  similarity  to  T.  tenue  tenue. 

Using  sequence  data,  I  have  been  able  to  show  that  there  are  two  distinct  clades 
of  Trichodesmium  represented  in  our  culture  collection.  Studying  the  metagenome  of 
these  two  clades  could  help  us  to  understand  potentially  important  differences  between 
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them  and  to  fine-tune  the  phylogenetic  relationships  of  Trichodesmium.  Finding  a  ge¬ 
netic  basis  for  cell  size  and  shape  could  help  to  shed  light  on  the  morphological  diversity 
found  within  each  clade.  Future  studies  examining  the  temperature  optima  and  chromatic 
adaptions  of  cultured  strains  could  illuminate  the  possibilities  for  niche  differentiation 
in  Trichodesmium.  The  culture  collection  does  not  incorporate  the  full  diversity  of  Tri¬ 
chodesmium  found  in  field  studies,  so  culturing  efforts  should  continue. 
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CHAPTER  3 


Development  of  a  quantitative  polymerase  chain 
reaction  (qPCR)  assay  to  distinguish  the  two 
clades  of  Trichodesmium 


ABSTRACT:  To  understand  the  role  Trichodesmium  has  in  the  global  nitrogen  cycle,  we 
must  understand  the  distribution  of  Trichodesmium  spp.  Previous  applications  of  real-time 
quantitative  polymerase  chain  reaction  (qPCR)  to  Trichodesmium  populations  have  been 
at  the  genus  level.  Investigating  ecological  differences  among  species  of  Trichodesmium 
requires  an  assay  with  higher  resolution,  but  designing  a  species-specific  assay  is  diffi¬ 
cult  due  to  the  high  degree  of  genetic  similarity  among  Trichodesmium  spp.  There  are  two 
distinct  clades  of  Trichodesmium ;  clade  I  contains  Trichodesmium  erythraeum  and  Tri¬ 
chodesmium  contortum,  and  clade  II  contains  Trichodesmium  thiebautii,  Trichodesmium 
tenue,  Trichodesmium  hildebrandtii,  and  Trichodesmium  pelagicum.  In  this  study,  I  de¬ 
veloped  a  qPCR  method  that  is  able  to  distinguish  between  the  two  major  clades  of  Tri¬ 
chodesmium ,  targeting  the  heterocyst  differentiation  gene  hetR.  Trichodesmium  spp.  pro¬ 
duce  compounds  which  can  inhibit  PCR,  making  quantitative  assays  troublesome,  so  I 
tested  a  variety  of  DNA  extraction  methods  to  minimize  inhibition.  This  method  was 
tested  in  pure  and  mixed  cultures  of  Trichodesmium  for  specificity  and  efficiency,  and 
can  be  applied  to  field  populations  to  identify  and  enumerate  the  two  clades. 
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Introduction 


Cyanobacteria  of  the  genus  Trichodesmium  contribute  significant  amounts  of  nitrogen  (N) 
to  tropical  and  subtropical  oligotrophic  oceans  through  N2  fixation  (Capone  et  al.,  1997). 
In  order  to  constrain  inputs  of  N  to  marine  ecosystems,  the  spatial  and  temporal  distri¬ 
butions  of  Trichodesmium  spp.  must  be  understood.  Cell  counts,  specifically  trichome 
(filament)  and  colony  counts,  can  be  effective  (Letelier  and  Karl,  1996;  Post  et  al.,  2002; 
Tyrrell  et  al.,  2003).  In  the  Gulf  of  Aqaba,  Red  Sea,  fusiform  “tuft”  colonies  were  found 
near  the  surface  while  spherical  “puff”  colonies  were  found  in  the  lower  half  of  the  photic 
zone  (Post  et  al.,  2002).  In  the  eastern  Atlantic,  Trichodesmium  was  found  at  highest  tri¬ 
chome  and  colony  densities  between  0  —  15°N,  with  a  complete  absence  south  of  30°S 
(Tyrrell  et  al.,  2003).  Trichomes  and  colonies  vary  in  size  and  their  counts  may  not  rep¬ 
resent  true  Trichodesmium  biomass.  Cell  counts  can  also  be  laborious,  and  identifying 
trichomes  to  species  based  on  morphology  can  be  confounding  (Chapt.  2). 

Trichodesmium  biomass  was  long-considered  to  be  concentrated  in  the  upper  50 
m  of  the  water  column  (Letelier  and  Karl,  1996;  Carpenter  et  al.,  2004).  Collection  of  Tri¬ 
chodesmium  via  net  tows  is  limited  to  surface  regions,  as  long  hauls  can  damage  the  frag¬ 
ile  trichomes.  A  non-invasive  study  of  Trichodesmium  distributions  using  a  video  plank¬ 
ton  recorder  (VPR)  reported  high  concentrations  of  Trichodesmium  colonies  as  deep  as 
130  m  (Davis  and  McGillicuddy,  2006).  The  VPR  also  revealed  disparate  depth  distribu¬ 
tions  between  colony  morphologies;  tufts  were  concentrated  near  the  surface  while  puffs 
were  distributed  throughout  the  upper  water  column  (Davis  and  McGillicuddy,  2006). 

The  VPR  provides  an  important  alternative  to  net  tows  and  cell  counts.  However,  the  VPR 
is  unable  to  distinguish  species  of  Trichodesmium  and  does  not  detect  single  trichomes, 
which  can  make  up  the  majority  of  biomass  (Letelier  and  Karl,  1996). 

Real-time  quantitative  polymerase  chain  reaction  (qPCR)  has  been  used  to  inves¬ 
tigate  population  dynamics  of  marine  cyanobacteria  on  spatial  and  temporal  scales.  Eco¬ 
types  of  Prochlorococcus  have  been  resolved  using  sequences  of  the  16S-23S  internal 
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transcribed  spacer  (ITS)  (Ahlgren  et  al.,  2006;  Zinser  et  al.,  2006;  Johnson  et  al.,  2006). 
Diazotrophs  including  heterocystous  Richelia  spp.,  unicellular  Crocosphaera  spp,  uni¬ 
cellular  “group  A,”  and  non-heterocystous  Trichodesmium  spp.  have  been  enumerated 
by  targeting  nifH,  one  of  the  genes  encoding  the  N2 -fixing  enzyme  nitrogenase  (Church 
et  al.,  2005;  Foster  et  al.,  2007;  Goebel  et  al.,  2007). 

The  nifH  qPCR  assay  targets  the  genus  Trichodesmium  as  a  whole.  Species  of 
Trichodesmium  are  highly  similar  genetically,  making  species-specific  primers  difficult  to 
design.  The  species  of  Trichodesmium  fall  into  two  clades.  Clade  I  includes  Trichodesmium 
erythraeum  and  Trichodesmium  contortum,  and  clade  II  includes  Trichodesmium  thiebau- 
tii,  Trichodesmium  tenue,  Trichodesmium  pelagicum  (formerly  Katagnymene  spp.),  and 
Trichodesmium  hildebrandtii  (Chapt.  2).  In  order  to  elucidate  important  ecological  differ¬ 
ences  between  the  clades,  a  qPCR  assay  must  target  each  clade  individually. 

The  purpose  of  this  study  was  to  develop  a  SYBR  green  qPCR  method  to  enu¬ 
merate  and  identify  the  two  major  clades  of  Trichodesmium.  SYBR  green  fluoresces  only 
when  bound  to  double-stranded  DNA.  As  PCR  progresses,  the  amplicon  increases  ex¬ 
ponentially,  and  the  SYBR  green  signal  increases.  Artifacts  such  as  primer-dimers,  sec¬ 
ondary  structures,  and  non-specific  products  can  increase  the  SYBR  green  signal  and 
overestimate  the  quantity  of  the  region  of  interest.  Melting  temperature  (T,„)  of  a  PCR 
product  is  affected  by  its  length  as  well  as  its  sequence,  so  non-specific  products  can  be 
detected  using  a  melt  curve  analysis.  For  an  effective  assay,  amplicons  should  be  <  200 
bp,  primers  should  not  bind  to  themselves  or  each  other  to  form  dimers,  the  amplicon 
should  not  have  secondary  structures  such  as  hairpins,  and  the  primers  should  not  am¬ 
plify  non-specific  products  (Templeton  and  Claas,  2003).  I  used  sequences  of  the  hete¬ 
rocyst  differentiation  protein  gene  hetR  from  Chapt.  2  to  design  primers  targeting  each 
clade  of  Trichodesmium.  I  tested  the  qPCR  primers  with  serial  dilutions  and  mixtures  of 
Trichodesmium  cultures  using  a  variety  of  DNA  extraction  methods  to  minimize  PCR 
inhibition.  The  method  is  able  to  accurately  enumerate  and  identify  the  two  clades  of  Tri¬ 
chodesmium  in  mixed  cultures. 
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Table  3.1:  Sequences  of  clade-specific  primers  used  in  qPCR 


Clade 

Primer 

Sequence  (5’ to  3’) 

I 

hetR-TeryF 

hetR-TeryR 

GCAGGAAGTAGAAGAAGCAC 

TTACCAGGAAGGCGAGAAC 

II 

hetR-ThiebF 

hetR-ThiebR 

TTCCTGGTGGTAATCTGACTAC 

CTTCCTGAGAACGCCTATGG 

Methods 

Primer  design  and  preparation  of  cloned  qPCR  standards.  Sequences  of  the  hete¬ 
rocyst  differentiation  gene  hetR  from  cultured  strains  of  Trichodesmium  (Chapt.  2)  were 
used  to  design  primers  for  SYBR  green  using  AllelelD  software  (PREMIER  Biosoft  In¬ 
ternational).  Primers  targeted  different  hetR  sequences  within  each  clade  of  Trichodesmium 
while  excluding  the  other  clade  as  well  as  marine  filamentous  cyanobacteria  Richelia 
(GenBank  accession  numbers  DQ366915  and  DQ366914)  and  Calothrix  (GenBank  ac¬ 
cession  number  DQ366916)  (Benson  et  al.,  2005;  Foster  and  Zehr,  2006).  Forward  and 
reverse  primers  for  each  clade  are  listed  in  Table  3.1  and  amplify  different  but  overlapping 
regions  of  hetR  for  each  clade.  AllelelD  was  also  used  to  design  SYBR  green  primers  tar¬ 
geting  ITS.  The  primer  pairs  were  tested  in  silico  using  Mac  Vector  on  sequences  from  all 
six  species  of  Trichodesmium. 

All  PCR  optimization  was  performed  on  an  iCycler  thermal  cycler  (BioRad)  with 
25  /uL  reactions.  The  primers  were  optimized  using  a  temperature  gradient  to  find  an  an¬ 
nealing  temperature  that  was  low  enough  for  efficient  PCR  in  the  targeted  clade  but  high 
enough  to  exclude  the  non-targeted  clade.  PCR  was  conducted  using  the  Mas  ter  Tag  Kit 
(Eppendorf,  5  Prime)  to  get  a  140  bp  fragment  for  Clade  I  and  a  138  bp  fragment  for 
Clade  II.  Reactions  consisted  of  1  x  Taq  Master  PCR  Enhancer,  1  x  Taq  Buffer,  200  //M 
dNTPs,  0.5  /jM  each  of  primers  hetR-TeryF  and  hetR-TeryR  for  clade  I  and  primers  hetR- 
ThiebF  and  hetR-ThiebR  for  clade  II  (Table  3.1),  1.25  U  Taq  DNA  polymerase,  and  4  ng 
genomic  DNA  from  cultured  strains.  Cycling  conditions  were  94°C  for  2  min;  30  cycles 
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of  94°C  for  30  s,  62°C  (clade  I)  or  61°C  (clade  II)  for  30  s,  and  72°C  for  30  s;  and  one 
cycle  of  72°C  for  5  min.  Both  primer  sets  were  tested  using  DNA  from  T.  erythraeum 
strains  IMS101,  St.  5-1A,  St.  6-1,  XII-13,  and  K-04  #20;  T.  contortion  strains  21-74  and 
20-70;  T.  thiebautii  strains  II-3,  VI-1,  St.  1-2,  St.  3-4,  and  St.  5-1JW;  T.  tenue  strains 
tenue  and  H9-4;  and  T.  pelagicum  20-71  (Chapt.  2).  Five  pL  of  each  PCR  product  was 
separated  on  a  2%  agarose  gel,  stained  with  ethidium  bromide,  and  photographed  using  a 
Gel  Logic  200  Imaging  System  and  Molecular  Imaging  Software  (Kodak).  The  clade  I  as¬ 
say  was  optimized  to  amplify  only  strains  of  clade  I,  and  the  clade  II  assay  was  optimized 
to  amplify  only  strains  of  clade  II. 

Standards  for  qPCR  were  derived  from  cloned  PCR  products  for  each  clade.  To 
provide  PCR  products  for  cloning,  T.  erythraeum  IMS  101  was  amplified  using  the  clade 
I  protocol  and  T.  thiebautii  II-3  was  amplified  using  the  clade  II  protocol  as  described 
above  with  the  modification  of  25  cycles  instead  of  30.  Products  were  cloned  and  lin¬ 
earized  as  described  in  Zinser  et  al.  (2006).  Briefly,  the  products  were  cloned  directly 
into  the  pCR4-TOPO  vector  and  introduced  into  Escherichia  coli  TOP10  cells  using  the 
TOPO-TA  cloning  kit  for  sequencing  (Invitrogen)  following  the  manufacturer’s  direc¬ 
tions.  Clones  were  screened  using  LBKan50  (50  pg  mL  1  kanamycin  in  Luria-Bertani 
medium).  Plasmids  from  the  E.  coli  cells  were  purified  using  the  QIAprep  Spin  Miniprep 
kit  (Qiagen)  following  manufacturer’s  instructions  and  stored  at  -80°C.  To  linearize  the 
purified  plasmids,  digestion  was  performed  using  a  PstI  restriction  enzyme  kit  (New  Eng¬ 
land  Biolabs)  in  30  pL  reactions  with  1  x  buffer  3,  1  x  bovine  serum  albumin  (BSA),  2 
pL  PstI,  and  3  pL  plasmid.  Reactions  were  incubated  at  37°C  for  3  h  and  products  were 
visualized  on  a  1%  agarose  gel  to  confirm  a  single  restriction  site.  Linearized  plasmids 
were  either  sequenced  as  described  in  Chapt.  2  by  the  Keck  Facility  at  the  Josephine  Bay 
Paul  Center,  Marine  Biological  Laboratory,  Woods  Hole,  MA,  or  purified  plasmids  were 
sent  to  Agencourt  Bioscience  Corporation,  Beverly,  MA.  Plasmids  whose  sequences  were 
identical  to  those  of  T.  erythraeum  IMS  101  and  T.  thiebautii  II-3  were  serially  diluted 
(10-fold),  quantitated  using  the  Quant-iT  PicoGreen  kit  (Molecular  Probe s/Invitrogen) 
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following  manufacturer’s  instructions,  and  stored  to  be  used  as  qPCR  standards.  Copies  of 
hetR  volume”1  for  each  standard  was  calculated  using  the  molecular  weight  of  the  plas¬ 
mid  plus  the  insert.  All  products  were  stored  at  -80°C. 

Serial  dilutions  and  cell  counts  of  Trichodesmium  cultures.  Cultures  of  T.  erythraeum 
IMS  101,  T.  contortum  21-74,  T.  thiebautii  II-3,  T.  tenue  tenue,  and  T.  tenue  H9-4  were 
grown  as  described  in  Chapt.  2.  For  dilution  series,  cultures  were  diluted  10-fold  up  to 
10, 000- fold.  For  qPCR,  5-100  mL  of  diluted  cultures  were  vacuum-filtered  on  47  mm, 

5  gm  polycarbonate  filters  (Millipore)  using  a  glass  funnel  and  fritted  base  (Millipore) 
and  filters  were  stored  in  liquid  N2  in  1.8  mL  cryo  tubes  (Nunc).  Dilutions  of  T.  ery¬ 
thraeum  IMS  101  were  mixed  with  dilutions  of  T.  tenue  tenue  and  dilutions  of  T.  con¬ 
tortum  were  mixed  with  dilutions  of  T.  thiebautii  II-3  to  test  the  ability  of  the  qPCR  as¬ 
say  to  distinguish  between  the  two  clades.  For  cell  counts,  1  -  300  mL  diluted  cultures 
were  vacuum-filtered  on  25  mm,  5  gm  polycarbonate  filters  (Millipore)  using  a  glass  fun¬ 
nel  and  fritted  base  (Millipore)  and  fixed  with  1  mL  10%  buffered  formalin  (Surgipath) 
for  2  min.  Fixed  filters  were  mounted  on  a  glass  microscope  slide  using  immersion  oil, 
covered  with  a  25  mm  square  glass  cover  slip,  and  frozen  at  -20°C  until  they  could  be 
viewed.  Slides  were  viewed  using  epifluorescence  (rhodamine  filter)  with  a  Zeiss  Ax- 
ioplan  2  microscope  under  either  a  5x  or  a  lOx  Plan-Neoflaur  objective  lens  (Fig.  3- 
1A).  At  least  20  micrographs  per  slide  were  taken  using  a  Zeiss  Axiocam  HRC  digital 
camera  and  Axiovision  4.6.3  software.  Cells  were  counted  by  taking  the  total  filament 
length  and  dividing  by  an  average  cell  length.  Total  filament  lengths  were  measured  us¬ 
ing  an  automated  image-processing  script  in  MATLAB  (Math Works)  using  scripts  from 
Kovesi  (1999).  In  the  automated  image-processing,  micrographs  were  filtered  to  reduce 
noise  using  a  threshold  value,  phase  congruency,  or  a  fast  Fourier  transform  (FFT),  and 
then  the  edges  were  detected  by  hysteresis  thresholding  (Kovesi,  1999).  The  result  is 
a  binary  image  (Fig.  3- IB).  The  area  of  the  filaments  was  calculated  as  the  total  num¬ 
ber  of  white  pixels  in  the  binary  image  and  total  filament  length  was  calculated  by  di- 
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Figure  3-1:  Image  processing  for  automated  cell  counts.  (A)  Greyscale  epifluorescence 
image  of  T.  thiebautii  II-3.  (B)  Filtered,  binary  image  of  (A).  Cell  counts  are  calculated 
from  the  white  area  in  the  processed  image.  White  bar  in  (A)  represents  200  /j m. 


viding  the  area  by  the  width  of  a  filament  in  the  binary  image.  The  automated  counting 
method  was  validated  by  manually  measuring  filament  lengths  in  98  micrographs  from  T. 
erythraeum  IMS  101  and  T.  thiebautii  II-3  using  ImageJ  software  (National  Institutes  of 
Health,  http://rsb.info.nih.gov/ij/). 

DNA  extraction.  Due  to  problems  with  PCR  inhibition,  a  variety  of  DNA  extraction 
protocols  were  tested.  Filters  were  aseptically  transfered  to  a  2  mL  bead  beating  tube 
(BioSpecs  or  MoBio)  with  two  glass  6.4  mm  beads  (BioSpecs),  four  stainless  steel  beads 
(MoBio),  or  garnet  beats  (MoBio).  To  shred  the  filter  and  homogenize  the  samples,  lysis 
buffer  was  added  as  described  below  and  tubes  were  beat  on  high  using  a  Mini-Beadbeater- 
8  (BioSpecs)  for  two  min,  iced  for  one  min,  beat  for  two  min,  iced  for  one  min,  and  beat 
for  one  min.  All  extracts  and  bead  tubes  with  filters  and  lysates  were  stored  at  -20°C. 

Hot  tris  method.  Samples  were  extracted  according  to  Zinser  et  al.  (2006)  with 
some  modifications.  Filters  were  transfered  to  bead  tubes  with  glass  beads  and  650  /jL 
pH  8.0,  10  mM  Tris-HCl  and  beat  as  described  above.  Avoiding  the  filter,  250  /jL  were 
transfered  to  1.5  mL  microcentrifuge  tubes  and  incubated  at  95°C  for  15  min. 
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Qiagen  tissue  method.  Samples  were  extracted  using  the  DNeasy  Blood  &  Tis¬ 
sue  Kit  (Qiagen)  following  the  protocol  for  animal  tissue.  Filters  were  transfered  to  bead 
tubes  with  glass  beads  and  650  pL  buffer  ATL,  and  beat  as  described  above.  250  pL  lysate 
were  transfered  to  a  new  1.5  mL  microcentrifuge  tube  and  extracted  following  the  manu¬ 
facturer’s  instructions.  DNA  was  eluted  from  the  Mini  Spin  column  sequentially  with  150 
pL  and  100  pL  buffer  AE,  combining  the  eluates. 

Qiagen  plant  methods.  Samples  were  extracted  using  the  DNeasy  Plant  Mini  Kit 
and  the  high-throughput  DNeasy  96  Plant  Kit  (Qiagen).  Filters  were  transfered  to  bead 
tubes  with  glass  beads  and  650  pL  buffer  API,  and  beat  as  described  above.  250  pL  lysate 
were  transfered  to  a  new  1.5  mL  microcentrifuge  tube  or  the  DNeasy  microtubes  and  ex¬ 
tracted  following  the  manufacturer’s  instructions.  DNA  was  eluted  from  the  Mini  Spin 
column  with  150  pL  and  100  pL  buffer  AE  or  from  the  DNeasy  96  plate  twice  with  100 
pL,  combining  the  eluates. 

MoBio  plant  method.  Samples  were  extracted  using  the  Power  Plant  Kit  (MoBio). 
Filters  were  transfered  to  bead  tubes  with  stainless  steel  beads  and  650  pL  Power  Plant 
Bead  Solution,  and  beat  as  described  above.  400  pL  lysate  were  transfered  to  a  new  2.0 
mL  microcentrifuge  tube  and  extracted  following  the  manufacturer’s  instructions.  DNA 
was  eluted  twice  from  the  spin  filter  with  50  pL  buffer  PB6,  combining  the  eluates. 

MoBio  soil  method.  Samples  were  extracted  using  the  Soil  Kit  (MoBio).  Filters 
were  transfered  to  bead  tubes  with  garnet  beads  and  750  pL  lysis  buffer,  and  beat  as  de¬ 
scribed  above.  250  pL  lysate  were  transfered  to  a  new  2.0  mL  microcentrifuge  tube  and 
extracted  following  the  manufacturer’s  instructions.  DNA  was  eluted  twice  from  the  spin 
filter  with  50  pL  buffer  C6,  combining  the  eluates. 

Xanthogenate  method.  Samples  were  extracted  using  the  xanthogenate-SDS  method 
described  by  Tillett  and  Neilan  (2000).  Filters  were  transfered  to  bead  tubes  with  glass 
beads  and  650  pL  TER,  beat  as  described  above,  and  extracted  (Tillett  and  Neilan,  2000). 
DNA  pellets  were  washed  three  times  with  75%  ethanol,  air-dried,  and  re-suspended  in 
200  pL  sterile  18.2  MO  water. 
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Quantitative  polymerase  chain  reaction  (qPCR). 


iCycler  method.  Reactions  were  performed  in  an  iCycler  (BioRad)  with  25  //L  reactions 
either  in  8-strip  0.2  mL  PCR  tubes  (Fisher)  with  optically  clear  flat  cap  strips  (BioRad)  or 
white  96- well  plates  (ABgene)  with  optically  clear  Microseal  ‘B’  film  (BioRad).  Each  re¬ 
action  contained  1  x  iQ  SYBR  Green  Supermix  (BioRad),  0.5  //M  each  of  either  primers 
hetR-TeryF  and  hetR-TeryR  for  clade  I  and  primers  hetR-ThiebF  and  hetR-ThiebR  for 
clade  II  (Table  3.1),  and  10  //L  template.  Cycling  conditions  were  one  cycle  at  95°C  for  3 
min,  40  cycles  of  95°C  for  30  s,  62°C  for  30s  for  clade  I  or  63°C  for  30s  for  clade  II,  and 
an  additional  step  for  clade  II  of  78°C  for  15  s  to  melt  primer  dimers;  followed  by  a  melt 
curve  analysis.  Samples  were  run  in  triplicate  and  cloned  standards  were  run  in  duplicate. 
The  threshold  cycle  ( Q )  was  determined  by  the  cycle  at  which  the  amplification  curve 
crossed  a  threshold  fluorescence  value,  calculated  by  the  best  fit  to  the  standard  curve  of 
diluted  cloned  standards  by  the  iCycler  software.  Sample  concentration  was  calculated 
using  the  sample  Ct  and  the  standard  curve.  PCR  efficiency  (E)  was  calculated  as: 

E  =  (10~™  —  1)  x  100%,  (3.1) 

where  m  is  the  slope  of  the  linear  regression  of  C,  to  logio(standard  concentration). 

LightCycler  method.  For  a  more  high-throughput  method,  reactions  were  performed  in 
a  LightCycler-480  thermal  cycler  (Roche)  with  20  /uh  reactions  in  white  3 84- well  plates 
sealed  with  optically  clear  sealing  foil  (Roche).  Each  reaction  contained  1  x  LightCycler 
SYBR  Green  I  Master  (Roche),  0.25  //M  each  of  HPLC-purified  primers  hetR-TeryF  and 
hetR-TeryR  for  clade  I  and  primers  hetR-ThiebF  and  hetR-ThiebR  for  clade  II,  and  2  /uL 
template.  Master  mix  and  DNA  templates  were  loaded  on  the  384-well  plate  using  a  Pre¬ 
cision  2000  liquid  handler  (Biotek  Instruments).  Cycling  conditions  were  95°C  for  15 
min;  45  cycles  of  95°C  for  20  s,  63°C  for  30,  72°C  for  5  s,  and  a  primer  dimer-melting 
step  for  clade  II  of  76°C  for  10  s;  followed  by  a  melt  curve  analysis.  Samples  and  stan- 
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dards  were  run  in  triplicate.  Crossing  point  (Cp)  was  determined  by  the  maximum  second 
derivative  method  by  the  LightCycler  software.  Sample  concentration  was  calculated  us¬ 
ing  the  sample  Cp  and  the  standard  curve.  PCR  efficiency  was  calculated  as  in  Equation 
3.1,  where  m  is  the  slope  of  the  linear  regression  of  Cp  to  logio(standard  concentration). 

Results 

qPCR  assay  development.  Cell  counts  done  using  the  automated  method  agreed  with 
counts  done  manually  (Fig.  3-2).  The  slope  of  the  linear  fit  of  manually  counted  total  fil¬ 
ament  length  versus  automated  total  filament  length  was  1.15  ±  0.03.  The  slope  for  the 
log io-log io  transformation  was  0.98  ±  0.02,  close  enough  to  1  to  not  transform  the  auto¬ 
mated  numbers  before  using  them. 


Automated  total  filament  length  (urn)  Log10(Automated  total  filament  length) 


Figure  3-2:  Comparison  of  the  performance  of  the  automated  cell  count  method  to  the 
manually  counted  method.  Solid  line  indicates  a  linear  fit  to  the  data.  Dashed  line  indi¬ 
cates  a  1:1  relationship.  (A)  Finear  fit  of  manual  total  filament  length  versus  automated 
filament  lengths:  y  =  (1.15  ±0.03)  x  x  +  10.9  ±  20.5  (R2  =  0.9263).  (B)  Finear  fit  of  the 
log  io  transformed  data:  log\o(y )  =  (0.98  ±0.02)  x  log\o(x)  ±0.10  ±0.01  (R2  =  0.9501). 

Representative  standard  curves  for  the  qPCR  assays  are  shown  in  Fig.  3-3.  PCR 
efficiencies  for  the  qPCR  assays  ranged  from  93  -  102  %.  Melting  temperatures  (Tm) 
were  slightly  higher  for  the  iCycler  method  than  for  the  FightCycler  method  (Fig.  3-4). 
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Tm  for  the  iCycler  (LightCycler)  clade  I  assay  were  80.5°C  (77.9°C)  for  the  standards, 
80.8°C  (78.1°C)  for  T.  erythraeum  IMS101,  and  79.9°C  (78.3°C)  for  T.  contortum  21- 
74  (Fig.  3-4A,  C).  The  Tm  for  T.  contortum  21-74  in  the  iCycler  method  was  slightly 
lower  than  that  of  the  standard  while  both  of  the  cultured  melt  curves  were  nearly  iden¬ 
tical  to  the  standard  in  the  LightCycler  method  (Fig.  3-4A,  C).  For  the  clade  II  assay,  iCy¬ 
cler  (LightCycler)  T,„’s  were  79.6°C  (77.04°C)  for  the  standards,  79.6°C  (77.3°C)  for  T. 
thiebautii  II-3,  79.0°C  (11. TC)  for  T.  tenue  tenue,  and  79.6°C  (77.4°C)  for  T.  tenue  H9-4 
(Fig.  3-4B,  D).  The  Tm  for  T.  tenue  tenue  was  lower  than  the  standard  Tm  for  the  iCycler 
method  and  higher  than  the  standard  for  the  LightCycler  method  while  the  melt  curves  for 
T.  thiebautii  II-3  and  T.  tenue  H9-4  were  identical  to  each  other  and  nearly  identical  to  the 
standard  (Fig.  3-4B,  D).  The  melt  curves  for  the  cultures  in  clade  II  showed  evidence  of 
primer  dimers  in  the  LightCycler  method  (Fig.  3-4D). 

DNA  extractions  and  serial  dilutions.  Serial  dilutions  of  T.  erythraeum  IMS  101  ampli¬ 
fied  well  when  extracted  using  the  hot  tris  method  described  in  Zinser  et  al.  (2006)  (data 
not  shown).  However,  extracts  of  T.  tenue  H9-4  showed  PCR  inhibition  (Fig.  3-5A).  Al¬ 
ternative  methods  of  DNA  extraction  were  tested  to  determine  a  suitable  method  that  did 
not  show  significant  inhibition.  The  MoBio  plant  method  yielded  the  highest  amplifica¬ 
tion  (Fig.  3-5).  For  high  concentrations  of  T.  tenue  H9-4,  the  hot  tris  and  Qiagen  Tissue 
methods  did  not  yield  a  PCR  product  (Fig.  3-5 A).  At  all  concentrations  of  H9-4,  extrac¬ 
tions  from  the  Qiagen  plant  methods  did  not  amplify  as  well  as  extractions  from  the  Mo- 
Bio  plant  method.  The  xanthogenate  and  MoBio  soil  extractions  did  not  amplify  as  well 
as  the  MoBio  plant  extractions  did  for  T.  tenue  H9-4  and  T.  erythraeum  IMS  101  (Fig.  3- 
5B). 

Serial  dilutions  of  Trichodesmium  cultures  extracted  with  the  MoBio  plant  method 
did  not  result  in  a  consistent  relationship  between  cell  counts  and  copies  of  hetR  as  mea¬ 
sured  by  the  qPCR  assay  (Fig.  3-6A-B).  Linear  fits  of  log io -transformed  copies  of  hetR 
rnL-1  to  logio-transformed  cells  mL”1  were  for  IMS101,  log io(y)  =  (1.40  ±  0.02)  x 
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Figure  3-3:  Standard  curves  from  the  qPCR  assays.  The  y-axis  represents  either  threshold 
cycle  (iCycler  method)  or  crossing  point  (LightCycler  method)  and  the  x-axis  repre¬ 
sents  logio(copies  of  hetR).  (A)  Clade  I,  iCycler  method  ( E  —  100.1%,  R2  =  1.000);  (B) 
clade  II,  iCycler  method  ( E  =  96.6%,  R2  =  0.999);  (C)  clade  I,  LightCycler  method  ( E  = 
98.5%,  R2  =  0.994);  (D)  clade  II,  LightCycler  method  ( E  =  98.9%,  R2  =  0.995). 

logw{x)  -0.66  ±0.02  (R2  =  0.9998,  Fig.  3-6A);  for  21-74,  logw(y)  =  (1.84  ±0.18)  x 
logio(x)  —  1.44  ±  0.15  (R2  =  0.9907,  Fig.  3-6A);  for  II-3,  logl0(y)  =  (1.68  ±0.23)  x 
logio(x)  —  1.48  ±0.18  (R2  =  0.9823,  Fig.  3-6B);  fortenue,  /ogio(y)  =  (1.21  ±0.48)  x 
logw{x)  -0.00 ±0.39  (R2  =  0.8627,  Fig.  3-6B);  andforH9-4,  logw(y)  =  (0.86 ±0.10)  x 
log io(x)  ±  1.11  ±0.08  (R2  =  0.9873,  Fig.  3-6B). 

In  mixed  culture  samples,  the  qPCR  method  was  able  to  accurately  measure  the 
components  of  the  mixture.  Dilutions  of  T.  erythraeum  IMS  101  were  mixed  with  dilu¬ 
tions  of  T.  tenue  tenue  and  T.  contortion  21-74  was  mixed  with  T.  thiebautii  II-3.  Fig.  3-7 
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Figure  3-4:  Melting  curves  from  the  qPCR  assays.  The  y-axis  represents  the  negative 
of  the  derivative  of  relative  fluorescence  unit  with  respect  to  temperature  and  the  x-axis 
represents  temperature  (°C).  (A)  Clade  I,  iCycler  method;  (B)  clade  II,  iCycler  method; 
(C)  clade  I,  LightCycler  method;  (D)  clade  II,  LightCycler  method. 


shows  a  linear  fit  of  the  logio  transformed  data  of  measured  qPCR  results  versus  expected 
qPCR  results  based  on  measurements  of  monocultures.  The  slope  of  the  fit  has  a  95% 
confidence  interval  of  0.70  -  1.00  (R2  =  0.9956),  indicating  a  one-to-one  relationship. 
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Dilution  of  H9-4 


H9-4  IMS101 

Strain 


Figure  3-5:  Comparison  of  DNA  extraction  methods  for  (A)  serial  dilutions  of  T.  tenue 
H9-4,  and  for  (B)  replicate  samples  of  T.  tenue  H9-4  and  T.  erythraeum  IMS  101.  “HT”  = 
hot  tris  (cross-hatched),  “QT”  =  Qiagen  DNeasy  Tissue  kit  (diagonal  stripes),  “QP-mini” 
=  Qiagen  DNeasy  Plant  mini  kit  (white),  “QP-96”  =  Qiagen  DNeasy  Plant-96  kit  (hori¬ 
zontal  stripes),  “MP”  =  MoBio  Plant  kit  (black),  “XG”  =  xanthogenate  (vertical  stripes), 
and  “MS”  =  MoBio  Soil  kit  (dotted).  Error  bars  represent  the  standard  deviation  of  2-3 
replicates. 
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Figure  3-6:  Linear  fit  of  qPCR  results  from  serial  dilutions  of  Trichodesmium  cultures 
using  the  LightCycler  method  in  logio(copies  of  hetR  rnL-1)  versus  logio(cells  mL_1). 
Dotted  lines  indicate  a  1:1  relationship.  (A)  Clade  I  strains:  T.  contortum  21-74  and  T. 
erythraeum  IMS  101.  (B)  Clade  II  strains:  T.  tenue  H9-4,  T.  thiebautii  II-3,  and  T.  tenue 
tenue. 


Figure  3-7:  Mixed  culture  qPCR  results. 
Measured  results  (logio(copies  of  hetR 
mL-1))  were  plotted  against  the  expected 
results  based  on  serial  dilutions  of  mono¬ 
cultures  (logio(copies  of  hetR  mL-1)). 
Dilutions  of  T.  contortum  21-74  were  mixed 
with  dilutions  of  T.  thiebautii  II-3,  and 
dilutions  of  T.  erythraeum  IMS  101  were 
mixed  with  dilutions  of  T.  tenue  tenue.  The 
solid  line  is  the  linear  regression  fit  through 
the  logio  -  transformed  data:  log io(y) 

(0.85  ±0.07)  x  logio(x)  ±0.70±0.09  (R2  = 
0.9556).  The  dashed  line  represents  a  1:1 
relationship. 
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Discussion 


Previous  qPCR  methods  have  targeted  Trichodesmium  at  the  genus  level.  In  this  study, 

I  developed  a  qPCR  assay  that  can  distinguish  the  two  major  clades  of  Trichodesmium 
by  targeting  the  heterocyst  differentiation  gene  hetR.  HetR  has  been  widely  studied  in 
Trichodesmium  (Janson  et  al.,  1999;  El-Shehawy  et  al.,  2003;  Lundgren  et  al.,  2005). 

HetR  is  more  variable  than  highly  conserved  regions  such  as  16S  rDNA,  making  it  a  bet¬ 
ter  candidate  for  distinguishing  between  clades  of  Trichodesmium  with  qPCR.  Species 
of  Trichodesmium  are  >97%  identical  in  the  16S  rDNA  region  while  the  inter-clade 
similarity  for  hetR  was  90-91%  (Chapt.  2).  Sequences  of  hetR  in  Trichodesmium  are 
also  more  variable  than  nifH  sequences,  which  were  93  -  94  %  similar  between  clades 
(Lundgren  et  al.,  2001).  ITS  has  been  used  as  a  target  for  qPCR  in  cyanobacteria  (Ahlgren 
et  al.,  2006).  ITS  and  hetR  have  similar  inter-  and  intra-clade  percent  identities  in  Tri¬ 
chodesmium  (89  -  90%  and  97  -  100%,  respectively).  However,  clade- specific  primers 
designed  for  ITS  by  AllelelD  were  rated  “poor"  by  the  software  and  were  incapable  of 
producing  a  PCR  product  when  tested  in  silico  by  Mac  Vector.  T.  erythraeum  IMS  101 
genome  has  two  copies  of  ITS,  complicating  enumeration.  HetR  is  specific  to  filamen¬ 
tous,  diazotrophic  cyanobacteria,  reducing  the  chances  of  amplifying  non -Trichodesmium 
DNA  in  the  qPCR  assay  compared  to  ITS  or  nifH.  The  initial  objective  of  this  study  was 
to  target  the  four  groups  of  Trichodesmium  found  in  Janson  et  al.  (1999)  and  Lundgren 
et  al.  (2005).  However,  in  Chapt.  2,  we  found  that  the  diversity  of  our  cultured  represen¬ 
tatives  spanned  only  two  major  clades,  and  species  within  each  clade  were  too  genetically 
similar  to  design  a  species-specific  assay. 

I  developed  an  effective  qPCR  assay  which  was  able  to  distinguish  between  the 
two  clades  of  Trichodesmium.  PCR  efficiency  was  high,  cross-amplification  was  low,  and 
melting  temperatures  of  cultured  templates  were  reasonably  close  to  the  melting  tempera¬ 
tures  of  standards.  The  primers  for  clade  II  occasionally  formed  primer  dimers.  To  reduce 
the  dimers,  I  used  HPLC-purified  primers  and  changed  to  a  hot-start  method  (LightCycler 


86 


method)  to  use  with  the  liquid  handler  since  reaction  mixtures  may  stay  at  room  temper¬ 
atures  for  an  extended  period  of  time.  To  ameliorate  the  effect  of  the  dimers  on  quantita¬ 
tion,  I  added  a  melting  step  to  the  clade  II  assay.  Reducing  the  effects  of  primer  dimers 
ensures  that  the  SYBR  green  qPCR  assay  is  accurate. 

The  relationships  between  cell  counts  and  copies  of  hetR  for  cultured  strains  of 
Trichodesmium  were  not  consistent.  There  were  multiple  copies  of  hetR  per  cell.  Slopes 
varied  from  0.86  ±  0.10  to  1.84  ±  0.18,  so  in  several  strains,  copies  of  hetR  cell-1  in¬ 
creased  as  cell  counts  increased.  There  are  several  possible  sources  of  the  disparity  be¬ 
tween  gene  copy  and  cell  counts.  (1)  Cyanobacteria  often  have  multiple  copies  of  their 
genome,  and  up  to  20  copies  per  cell  have  been  observed  in  Anabaena  sp.  strain  PCC 
7120  (Herdman  et  al.,  1979;  Cohen  and  Gurevitz,  2006).  Trichodesmium  spp.  may  also 
contain  several  chromosomes  per  cell.  (2)  The  cells  in  the  variety  of  Trichodesmium  cul¬ 
tures  sampled  may  have  been  in  different  cell  cycle  stages.  While  Trichodesmium  cell 
division  peaks  early  in  the  dark  period,  cells  divide  throughout  the  light-dark  cycle  (Sandh 
et  al.,  2009).  (3)  T.  tenue  H9-4  had  colonies,  and  counting  cells  within  colonies  using  mi¬ 
croscopy  can  be  erroneous.  (4)  DNA  extraction  at  low  cell  concentrations  may  be  ineffi¬ 
cient  due  to  interaction  of  the  trichomes  with  the  filter  and  the  beads.  (5)  Cell  count  slides 
were  prepared  using  a  25  mm  glass  filter  tower.  The  meniscus  from  this  small  tower  inter¬ 
acted  with  trichomes  and  colonies,  pulling  them  toward  the  glass.  In  spite  of  efforts  to  dis¬ 
tribute  cells  evenly  throughout  the  filter,  cells  tended  to  be  more  concentrated  around  the 
outside  of  the  field  than  in  the  middle.  Cell  counts  were  underestimated,  and  this  imbal¬ 
ance  of  distribution  can  increase  with  increasing  cell  numbers.  In  the  future,  the  meniscus 
effect  can  be  reduced  by  using  a  47  mm  filter  tower.  Due  to  variability  in  relationships  be¬ 
tween  cells  and  gene  copy,  results  from  qPCR  will  be  reported  in  terms  of  copies  of  hetR 
rather  than  cells. 

Some  Trichodesmium  cultures  showed  evidence  of  PCR  inhibition,  so  I  tested  a 
variety  of  DNA  extraction  methods  to  determine  which  method  minimized  inhibition. 
Cyanobacteria  produce  polysaccharides  and  phenolics  that  can  inhibit  PCR  (Tillett  and 
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Neilan,  2000;  Bessetti,  2007;  Radstrom  et  al.,  2008).  The  xanthogenate  method  was  de¬ 
signed  to  extract  DNA  from  cyanobacteria  (Tillett  and  Neilan,  2000),  but  has  been  shown 
to  inhibit  qPCR  (Foster  et  al.,  2007).  The  MoBio  soil  method,  which  specializes  in  re¬ 
moving  compounds  that  inhibit  PCR,  also  resulted  in  low  amplification.  Low  amplifica¬ 
tion  may  be  due  to  a  combination  of  DNA  extraction  efficiency  and  PCR  inhibition.  The 
MoBio  plant  method  gave  the  best  results  in  our  qPCR  assay  when  compared  to  all  the 
other  methods  tried  and  will  be  used  in  future  field  studies. 
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CHAPTER  4 


The  Distribution  of  Trichodesmium  field 
populations 


ABSTRACT:  A  real-time  quantitative  polymerase  chain  reaction  (qPCR)  assay  targeting 
two  different  clades  of  Trichodesmium  was  applied  to  depth  profiles  in  the  Atlantic  and 
Pacific  Oceans.  Clade  II  ( Trichodesmium  thiebautii,  Trichodesmium  tenue,  Trichodesmium 
hildebrandtii,  and  Trichodesmium  pelagicum)  was  typically  more  numerous  in  the  open 
ocean  than  clade  I  ( Trichodesmium  erythraeum  and  Trichodesmium  contortum).  Concen¬ 
trations  of  both  clades  were  higher  in  the  Atlantic  than  in  the  Pacific  during  the  timeframe 
of  this  study.  Neither  clade  correlated  with  surface  phosphorus  (P)  or  iron  (Fe)  concen¬ 
trations,  but  clade  I  had  deeper  and  colder  distributions  than  clade  II.  Some  sequences 
of  qPCR  products  matched  varieties  of  Trichodesmium  previously  found  in  the  field  but 
not  in  culture.  Numerical  models  parameterizing  N2  fixation  are  often  based  on  clade  I, 
while  clade  II  was  more  numerically  dominant  in  this  study.  Model  results  and  our  under¬ 
standing  of  Trichodesmium  and  global  N2  fixation  may  be  distorted  if  there  are  important 
ecological  differences  between  clades  I  and  II. 
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Introduction 


Cyanobacteria  of  the  genus  Trichodesmium  are  responsible  for  a  large  fraction  of  the  N2 
fixation  in  oligotrophic  oceans  (Capone  et  al.,  1997).  There  are  two  major  clades  of  Tri¬ 
chodesmium:  clade  I  comprising  of  Trichodesmium  erythraeum  and  Trichodesmium  con- 
tortum,  and  clade  II  comprising  of  Trichodesmium  thiebautii,  Trichodesmium  tenue,  Tri¬ 
chodesmium  hildebrandtii,  and  Trichodesmium  pelagicum  (Chapt.  2). 

The  distribution  of  Trichodesmium  spp.  can  be  limited  by  a  variety  of  physical, 
chemical,  and  biological  factors.  Trichodesmium  thrives  in  stratified,  warm  water  with 
high  light  (Capone  et  al.,  1997;  LaRoche  and  Breitbarth,  2005).  Nutrients  such  as  phos¬ 
phorus  (P)  and  iron  (Fe)  affect  Trichodesmium  growth  and  N2  fixation  rates  (Wu  et  al., 
2000;  Moore  et  al.,  2004).  Factors  affecting  the  mortality  of  Trichodesmium  are  not  as 
well  understood.  Grazing  on  Trichodesmium  is  considered  to  be  low  (Carpenter  et  al., 
1987),  and  mortality  due  to  viral  lysis  may  be  important  but  is  not  well-constrained  (Hew- 
son  et  al.,  2004).  Evidence  for  apoptosis  has  been  found  in  both  culture  and  field  popula¬ 
tions  of  Trichodesmium  in  response  to  nutrient  stress  and  high  light  (Berman-Frank  et  al., 
2004,  2007),  but  the  function  of  this  method  of  mortality  and  its  extent  in  field  popula¬ 
tions  are  not  well-understood. 

Trichodesmium  spp.  often  form  extensive  “blooms”  that  are  visible  from  space 
(Capone  et  al.,  1997).  These  regions  of  high  density  occur  in  warm  water  (>  26°C)  with 
low  windspeed  (Carpenter  and  Capone,  1992).  Rather  than  being  a  product  of  high  growth 
and  low  mortality,  these  surface  slicks  of  Trichodesmium  arise  from  the  accumulation  of 
biomass  due  to  periods  of  low  mixing  followed  by  Langmuir  circulation.  Trichodesmium 
cells  are  viable  during  the  early  stages  of  the  slick  (Carpenter  and  Capone,  1992),  but  be¬ 
come  inactive,  possibly  due  to  photooxidative  stress  and  autolysis  (Berman-Frank  et  al., 
2004). 

Field  populations  of  Trichodesmium  have  previously  been  enumerated  using  cell 
counts  (Letelier  and  Karl,  1996;  Post  et  al.,  2002;  Tyrrell  et  al.,  2003),  the  video  plankton 
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recorder  (VPR)  (Davis  and  McGillicuddy,  2006),  and  real-time  quantitative  polymerase 
chain  reaction  (qPCR)  targeting  the  nitrogenase  enzyme  gene  nifH  (Church  et  al.,  2005; 
Foster  et  al.,  2007;  Goebel  et  al.,  2007).  In  a  depth  profile  measured  with  qPCR  at  Sta. 
ALOHA,  Church  et  al.  (2005)  found  that  Trichodesmium  spp.  abundance  was  high  at 
the  surface  with  a  subsurface  peak  near  the  mixed  layer  depth  (MLD).  During  a  2-year 
study  at  Sta.  ALOHA  coupling  qPCR  and  modeling,  Goebel  et  al.  (2007)  found  that  Tri¬ 
chodesmium  spp.  dominated  N2  fixation  in  the  summer  and  autumn,  and  that  abundance 
was  highly  variable.  Foster  et  al.  (2007)  found  that  Trichodesmium  spp.  abundances  in  the 
western  tropical  Atlantic  Ocean  were  high  in  regions  of  low  surface  nutrients,  high  salin¬ 
ity,  and  high  temperatures.  All  of  these  previous  qPCR  studies  have  targeted  the  genus 
Trichodesmium  as  a  whole  using  nifH.  I  developed  a  method  to  distinguish  the  two  ma¬ 
jor  clades  of  Trichodesmium  by  targeting  the  gene  for  the  heterocystous  differentiation 
protein  hetR  (Chapt.  3). 

In  this  study,  I  used  the  qPCR  assay  developed  in  Chapt.  3  to  enumerate  the  two 
clades  of  Trichodesmium  in  field  samples  collected  from  depth  profiles  in  the  west  Pacific 
warm  pool  from  Hawai’i  to  Australia  (Fig.  4-1  A),  the  Equatorial  Atlantic  Ocean  from 
Cape  Verde  to  Barbados  (Fig.  4-1B),  and  the  South  Pacific  Ocean  between  Australia  and 
New  Caledonia  (Fig.  4- 1C).  I  found  that  clade  II  dominates  open  ocean  populations  and 
that  the  two  clades  have  different  depth  and  temperature  distributions. 

Methods 

Sample  collection,  storage,  and  extraction.  Field  samples  were  collected  during  an 
east-west  transect  of  the  Equatorial  Atlantic  Ocean  in  July,  2006  on  the  R/V  Seward- 
Johnson  (SJ0609),  on  a  north-south  transect  of  the  western  Pacific  Ocean  in  January- 
February,  2007  as  part  of  the  West  Pacific  Warm  Pool  cruise  (WP2,  KM0701)  on  the  R/V 
Kilo  Moana,  and  in  the  South  Pacific  Ocean  east  of  Australia  in  March,  2007  on  the  R/V 
Kilo  Moana  (KM0703).  Sea  water  was  collected  at  specific  depths  from  12  L  Niskin  bot¬ 
tles  mounted  on  a  rosette  equipped  with  a  conductivity,  temperature,  and  depth  sensor 
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Longitude 


Figure  4-1:  Map  of  cruise  tracks  for  (A)  west  Pacific  warm  pool  cruise  KM0701  (Jan¬ 
uary  -  February,  2007),  (B)  Equatorial  Atlantic  cruise  SJ0609  (July,  2006),  and  (C)  South 
Pacific  cruise  KM0703  (March,  2007). 


(CTD).  Mixed  layer  depth  (MLD)  for  each  profile  was  determined  using  potential  den¬ 
sity.  Samples  from  four  depths  down  to  the  deep  chlorophyll  maximum  (DCM)  were 
collected  for  SJ0609  and  KM0703.  Six  depths  down  to  150  m  and  including  the  DCM 
were  collected  from  profiles  of  KM0701.  During  a  Trichodesmium  bloom  between  Stas. 
KM0701-16  and  17,  a  sample  was  taken  from  2  m  using  a  hand-deployed  Niskin  bottle. 

At  the  center  of  the  warm  pool  (Sta.  KM0701-14),  depth  profiles  were  taken  every  4  h  for 
24  h  during  an  Eulerian  diel  station.  Eight  to  twelve  L  of  sea  water  was  gravity-filtered 
through  a  47  mm,  5  /an  polycarbonate  filter  (Millipore).  The  filters  were  contained  in 
polypropylene  in-line  filter  holders  (Advantec  MFS)  and  silicone  tubing  that  were  washed 
with  10%  bleach  and  rinsed  with  distilled  water  between  each  deployment.  The  filtrate 
was  measured  and  the  filters  were  placed  in  1.8  mL  cryo  tubes  (Nunc)  and  frozen  in  liquid 
No.  Samples  were  homogenized  by  beadbeating  and  DNA  was  extracted  with  the  Power- 
Plant  DNA  Isolation  Kit  (MoBio)  using  the  methods  described  in  Chapt.  3.  In  addition  to 
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the  extraction,  50  /jL  of  eluate  was  used  in  the  optional  clean-up  protocol  following  man¬ 
ufacturer’s  instructions  (MoBio),  and  the  final  pellet  was  resuspended  in  50  fj L  of  buffer 
PB6  and  50  /uh  of  molecular-grade  water  (5Prime).  West  Pacific  warm  pool  samples  from 
Stas.  KM0701-2  -  14  were  first  extracted  using  the  DNeasy  96  Plant  Kit  (Qiagen)  as  de¬ 
scribed  in  Chapt.  3,  and  then  stored  lysates  were  extracted  using  the  MoBio  plant  method. 
All  extractions  and  bead  tubes  containing  filters  and  lysates  were  stored  at  -20°C. 

During  the  Atlantic  cruise,  10  L  were  gravity-filtered  on  47  mm,  5  /urn  polycarbon¬ 
ate  filters  (Millipore),  fixed  with  formalin,  mounted  on  slides,  and  stored  at  -20°C  for  cell 
counts.  Slides  were  generously  provided  by  the  Zehr  laboratory  (University  of  Calif omia- 
Santa  Cruz).  Slides  were  viewed  using  epifluorescence  (CY3LP  filter)  with  a  Zeiss  Ax- 
ioplan  2  microscope  under  a  lOx  Plan-Neoflaur  objective  lens.  Eighty  micrographs  per 
slide  were  taken  over  the  full  area  of  the  filter  using  a  Zeiss  Axiocam  HRC  digital  camera 
and  Axiovision  4.6.3  software.  Cells  were  counted  as  described  in  Chapt.  3  either  man¬ 
ually  using  ImageJ  software  (http://rsbweb.nih.gov/ij)  or  using  the  automated  method  if 
there  was  little  epifluorescence  from  other  cyanobacteria. 

Nutrient  analyses  from  SJ0709  and  KM0701  were  sent  to  the  College  of  Oceanic 
and  Atmospheric  Sciences,  Oregon  State  University,  and  dissolved  inorganic  phosphorus 
(DIP)  was  analyzed  using  a  Technicon  AutoAnalyzer  II  by  J.  Jennings  with  a  detection 
level  of  6  nmol  L  1 .  KM0703  nutrients  were  analyzed  by  R.  Horak  at  Georgia  Tech  using 
an  AutoAnalyzer.  Dissolved  nutrient  concentrations  from  SJ0609  were  provided  by  P. 

D.  Chappell  (Chappell,  2009),  from  KM0701  were  provided  by  E.  Webb  (Hynes  et  al., 
2009),  and  from  KM0703  were  provided  by  R.  Horak.  Dissolved  Fe  concentrations  for 
SJ0609  and  KM0701  were  provided  by  P.  D.  Chappell  (Chappell,  2009). 

Quantitative  polymerase  chain  reaction  (qPCR).  Reactions  were  performed  in  tripli¬ 
cate  as  described  in  Chapt.  3  in  a  LightCycler-480  (Roche)  with  20  /jL  reactions  in  white 
384-well  plates  sealed  with  optically  clear  sealing  foil  and  the  LightCycler  SYBR  Green 
Master  (Roche).  Primers  targeting  hetR  were  clade-specific  and  standard  curves  were  de- 
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rived  from  cloned  standards  (Chapt.  3).  PCR  efficiency  of  each  sample  was  determined 
as  described  by  Short  et  al.  (2004)  and  Foster  et  al.  (2009),  using  a  fourth  replicate  of 
samples  spiked  with  0.07  pg  of  cloned  standard.  Samples  with  efficiency  <  95%  were 
considered  to  be  inhibited  and  the  compensated  concentration  was  calculated.  Selected 
qPCR  products  were  purified  using  the  QIAquick  PCR  Purification  Kit  (Qiagen)  and  sent 
to  North  woods  DNA,  Inc.  for  sequencing.  Linear  correlations  between  qPCR  results  and 
ancillary  data  were  tested  for  significance  using  a  t-test  on  the  correlation  coefficient  R  (a 
=  0.05). 

Results 

Transect  hydrography,  nutrients,  and  qPCR.  Open  ocean  Trichodesmium  populations 
were  higher  in  the  Atlantic  than  in  the  Pacific  Ocean  and  generally  dominated  by  clade  II 
(Figs.  4-2  -  4-4).  Stations  with  populations  made  up  almost  exclusively  of  clade  II  were 
Stas.  SJ0609-18  (Fig.  4-2A,  B)  and  KM0703-6  and  12  (Fig.  4-4A,  B).  Notable  excep¬ 
tions  to  clade  II  dominance  were  Stas.  SJ0609-14  (Fig.  4-2A,  B),  KM0701-5  -  12  and 
21  -  23  (Fig.  4-3A,  B),  and  the  bloom  station  during  KM0701,  where  concentrations  of 
clade  I  were  higher  than  clade  II.  Cell  counts  for  the  Atlantic  stations  resulted  in  a  linear 
fit  of  logio(total  copies  hetR  L”1)  to  logio(total  cells  L”1)  of  logioiy )  =  (0.86  ±  0.28)  x 
log \o(x)  +  1.66  ±0.14  (R2  =  0.4104)  (Fig.  4-5). 

Maximum  transect  concentrations  of  Trichodesmium  for  clade  I  and  clade  II  in 
the  Equatorial  Atlantic  were  8. 1  x  104  and  1 .5 x  106  copies  hetR  L  1 ,  respectively  (Fig. 
4-2A,  B).  Ratios  of  clade  ILclade  I  non-zero  concentrations  ranged  from  0.007  -  427.9 
with  only  one  ratio  <  1.  The  MLD  of  the  Equatorial  Atlantic  cruise  varied  from  17-85 
m  with  shallower  depths  at  the  east  and  west  ends  of  the  transect  (Fig.  4-2C-E).  Surface 
temperature  ranged  from  26.8  -  29.1°C  (Fig.  4-2C),  and  surface  salinity  ranged  from  13.2 
-  36.2  PSU  with  a  lens  of  fresher  surface  water  at  the  west  end  of  the  transect  (Fig.  4-2D). 
Surface  DIP  concentrations  ranged  from  0.014  -  0.080  /./mol  F”1  (Fig.  4-2E),  and  surface 
dissolved  Fe  concentrations  ranged  from  0.09  -  1.89  nmol  F^^Chappell,  2009). 
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During  the  west  Pacific  warm  pool  cruise,  maximum  transect  concentrations  for 
clade  I  and  II  were  4. Ox  104  and  6.4x  105  copies  hetR  L_1,  respectively  (Fig.  4-3A,  B). 

At  a  surface  slick  of  Trichodesmium  between  Stas.  KM0701-16  and  17,  surface  con¬ 
centrations  were  8.0x  106  copies  hetR  L_1  for  clade  I  and  1.9  x  106  copies  hetR  L~ 1  for 
clade  II.  Ratios  of  non-zero  clade  ILclade  I  concentrations  ranged  from  0.06  to  392.2,  and 
a  majority  of  the  ratios  were  >  1.  The  MLD  ranged  from  11-  118  m,  becoming  more  shal¬ 
low  as  the  transect  moved  south  (Fig.  4-3C,  D).  Surface  temperature  ranged  from  20.8  - 
30.8°C  with  the  warmest  temperature  at  Sta.  KM0701-14  and  the  coolest  temperature  at 
Sta.  KM0701-24,  the  southernmost  station  and  turning  point  (Fig.  4-3C).  Surface  salinity 
ranged  from  33.9  -  35.9  PSU  with  regions  of  high  salinity  near  the  Equator  and  south  of 
New  Caledonia  (Fig.  4-3D).  Surface  DIP  concentrations  ranged  from  0.050  -  0.324  /rmol 
L'1  with  maximum  [DIP]  at  the  Equator  (Fig.  4-3E)  (Hynes  et  al.,  2009),  and  surface 
dissolved  Fe  concentrations  ranged  from  0.09  -  1.40  nmol  L^1  (Chappell,  2009). 

Maximum  Trichodesmium  concentrations  for  the  South  Pacific  cruise  were  3.1  x  103 
and  5.2x  104  copies  hetR  L_1  for  clade  I  and  clade  II,  respectively  (Fig.  4-4A,  B).  Ratios 
of  clade  ILclade  I  where  concentrations  of  both  clades  were  non-zero  ranged  from  0.2  - 
125.6.  For  the  majority  of  samples,  clade  ILclade  I  was  >  1.  MLD  spanned  21  -  61  m, 
the  deepest  mixed  layers  occurring  at  the  southeastern  end  of  the  transect  (Fig.  4-4C-E). 
Surface  temperature  ranged  from  24.2  -  29.2°C  (Fig.  4-2C),  and  surface  salinity  ranged 
from  34.4  -  35.9  PSU  (Fig.  4-2D).  There  was  a  region  of  relatively  cool,  salty  water  off 
the  east  coast  of  Australia  at  the  northern  end  of  the  transect.  Surface  DIP  concentrations 
ranged  from  below  detection  level  (BDL)  -  0.554  /rmol  L-1  (Fig.  4-2E)  (R.  Horak  et  al, 
unpublished  data). 

Across  all  three  cruises,  there  was  a  low  but  significant  correlation  between  the 
concentrations  of  clade  I  and  clade  II  (Fig.  4-6).  A  linear  fit  through  the  logio-logio  trans¬ 
formed  non-zero  data  of  clade  II  (copies  hetR  L  1 )  versus  clade  I  (copies  hetR  L  1 )  of 
the  combined  data  had  the  form  /ogio(y)  =  0.75  x  log\o(x)  +  1.70  (R2  =  0.27,  p  = 

1.3  x  10  7).  This  relationship  was  not  consistent  for  each  cruise.  The  Equatorial  Atlantic 
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Figure  4-2:  Data  from  the  Equatorial  Atlantic  cruise  SJ0609.  The  x-axes  represent  dis¬ 
tance  along  the  transect  in  km  and  have  been  reversed  to  show  the  data  in  a  west-east 
orientation.  Station  numbers  are  indicated  at  the  top.  The  y-axes  represent  depth  (m). 
Black  circles  connected  by  a  black  line  indicate  mixed  layer  depth  (MLD).  (A)  Plot 
of  qPCR  data  for  Trichodesmium  clade  I.  Colored  squares  represent  logio(copies  hetR 
mL  1).  Black  dots  indicate  the  location  of  the  qPCR  samples.  (B)  Plot  as  described  in 
(A)  for  Trichodesmium  clade  II.  (C)  Potential  temperature  in  °  C.  (D)  Salinity  in  PSU.  (E) 
logio([DIP]  in  /rmol  L-1).  Black  dots  indicate  location  of  samples  taken. 
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Figure  4-3:  Data  from  the  west  Pacific  warm  pool  cruise  KM0701.  The  x-axes  represent 
distance  along  the  transect  in  km.  Station  numbers  are  indicated  at  the  top.  The  y-axes  in 
A  -  D  represent  depth  (m).  Black  circles  connected  by  a  black  line  indicate  mixed  layer 
depth  (MLD).  (A)  Plot  of  qPCR  data  for  Trichodesmium  clade  I.  Colored  squares  repre¬ 
sent  logio(copies  hetR  m L  '  ) .  Black  dots  indicate  the  location  of  the  qPCR  samples.  (B) 
Plot  as  described  in  (A)  for  Trichodesmium  clade  II.  (C)  Potential  temperature  in  °  C.  (D) 
Salinity  in  PSU.  (E)  Bar  graph  of  surface  [DIP]  in  /rmol  L_1). 
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Figure  4-4:  Data  from  the  South  Pacific  cruise  KM0703.  The  x-axes  represent  distance 
along  the  transect  in  km  and  station  numbers  are  indicated  at  the  top.  The  y-axes  rep¬ 
resent  depth  (m).  Black  circles  connected  by  a  black  line  indicate  mixed  layer  depth 
(MLD).  (A)  Plot  of  qPCR  data  for  Trichodesmium  clade  I.  Colored  squares  represent 
logio(copies  hetR  mtA1).  Black  dots  indicate  the  location  of  the  qPCR  samples.  (B)  as 
described  in  (A)  for  Trichodesmium  clade  II.  (C)  Potential  temperature  in  °  C.  (D)  Salinity 
in  PSU.  (E)  logio([DIP]  /./mol  L  1 ).  Black  dots  indicate  location  of  samples  taken. 
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Figure  4-5:  Total  qPCR  results  versus  cell 
counts  of  SJ0609.  The  y-axis  represents 
logio(total  copies  of  hetR  L-1)  and  the  x- 
axis  represents  logio(total  cells  L  1 ).  Linear 
fit  of  the  logio-transformed  data  is  repre¬ 
sented  by  the  solid  line.  The  dashed  line 
represents  a  1:1  relationship. 


Figure  4-6:  Comparison  of  concentrations  of 
the  two  Trichodesmium  clades.  The  x-axis 
represents  logio(copies  hetR  L  1 )  for  clade 
I.  The  y-axis  represents  logio(copies  hetR 
L”1)  for  clade  II.  The  grey  line  indicates  a 
one-to-one  relationship.  Black  lines  repre¬ 
sent  linear  fits  of  the  logio-logio  transformed 
data.  The  solid  line  was  the  fit  through 
the  points  for  all  three  cruises.  The  dotted 
line  was  fit  through  SJ0609  (asterisks),  the 
dashed  line  was  fit  through  KM0701  (cir¬ 
cles),  and  the  dash-dot  line  was  fit  through 
KM0703  (triangles). 


cruise  did  not  have  a  significant  relationship  between  the  concentrations  of  the  two  clades: 
log\o(y)  =  0.004  x  log\o(x)  +  5.07  (R2  =  1.49  x  10  5,  p  =  0.492).  The  correlation  be¬ 
tween  the  two  clades  was  significant  in  the  west  Pacific  warm  pool:  Zogio(y)  =  0.77  x 
Zogio(x)  +  1.27  (R2  =  0.32,  p  =  4.2  x  10~5),  but  the  fit  through  the  South  Pacific  cruise 
data  was  not  significant  at  a  95%  confidence  level:  Zogio(y)  =  0.52  x  log io(x)  4-  2.03 
(R2  =  0.16,  p  =  0.07). 
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Figure  4-7:  Plot  of  Fe  (nmol  L-1)  versus  DIP  (/rmol  L_1)  represented  by  black  dots. 
Circle  color  and  size  represent  Trichodesmium  concentrations  (logio(copies  hetR  L  1 )). 
Colored  circles  along  the  x-axis  represent  samples  for  which  Fe  concentrations  were  not 
available.  (A)  Clade  I,  Pacific  samples;  (B)  clade  II,  Pacific  samples;  (C)  clade  I,  Atlantic 
samples;  (D)  clade  II,  Atlantic  samples. 


Nutrients  and  Trichodesmium.  There  were  no  distinct  correlations  between  concentra¬ 
tions  of  either  clade  of  Trichodesmium  and  DIP  or  Fe  (Fig.  4-7).  Each  ocean  basin  had 
similar  ranges  of  Fe,  but  the  Pacific  samples  (Fig.  4-7A-B)  had  higher  concentrations  of 
DIP  than  the  Atlantic  (Fig.  4-7C-D).  Ratios  of  clade  ILclade  I  were  lower  and  more  vari¬ 
able  in  the  Pacific  (Fig.  4-8A)  than  in  the  Atlantic  (Fig.  4-8B).  There  was  no  correlation 
between  Trichodesmium  clade  ratios  and  Fe  and  high  clade  ILclade  I  ratios  appear  to  be 
associated  with  low  DIP  (Fig.  4-8A-B). 
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Figure  4-8:  Ratio  of  clade  IFclade  I  concentrations  in  nutrient  space.  Y-axis  represents 
Fe  concentration  (nmol  L“  1 )  and  x-axis  represents  DIP  concentration  (/./mol  L  1 ).  Circle 
color  represents  logio-transformed  ratios  of  clade  II  concentrations  (copies  hetR  L  “ 1 )  to 
clade  I  concentrations  (copies  hetR  L“  1 ).  (A)  Pacific  samples,  (B)  Atlantic  samples. 


Temperature,  depth,  and  Trichodesmium.  Clade  I  had  broader  temperature  and  depth 
distributions  than  clade  II  from  all  three  cruises  (Fig.  4-9-4-10).  Neither  clade  amplified 
in  the  few  samples  collected  at  temperatures  <  16°C  (Fig.  4-9A-B,  Fig.  4-10A,  C,  E). 
Clade  I  amplified  in  colder  and  deeper  areas  in  the  Pacific  than  clade  II  (Fig.  4-9A-B).  A 
vast  majority  of  the  Atlantic  samples  were  taken  from  the  mixed  layer,  so  the  ranges  of 
temperature  and  depth  were  limited  (Fig.  4-9C-D).  High  concentrations  of  clade  I  were 
found  throughout  the  sample  space  for  both  ocean  basins  (Fig.  4-9A,  C),  while  concentra¬ 
tions  of  clade  II  were  higher  in  shallow,  warm  water  (Fig.  4-9B,  D).  For  the  samples  from 
16  -  18°C,  18  -  20°C,  and  20  -  22°C  the  proportions  of  clade  I  samples  amplifying  were 
0.33,  0.80,  and  0.85,  respectively  (Fig.  4- 10C).  In  comparison,  the  proportions  of  clade  II 
samples  amplifying  at  the  above  temperatures  were  0.17,  0,  and  0.15,  respectively  (Fig. 
4-10E).  As  depth  increased,  the  proportion  of  clade  I  samples  that  amplified  remained 
above  0.5  while  the  proportion  of  clade  II  samples  that  amplified  decreased  below  0.5 
when  depth  >  30  m  (Fig.  4-10D,  F).  Ratios  of  clade  II:  clade  I  concentrations  were  lower 
in  the  Pacific  than  the  Atlantic  (Fig.  4-11A-B).  In  both  basins,  clade  IFclade  I  ratios  were 
low  at  temperatures  less  than  25°C  and  depths  greater  than  50  m  (Fig.  4-1 1A-B). 
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Figure  4-9:  Plot  of  depth  (m)  versus  temperature  (°C)  represented  by  black  dots.  Circle 
color  and  size  represent  Trichodesmium  concentrations  (logio(copies  hetR  L  ')).  (A) 
Clade  I,  Pacific  samples;  (B)  clade  II,  Pacific  samples;  (C)  clade  I,  Atlantic  samples;  (D) 
clade  II,  Atlantic  samples. 
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Figure  4-10:  Histograms  of  temperatures  (A,  C,  E)  and  depths  (B,  D,  F)  of  qPCR  samples 
from  all  three  cruises.  (A)  Histogram  of  potential  temperature  (°C)  for  all  qPCR  samples. 
(B)  Histogram  of  depths  (m)  for  all  qPCR  samples.  (C)  Clade  I  stacked  bar  graph  of  pro¬ 
portions  of  samples  at  each  temperature  for  positive  qPCR  results  (black)  and  negative 
qPCR  results  (white).  (D)  Clade  I  stacked  bar  graph  of  proportions  of  samples  at  each 
depth  for  positive  qPCR  results  (black)  and  negative  qPCR  results  (white).  (E)  Clade  II 
stacked  bar  graph  of  proportions  of  samples  at  each  temperature  for  positive  qPCR  results 
(black)  and  negative  qPCR  results  (white).  (F)  Clade  II  stacked  bar  graph  of  proportions 
of  samples  at  each  depth  for  positive  qPCR  results  (black)  and  negative  qPCR  results 
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Figure  4-11:  Ratio  of  clade  IFclade  I  concentrations  in  temperature-depth  space.  Y-axis 
represents  depth  (m)  and  x-axis  represents  temperature  (°C).  Circle  color  represents 
logio-transformed  ratios  of  clade  II  concentrations  (copies  hetR  L_1)  to  clade  I  concentra¬ 
tions  (copies  hetR  L  ').  (A)  Pacific  samples,  (B)  Atlantic  samples. 


The  depths  of  maximum  concentration  (zmax)  for  clade  I  were  often  found  deeper 
as  well  as  shallower  than  the  MLD  (Fig.  4-12A,  C),  while  zmax\  for  clade  II  were  typi¬ 
cally  shallower  than  the  MLD  (Fig.  4-12B,  D).  During  the  west  Pacific  Warm  pool  cruise, 
the  zmax  for  all  Trichodesmium  was  mostly  in  the  upper  30  m;  however,  clade  I  had  more 
profiles  with  zmax  >  30  m  than  clade  II  (Fig.  4-13).  Several  profiles  had  double  Trichodesmium 
local  maxima.  For  clade  I,  double  maxima  were  observed  in  Stas.  SJ0609-14,  15,  and 
23;  KM0701-3,  5,  7-10,  13,  18,  19,  21  -  23,  and  30;  and  KM0703-4  (Figs.  4-2A  -  4-4A). 
Double  maxima  for  clade  II  were  observed  in  Stas.  SJ0609-16,  17,  and  20;  KM0701-3, 

16,  and  18;  and  KM0703-11  and  12  (Figs.  4-2B  -  4-4B). 
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Figure  4-12:  Plot  of  depth  of  maximum  concentration  (zmax)  versus  mixed  layer  depth 
(MLD),  represented  by  black  dots.  Circle  color  and  size  represent  maximum  Tri- 
chodesmium  concentrations  for  each  profile  (logio(copies  hetR  L  ')).  Black  line  de¬ 
marcates  zmax  =  MLD.  (A)  Clade  I,  Pacific  samples;  (B)  clade  II,  Pacific  samples;  (C) 
clade  I,  Atlantic  samples;  (D)  clade  II,  Atlantic  samples. 
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Figure  4-13:  Histograms  of  depth  of  maximum  concentration  (zmax)  for  depth  profiles 
from  all  three  transects.  (A)  Clade  I  zmax.  (B)  Clade  II  zmax. 
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Figure  4-14:  Melt  curves  of  qPCR  products.  (A)  Clade  I  melt  curve  of  a  cloned  standard 
and  selected  field  samples.  (B)  Clade  II  melt  curve  of  standard  and  selected  field  sample. 


Melting  types,  temperature,  and  depth.  Melt  curve  analysis  of  qPCR  products  re¬ 
vealed  a  variety  of  melting  temperatures  (Tm)  of  clade  I  products  (Fig.  4- 14 A).  Cloned 
standards  had  Tm  «  78°C.  Many  field  samples  matched  the  standard  Tm  (1-78),  but  many 
had  double  peaks  or  a  single  peak  with  Tm  m  80°C  (1-80).  All  Clade  II  field  products  had 
melt  curves  resembling  cloned  standards  with  T„,  m  78°C,  although  many  had  primer 
dimers  as  indicated  by  a  Tm  ~  73°C  (Fig.  4-14B).  1-78  sequences  closely  resembled  T. 
erythraeum  IMS  101  while  sequences  of  1-80  closely  resembled  the  T.  contortum-T.  tenue 
clade  in  Janson  et  al.  (1999)  (Fig.  4-15A).  Some  clade  II  product  sequences  closely  re¬ 
sembled  the  field  sample  T.  aureum  reported  in  Lundgren  et  al.  (2005)  (Fig.  4-15B).  1-80 
products  from  the  west  Pacific  warm  pool  were  found  in  almost  the  full  range  of  tempera¬ 
ture  (Fig  4-16C)  and  were  more  common  in  deep  than  in  shallow  samples  (Fig.  4-16D). 
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Figure  4-15:  Parsimony  tree  of  qPCR  products  from  select  field  samples,  cultured  Tri- 
chodesmium  strains,  and  published  field  samples.  Each  tree  was  bootstrapped  2000  x.  (A) 
Clade  I  qPCR  products,  cultured  T.  erythraeum  strains,  and  field  samples  from  Janson 
et  al.  (1999),  indicated  by  GenBank  numbers  in  parentheses.  The  tree  was  rooted  with  T. 
thiebautii  II-3.  (B)  Clade  II  qPCR  products,  cultured  T.  thiebautii  strains,  and  field  sam¬ 
ples  from  Lundgren  et  al.  (2005),  indicated  by  GenBank  numbers  in  parentheses.  The  tree 
was  rooted  with  T.  erythraeum  IMS  101. 
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Figure  4-16:  Histograms  of  temperature  (A,  C)  and  depth  (B,  D)  for  clade  I  qPCR  prod¬ 
ucts  of  the  west  Pacific  warm  pool  cruise  KM0701,  separated  by  melting  temperature 
(T m).  (A)  Histogram  of  potential  temperature  (°C)  for  qPCR  samples  with  positive  results 
for  clade  I.  (B)  Histogram  of  depth  (m)  for  qPCR  samples  with  positive  results  for  clade  I. 
(C)  Stacked  bar  graph  of  proportions  of  samples  at  each  temperature  for  T„,  ~  80°C  (1-80, 
black)  and  T,„  ~  78°C  (1-78,  white).  (D)  Stacked  bar  graph  of  proportions  of  samples  at 
each  depth  for  1-80  (black)  and  1-78  (white).  1-80  sequences  resemble  the  T.  contortum-T. 
tenue  clade  of  Janson  et  al.  (1999)  and  1-78  sequences  resemble  clade  I). 
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West  Pacific  warm  pool  diel  cycle.  The  diel  station  in  the  center  of  the  west  Pacific 
warm  pool  (Sta  KM0701-14)  was  an  Eulerian  station,  fixed  in  space.  The  maximum  con¬ 
centration  for  clade  I  was  5.4x  104  copies  hetR  L_1  and  for  clade  II  was  4.8 x  105  copies 
hetR  L-1  (Fig.  4-17A,  B).  Ratios  of  non-zero  clade  ITclade  I  concentrations  ranged  from 
0.005  to  243.9  with  ratios  >  1  near  the  surface  and  ratios  <  1  at  depth.  MLD  ranged  from 
28  -  41  m  with  the  maximum  MLD  occurring  at  about  18:30  (Fig.  4-17C,  D).  Tempera¬ 
ture  and  salinity  remained  relatively  constant  throughout  the  cycle  (Fig.  4-17C,  D).  Pho- 
tosynthetically  active  radiation  (PAR)  increased  after  7:00  h  and  decreased  after  18:00 
h  (Fig.  4-17E).  During  the  light  period,  concentrations  of  clade  II  at  100  m  and  150  m 
decreased,  and  increased  during  the  dark  period  (Fig.  4-17B).  This  pattern  was  not  as 
distinct  for  concentrations  of  clade  I  (Fig.  4-17A).  Clade  I  had  two  populations  during 
the  day  at  5  m  and  75  m,  and  one  population  at  75  m  at  night  (Fig.  4- 17 A).  As  PAR  in¬ 
creased,  the  zmax  for  both  clades  shallowed  and  then  deepened  again  as  PAR  decreased 
(Fig.  4-17E). 
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Figure  4-17:  Data  from  the  diel  cycle  at  Sta  14  of  cruise  KM0701.  The  x-axes  represent 
local  time.  The  y-axes  represent  depth  (m).  Black  circles  connected  by  a  black  line  in¬ 
dicate  mixed  layer  depth  (MLD).  (A)  Pseudocolor  plot  of  qPCR  data  for  Trichodesmium 
clade  I.  Colors  represent  logio(copies  hetR  m L  1 ) .  Black  dots  indicate  the  location  of  the 
qPCR  samples.  (B)  Pseudocolor  plot  as  described  in  (A)  for  Trichodesmium  clade  II.  (C) 
Potential  temperature  in  °  C.  (D)  Salinity  in  PSU.  (E)  logio(PAR  in  /./mol  photons  m  2 
s"1).  Clade  I  zmax  is  indicated  by  squares  connected  by  a  solid  line  and  clade  II  zmax  is 
indicated  by  asterisks  and  a  dashed  line. 
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Discussion 


This  study  is  the  first  attempt  at  quantifying  different  clades  of  Trichodesmium  using 
qPCR.  Concentrations  of  the  two  clades  were  measured  across  horizontal  and  vertical 
gradients  in  the  Atlantic  and  Pacific  Oceans.  The  Pacific  Ocean  has  been  under-sampled 
with  respect  to  Trichodesmium  (LaRoche  and  Breitbarth,  2005),  and  this  study  provided 
a  snapshot  of  field  populations  in  the  Pacific  along  a  large  horizontal  gradient.  Concen¬ 
trations  of  clade  II  were  higher  than  those  of  clade  I  in  all  three  transects,  and  both  clades 
reached  higher  concentrations  in  the  Atlantic  than  in  the  Pacific.  Where  both  clades  were 
found,  clade  II  was  up  to  two  orders  of  magnitude  greater  than  clade  I.  On  the  rare  occa¬ 
sion  that  clade  I  was  more  dominant,  it  was  up  to  two  orders  of  magnitude  greater  than 
clade  II.  The  correlation  found  between  clade  I  and  clade  II  concentrations  was  driven  by 
samples  from  the  west  Pacific  warm  pool.  Concentrations  of  the  two  clades  were  corre¬ 
lated  in  the  Pacific  Ocean,  indicating  that  they  may  be  limited  by  the  same  factors  where 
they  coexist  in  the  Pacific  Ocean.  In  contrast,  the  two  clades  were  not  correlated  in  the 
Equatorial  Atlantic.  The  disparate  distributions  in  the  Atlantic  indicate  that  the  two  clades 
were  limited  by  different  factors. 

The  relationship  between  cell  counts  and  qPCR  results  from  the  Equatorial  At¬ 
lantic  was  not  one-to-one,  but  the  linear  fit  of  the  logio  transformed  data  was  within  the 
range  found  in  cultured  strains  (Chapt.  3).  Cyanobacteria  are  known  to  contain  several 
genome  copies  per  cell  (Cohen  and  Gurevitz,  2006).  For  example,  the  filamentous,  het- 
erocystous  cyanobacteria  Anabaena  spp.  have  been  estimated  to  have  from  five  to  eight 
copies  per  cell  in  Anabaena  variabilis  (Chapman  and  Meeks,  1987)  and  up  to  10  -  20 
copies  per  cell  in  Anabaena  sp.  PCC  7120  (Herdman  et  al.,  1979;  Sakr  et  al.,  2006).  Esti¬ 
mating  cell  numbers  in  Trichodesmium  colonies  is  not  trivial  given  the  diversity  of  colony 
diameter  and  filament  density  within  a  colony.  Cell  counts  within  dense  colonies  can  be 
greatly  underestimated  using  microscopy.  Multiple  genome  copies  per  cell  and  under¬ 
estimation  of  cells  within  colonies  can  explain  the  1-2  orders  of  magnitude  difference 
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between  cell  counts  and  copies  of  hetR  per  volume. 

Trichodesmium  spp.  are  able  to  fulfill  their  need  for  N  through  N2  fixation,  so  they 
tend  to  be  limited  by  other  nutrients  such  as  P  and  Fe.  I  expected  to  see  a  correlation  be¬ 
tween  nutrients  and  Trichodesmium  concentrations.  However,  there  was  no  correlation 
found  between  concentrations  of  Trichodesmium  and  DIP  or  Fe,  for  neither  individual 
transects  nor  the  combined  data.  In  the  case  of  P,  Trichodesmium  spp.  are  able  to  access 
dissolved  organic  phosphorus  (DOP)  such  as  phosphomonoesters  (Sohm  and  Capone, 
2006;  Mather  et  al.,  2008)  and  phosphonates  (Dyhrman  et  al.,  2006).  While  populations 
may  show  P  stress,  cells  can  alleviate  P  limitation  by  the  use  of  DOP  instead  of  DIP. 

While  DIP  concentrations  were  lower  in  the  Atlantic  than  in  the  Pacific  in  the  study,  pop¬ 
ulations  of  Trichodesmium  in  the  Atlantic  may  be  obtaining  their  P  from  organic  sources. 
Measurements  of  Fe  had  the  same  range  in  both  the  Atlantic  and  the  Pacific  (Chappell, 
2009),  but  we  don’t  know  the  flux  of  Fe  from  sources  to  each  region  or  the  amount  of 
Fe  storage  in  Trichodesmium  cells.  Aeolian  input  of  Fe  from  dust  is  greater  in  the  At¬ 
lantic  than  the  Pacific  Ocean  (Wu  et  al.,  2000;  Moore  and  Braucher,  2008),  but  input  of  Fe 
from  continental  margins  is  important  in  both  regions  (Moore  and  Braucher,  2008).  Stas. 
KM0701-17  and  18  were  close  to  the  Vanuatu  Islands,  and  Sta.  KM0701-19  was  close  to 
New  Caledonia.  Chappell  (2009)  found  that  Sta.  KM0701-17  had  high  dissolved  Fe,  pos¬ 
sibly  due  to  advection  of  Fe  from  Aniwa  Island  of  the  Vanuatu  archipelago.  The  potential 
of  these  islands  as  sources  of  nutrients  and  their  influence  on  Trichodesmium  populations 
cannot  be  ignored.  Nutrient  history  and  cell  quotas  may  be  more  important  than  environ¬ 
mental  nutrient  concentrations  to  the  growth  of  Trichodesmium  cells,  as  ambient  nutrient 
concentrations  may  not  reflect  the  nutrient  status  of  the  cells. 

Nutrient  limitation  can  define  phytoplankton  niches.  According  to  classic  resource 
competition  theory,  species  with  different  limiting  resources  can  coexist  (Tilman,  1977). 
Coexistence  of  the  two  clades  of  Trichodesmium  could  be  possible  through  differential 
limitation  by  P  and  Fe.  The  highest  clade  ILclade  I  ratios  were  found  in  areas  of  low  DIP, 
which  is  consistent  with  the  hypothesis  that  clade  I  may  be  more  limited  by  P  while  clade 
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II  may  be  more  limited  by  Fe.  For  example,  Pacific  Sta.  KM0701-10  had  the  highest 
[DIP]  of  the  transect,  high  clade  I  concentration,  and  no  amplification  of  clade  II.  How¬ 
ever,  the  majority  of  the  low  DIP  stations  were  in  the  Atlantic,  and  Atlantic  concentrations 
of  clade  I  as  well  as  clade  II  exceeded  those  of  the  Pacific.  Factors  that  differ  between  the 
Atlantic  and  the  Pacific  other  than  nutrients  might  be  influencing  Trichodesmium  distribu¬ 
tions.  Niche  differentiation  of  Trichodesmium  spp.  is  not  as  simple  as  differential  nutrient 
limitation,  but  further  investigations  are  required  to  understand  the  influence  P  and  Fe 
have  on  the  growth  and  distributions  of  different  Trichodesmium  spp. 

While  there  was  not  a  distinct  spatial  separation  between  the  clades  of  Trichodesmium , 
clade  I  populations  had  colder  and  deeper  ranges  than  clade  II.  The  presence  of  Trichodesmium 
cells  in  deep  or  cold  waters  does  not  necessarily  indicate  growth  under  those  conditions; 
cells  may  have  been  advected  there.  However,  these  temperature  and  depth  ranges  may 
indicate  a  higher  tolerance  for  low  light  and  cold  conditions  in  clade  I  than  in  clade  II. 

Chappell  (2009)  found  that  while  cultured  strains  of  clade  I  had  higher  optimum  tem¬ 
peratures  (28°C)  than  a  clade  II  strain  (26°C),  the  clade  I  strains  had  higher  growth  rates 
across  all  temperatures  measured.  Clade  I  had  both  surface  and  deep  populations,  dis¬ 
cernible  by  T,„  as  1-78  and  1-80.  1-80  may  be  the  T.  contortum  -  T.  tenue  group  sequenced 
in  Janson  et  al.  (1999).  Clade  II  products  also  had  sequences  similar  to  T.  aureum  de¬ 
scribed  in  Lundgren  et  al.  (2005).  The  range  of  depth  for  clade  I  may  be  driven  by  the 
1-80  population,  which  was  the  prominent  melting  type  of  clade  I  at  depth.  The  T.  tenue 
in  Janson  et  al.  (1999)  was  collected  from  the  central  Atlantic  Ocean  at  75  m,  while  our 
culturing  efforts  have  been  mostly  focused  at  the  surface.  Culturing  efforts  must  extend 
deeper  in  order  to  capture  the  full  diversity  of  Trichodesmium  spp.  Having  cultured  rep¬ 
resentatives  of  these  deep  clades  would  enable  the  development  of  more  specific  qPCR 
primers.  The  1-80  population  did  not  have  a  particularly  colder  range  than  1-78,  so  the 
depth  distribution  of  1-80  may  be  driven  by  light  rather  than  temperature.  Depth  analysis 
of  1-80  could  be  done  only  on  the  KM0701  profiles,  which  were  sampled  consistently  to 
150  m.  Even  at  150  m,  the  bottom  of  the  Trichodesmium  layer  was  not  always  reached. 
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In  addition,  distributions  of  Trichodesmium  populations  are  patchy.  Future  studies  should 
aim  for  higher  resolution  and  deeper  samples  than  this  or  previous  studies,  for  example 
using  the  VPR  in  concert  with  qPCR. 

The  deep  range  of  clade  I  was  surprising.  Cultured  strains  of  clade  I  tend  to  be 
highly  buoyant.  In  field  populations,  tuft  colonies,  which  are  typically  associated  with  T. 
erythraeum,  were  found  to  be  concentrated  at  the  surface  while  puff  colonies,  presum¬ 
ably  T.  thiebautii,  were  distributed  at  depth  (Post  et  al.,  2002;  Davis  and  McGillicuddy, 
2006).  With  the  exception  of  Sta.  KM0701-10,  the  majority  of  colonies  during  the  Equa¬ 
torial  Atlantic  and  west  Pacific  warm  pool  cruises,  both  puffs  and  tufts,  belonged  to  clade 
II  (Chappell,  2009).  Most  of  the  biomass  of  clade  I  during  this  study  was  probably  in 
the  form  of  single  trichomes,  which  would  not  have  the  same  buoyancy  as  colonies  and 
would  not  be  detected  using  the  VPR  (Davis  and  McGillicuddy,  2006).  Composition  of 
the  photosynthetic  accessory  pigment  phycoerythrin  (PE)  of  the  two  strains  also  predicts 
clade  I  to  be  a  surface-dweller  and  clade  II  to  be  depth- adapted.  Cultured  strains  of  clade  I 
had  high  amounts  of  the  PE  subcomponent  phycoerythrobilin  (PEB),  and  strains  of  clade 
II  had  high  amounts  of  the  subcomponent  phycourobilin  (PUB)  (Chapt.  2).  PEB  absorbs 
green  wavelengths  of  light  (around  545  nm)  while  PUB  absorbs  blue-green  light  (around 
495  nm).  Since  blue  wavelengths  penetrate  deeper  in  the  water  column  than  green  wave¬ 
lengths,  PUB  should  be  the  preferred  pigment  at  depth  and  clade  II  should  be  more  depth- 
adapted.  However,  clade  I  cells  have  higher  PE  content  than  clade  II  (Carpenter  et  al., 
1993),  which  might  be  an  adaptation  to  low  light.  The  T.  contortum  described  by  Janson 
et  al.  (1999)  was  straw-colored  while  the  T.  tenue  described  was  deeply  pigmented,  so 
the  pigmentation  of  the  deep  1-80  population  of  clade  I  might  have  been  diverse.  Clade  II 
distributions  indicate  a  high  light  requirement  or  at  least  a  preference  for  shallow  waters. 
The  z„MX’s  for  clade  II  were  consistently  shallower  than  the  MLD.  During  the  west  Pacific 
warm  pool  cruise,  clade  II  was  not  found  in  the  warm  but  deep  mixed  layers  of  stations 
north  of  the  Equator,  but  was  found  in  the  cool  but  shallow  mixed  layer  of  the  southern¬ 
most  end  of  the  transect. 
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The  difference  in  Trichodesmium  clade  concentrations  might  be  due  to  differ¬ 
ential  mortality.  Clade  II  is  known  to  produce  more  toxins  than  clade  I,  and  clade  I  is 
preferentially  grazed  by  the  harpacticoid  copepod  Macrosetella  gracilis  (O’Neil  and  Ro¬ 
man,  1994).  If  grazing  mortality  is  significant,  standing  stocks  of  clade  I  could  be  lower 
than  clade  II  even  if  clade  I  growth  is  higher.  Virus-like  particles  (VLP)  found  in  clade 
II  field  samples  from  the  North  Pacific  were  larger  than  the  VLP  found  in  clade  I  (Hew- 
son  et  al.,  2004),  but  rates  of  viral  lysis  for  the  two  clades  are  not  known.  Autolysis  in  the 
two  clades  might  be  triggered  by  different  stressors  or  ranges  of  stressors  such  as  light  or 
nutrient  levels. 

During  the  diel  cycle  at  KM0701-14,  depth  distributions  of  Trichodesmium  were 
more  shallow  at  daytime  than  nighttime.  The  effect  was  much  more  noticeable  in  clade  II 
populations  than  in  clade  I.  This  was  an  Eulerian  station,  so  the  populations  sampled  were 
not  from  the  same  water  mass,  and  we  cannot  ignore  the  possibility  that  the  nighttime 
shallowing  pattern  observed  may  be  from  the  advection  of  a  water  mass  with  different 
properties  and  populations  rather  than  a  response  to  a  diel  cycle.  However,  the  possible 
movement  of  Trichodesmium  from  the  surface  to  below  the  phosphocline  at  night  would 
support  the  P  mining  hypothesis  (Karl  et  al.,  1992;  White  et  al.,  2006b). 

Most  of  the  information  available  about  the  physiology  of  Trichodesmium  is  based 
on  clade  I,  particularly  T.  erythraeum  IMS  101.  Culture  studies  investigating  Trichodesmium 
physiology  in  response  to  nutrients  (Mulholland  et  al.,  2002;  Holl  and  Montoya,  2005), 
trace  metals  (Tuit  et  al.,  2004),  pC02  (Hutchins  et  al.,  2007),  light  (Bell  and  Fu,  2005), 
and  temperature  (Breitbarth  et  al.,  2007)  have  all  been  based  on  cultures  of  T.  erythraeum. 
Consequently,  models  investigating  N2  fixation  incorporate  parameters  based  on  T.  ery¬ 
thraeum  (Moore  et  al.,  2004;  Rabouille  et  al.,  2006;  Goebel  et  al.,  2007).  In  this  study, 
clade  I,  which  includes  T.  erythraeum,  has  been  shown  to  be  less  prominent  than  clade 
II.  If  the  physiology  of  clade  I  and  clade  II  are  significantly  different,  our  understanding 
of  global  N2  fixation  based  on  parameters  from  clade  I  may  be  skewed.  This  study  has 
taken  a  step  toward  examining  niche  differentiation  in  Trichodesmium  and  getting  a  better 
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understanding  of  the  global  distribution  of  phylotypes  of  this  important  N2  fixer. 
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CHAPTER  5 


Cross-basin  comparison  of  phosphorus  stress 
and  nitrogen  fixation  in  Trichodesmium 
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ABSTRACT: 

We  investigated  the  phosphorus  (P)  status  and  N2  fixation  rates  of  Trichodesmium 
populations  from  the  North  Pacific,  western  South  Pacific,  and  western  North  Atlantic. 
Colonies  of  Trichodesmium  were  collected  and  analyzed  for  endogenous  alkaline  phos¬ 
phatase  (AP)  activity  using  enzyme-labeled  fluorescence  (ELF)  and  for  nitrogenase  activ¬ 
ity  using  acetylene  reduction.  AP  hydrolyzes  dissolved  inorganic  phosphate  (DIP)  from 
dissolved  organic  phosphorus  (DOP)  and  is  active  in  Trichodesmium  colonies  experienc¬ 
ing  P  stress.  Across  multiple  stations  in  the  subtropical  North  and  South  Pacific,  there 
was  low  to  moderate  ELF  labeling  in  Trichodesmium,  although  labeling  was  present  in 
other  taxa.  In  contrast,  Trichodesmium  ELF  labeling  in  the  North  Atlantic  ranged  from 
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low  to  high.  Low  ELF  labeling  corresponded  with  high  DIP  concentrations  while  high 
ELF  labeling  occurred  only  at  North  Atlantic  stations  with  DIP  concentrations  <  40  nmol 
L  " 1 ,  indicating  that  Trichodesmium  was  not  experiencing  dramatic  P  stress  in  the  Pacific 
Ocean  while  P  stress  was  evident  in  the  western  North  Atlantic.  However,  nitrogenase 
activity  was  significantly  higher  in  the  P- stressed  western  North  Atlantic  than  in  the  Pa¬ 
cific  Ocean  (0.40  -  1.30  compared  to  0.01  -  0.46  nmol  C2H4  h-1  colony-1).  These  data 
underscore  the  differential  basin-level  importance  of  P  availability  to  Trichodesmium  and 
suggest  that  factors  other  than  P  are  constraining  their  N2  fixation  rates  in  the  Pacific. 

Introduction 

New  nitrogen.  Fixed  nitrogen  (N)  is  considered  to  be  the  proximal  limiting  nutrient 
in  ocean  ecosystems.  For  a  system  in  steady-state,  new  production  is  the  excess  pro¬ 
duction  available  for  export  and  is  a  function  of  new  N  entering  the  system  and  stimu¬ 
lating  production  (Eppley  and  Peterson,  1979).  New  N  comes  from  eddy  diffusion  pro¬ 
cesses,  seasonal  deep  mixing,  atmospheric  deposition,  lateral  advection,  and  N2  fixation 
by  diazotrophs  (Dugdale  and  Goering,  1967).  Upwelled  water  from  the  deep  ocean  con¬ 
tains  both  elevated  dissolved  inorganic  nitrogen  (DIN)  and  dissolved  inorganic  carbon 
(DIC),  so  export  production  from  upwelled  DIN  does  not  generally  result  in  a  significant 
net  removal  of  DIC  from  surface  water.  However,  DIN  from  N2  fixation  is  not  related 
to  the  DIC  supply  from  deep  water  and  can  drive  uptake  of  CO2  and  export  production 
(Hood  et  al.,  2000).  In  global  marine  models,  the  presence  of  N2-fixers  increased  primary 
production  by  diatoms  and  small  phytoplankton  through  the  excretion  of  dissolved  or¬ 
ganic  nitrogen  (DON)  in  N-limited  areas  (Moore  et  al.,  2004).  In  the  subtropical  Atlantic, 
model  results  show  that  N2  fixation  increases  DIN,  phytoplankton  biomass,  primary  pro¬ 
duction,  and  export  flux  (Coles  et  al.,  2004). 

The  cyanobacterium  Trichodesmium  is  an  important  source  of  new  N  in  tropical 
and  subtropical  oceans  (Capone  et  al.,  1997).  In  the  North  Atlantic,  N2  fixation  rates  of 
Trichodesmium  can  equal  or  exceed  the  vertical  flux  of  nitrate  into  surface  waters  (Capone 
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et  al.,  2005),  while  in  the  North  Pacific,  N?  fixation  is  the  source  of  up  to  half  of  the  new 
N  (Karl  et  al.,  1997).  Global  ecosystem  models  use  parameters  from  Trichodesmium  to 
model  No  fixation  (Coles  et  al.,  2004;  Moore  et  al.,  2004;  Moore  and  Doney,  2007).  Elu¬ 
cidating  the  factors  limiting  the  growth  and  N2  fixation  rates  of  Trichodesmium  is  vital  to 
understanding  the  N  cycle  and  its  link  to  the  carbon  cycle. 

Phosphorus  and  iron.  No  fixation  by  Trichodesmium  can  be  limited  by  a  variety  of  fac¬ 
tors  including  light,  temperature,  mixed  layer  depth,  dissolved  oxygen,  and  nutrients  such 
as  phosphorus  (P)  and  iron  (Fe)  (Wu  et  al.,  2000;  Moore  et  al.,  2004).  Like  N,  new  P  can 
reach  the  surface  by  mixing,  upwelling  from  deep  waters,  lateral  advection,  or  deposition; 
however,  there  is  no  process  analogous  to  N2  fixation  to  act  as  a  source  for  bioavailable 
P  (Karl,  2002).  Nitrogenase,  the  enzyme  responsible  for  N2  fixation,  requires  multiple 
Fe  atoms  per  molecule.  Sources  of  Fe  to  open-ocean  environments  are  dust  deposition, 
advection  from  continental  margins,  and  upwelling  (Fung  et  al.,  2000).  Due  to  these  obli¬ 
gate  requirements  and  low  oceanic  concentrations,  P  and  Fe  bioavailability  may  influence 
the  growth  and  N2  fixation  rates  of  Trichodesmium  field  populations  (Berman-Frank  et  al., 
2001;  Sanudo-Wilhelmy  et  al.,  2001).  In  the  central  Atlantic,  N2  fixation  positively  cor¬ 
related  to  the  P  content  of  Trichodesmium  (Sanudo-Wilhelmy  et  al.,  2001)  while  in  the 
eastern  tropical  Atlantic,  N2  fixation  was  co-limited  by  P  and  Fe  (Mills  et  al.,  2004).  In 
the  North  Pacific  Subtropical  Gyre  (NPSG),  P  supply  may  be  the  ultimate  control  on  new 
and  export  production  at  Sta.  ALOHA  (Karl  et  al.,  1997),  but  Fe  rather  than  P  is  likely  to 
limit  N2  fixation  (Wu  et  al.,  2000). 

While  subtropical  gyres  of  both  the  North  Pacific  and  the  North  Atlantic  are  olig- 
otrophic,  the  NPSG  has  higher  dissolved  inorganic  phosphate  (DIP)  concentrations:  3.0 
-  191.3  nmol  L  1  in  the  NPSG  compared  to  0.2  -  1.0  nmol  L  1  in  the  North  Atlantic 
(Hawaii  Ocean  Time-series  (HOT),  http://hahana.soest.hawaii.edu/hot/hotjgofs.html; 

(Wu  et  al.,  2000).  Dissolved  organic  and  particulate  pools  of  P  are  also  higher  in  the 
NPSG  than  the  North  Atlantic  (Ammerman  et  al.,  2003).  In  the  western  South  Pacific, 
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DIP  concentrations  ranged  from  80  to  218  nmol  IA 1  in  oceanic  regions  (Campbell  et  al., 
2005)  and  from  <  30  to  63  nmol  L  1  near  New  Caledonia  (Van  Den  Broeck  et  al.,  2004). 
P  limitation  may  be  related  to  inputs  of  Fe  (Wu  et  al.,  2000).  The  NPSG  has  lower  aeo- 
lian  Fe  input  than  the  subtropical  North  Atlantic:  0.08  -  0.16  //mol  Fe  m  2  d  1  compared 
to  0.2  -  0.8  fj mol  Fe  nA2  d~\  respectively  (Wu  et  al.,  2000).  Open  ocean  dissolved  Fe 
concentrations  in  the  upper  100  m  averaged  720  pmol  L”1  in  the  North  Atlantic,  200 
pmol  L  1  in  the  North  Pacific,  and  110  pmol  IA1  in  the  Equatorial  Pacific  (Moore  and 
Braucher,  2008).  During  a  cruise  in  the  western  South  Pacific,  open  ocean  dissolved 
Fe  concentrations  ranged  from  <  60  to  760  pmol  IA1  while  coastal  dissolved  Fe  con¬ 
centrations  ranged  from  <  60  to  1000  pmol  IA1  (Campbell  et  al.,  2005).  In  addition  to 
low  concentrations,  bioavailability  of  Fe  and  P  can  limit  abundance  and  activity  of  Tri- 
chodesmium,  and  differences  in  the  availability  of  these  nutrients  across  ocean  basins  cre¬ 
ate  a  variety  of  limitation  scenarios. 

Alkaline  phosphatase.  Herein  we  define  nutrient  stress  as  the  exposure  of  a  cell  to  low 
inorganic  nutrient  conditions  which  stimulates  a  response  from  the  cell.  Phytoplankton 
can  launch  a  physiological  response  under  conditions  of  P  stress  to  ameliorate  potential 
limitation.  DIP  is  the  most  readily  available  form  of  phosphorus.  The  utilization  of  dis¬ 
solved  organic  phosphorus  (DOP)  can  prevent  limitation  of  Trichodesmium  growth  and 
N?  fixation  rates  in  low  DIP  environments.  However,  accessing  the  DOP  pool  requires  the 
cell  to  invest  energetic,  carbon,  and  nitrogen  resources  in  enzymes  such  as  alkaline  phos¬ 
phatase  (AP).  We  define  this  as  distinct  from  P  limitation,  which  is  the  reduction  of  either 
growth  or  N2  fixation  rates  as  a  result  of  persistent  stress.  AP  hydrolyzes  DIP  from  DOP 
and  is  commonly  induced  in  phytoplankton  species  experiencing  P  stress  (Stihl  et  al., 
2001;  Dyhrman  et  al.,  2002).  DOP  is  likely  a  significant  source  of  P  for  Trichodesmium 
spp.  (Sohm  and  Capone,  2006;  Mather  et  al.,  2008),  and  the  presence  of  AP  activity  can 
indicate  Trichodesmium  field  populations  responding  to  low  levels  of  DIP  (Stihl  et  al., 
2001).  AP  activity  can  be  measured  using  colorimetric  and  fluorometric  assays  or  de- 
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tected  through  molecular  diagnostic  tools.  Enzyme-labeled  fluorescence  (ELF)  is  a  com¬ 
mercially  available,  cell-specific  molecular  diagnostic  tool  that  detects  AP  activity.  The 
phosphatase  substrate  precipitates  when  cleaved  by  AP  and  fluoresces  green  under  ultra¬ 
violet  light  (Gonzalez-Gil  et  al.,  1998).  In  the  laboratory,  P-depleted  Trichodesmium 
erythraeum  IMS  101  batch  cultures  were  ELF-labeled  while  replete  cultures  were  not 
(Dyhrman  et  al.,  2002).  In  ship-board  incubations,  addition  of  DIP  to  ELF-positive  pop¬ 
ulations  of  Trichodesmium  reduced  ELF  labeling  within  48  h  (Webb  et  al.,  2007).  We 
believe  that  ELF  labeling  of  P-stressed  Trichodesmium  gives  valuable  insight  into  the  in 
situ  physiological  status  of  the  cells. 

In  this  study  we  compare  Trichodesmium  P  stress  and  nitrogenase  activity  among 
tropical  and  subtropical  regions  of  the  western  North  Atlantic,  the  North  Pacific,  and  the 
western  South  Pacific  (Fig.  5-1).  We  collected  colonies  of  Trichodesmium  from  four 
cruises  and  analyzed  endogenous  AP  activity  using  ELF,  measured  nitrogenase  activity 
using  acetylene  reduction,  and  sampled  surface  water  nutrients  to  assess  the  P  status  of 
Trichodesmium  populations.  Our  results  show  that  populations  in  the  western  North  At¬ 
lantic  were  more  P  stressed  but  had  higher  nitrogenase  activities  than  populations  in  the 
North  and  South  Pacific. 

Methods 

Hydrological  context.  Samples  were  collected  along  a  north-south  transect  of  the  Pa¬ 
cific  Ocean  as  part  of  the  Western  Pacific  Warm  Pool  cruise  (WP2,  KM0701)  on  the  R/V 
Kilo  Moana  in  January-February  2007  (Fig.  5-1  A).  As  part  of  the  Marine  Nitrogen  Trans¬ 
formation  (MANTRA)  Biocomplexity  Project,  samples  were  taken  from  south  of  Hawaii 
during  MP09  in  August  2003,  on  the  R/V  Revelle  (Fig.  5-1B).  North  Atlantic  field  sam¬ 
ples  were  collected  near  the  Bermuda  Atlantic  Time-series  Station  (BATS)  in  June  and 
October  2001,  aboard  the  R/V  Endeavor  (EN355  and  EN361,  respectively)  (Fig.  5-1C). 
Temperature,  salinity,  pressure,  and  density  data  from  each  cruise  was  obtained  from  a 
CTD  system.  Mixed  layer  depths  (MLD)  were  estimated  from  potential  density  profiles. 
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Longitude  Longitude 


Figure  5-1:  Maps  of  cruise  tracks:  (A)  KM0701  (circles)  during  January  and  February 
2007,  across  the  west  Pacific  warm  pool;  (B)  MP09  (diamonds)  during  August  2003,  in 
the  North  Pacific;  (C)  EN355  (squares)  during  June  2001,  and  EN361  (triangles)  during 
October  2001,  in  the  western  North  Atlantic.  The  highest  level  of  ELF  labeling  for  all 
colony  types  at  each  station  is  indicated  by  color:  black  symbols  =  ‘+  +’,  gray  symbols  = 
‘+  open  symbols  =  small  dots  =  no  ELF  data. 


Sea  surface  temperature  (SST)  and  salinity  were  obtained  from  the  ships’  underway  data. 

Collection  of  colonies.  A  0.5  m  diameter,  1 30.0-qm  mesh  phytoplankton  net  with  a  30- 
m  line  was  hand  towed  for  10  -  20  minutes.  Tows  for  ELF  were  typically  done  in  the 
morning  while  tows  for  nitrogenase  activity  were  done  around  noon,  the  peak  of  N2  fix¬ 
ation  in  Trichodesmium,  targeting  local  times  of  11:00  h,  12:00  h,  and  13:00  h.  Specimens 
were  immediately  shaded  and  taken  into  an  air-conditioned  laboratory.  Colonies  of  Tri¬ 
chodesmium  spp.  were  collected  with  polyethylene  bulb  pipets  or  plastic  bacteriological 
transfer  loops  (Fisher  Scientific);  separated  by  colony  morphology  into  spherical  ‘puffs,’ 
fusiform  ‘tufts,’  and  occasionally  ‘bowties’;  and  washed  four  times  in  sterilized  seawa- 
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ter  (0.2  -  0.4  //in- filtered  and  microwaved  to  a  rolling  boil).  The  colonies  were  processed 
immediately  for  AP  or  nitrogenase  activity.  Stations  with  no  data  reported  had  too  little 
Trichodesmium  biomass  to  run  the  assays. 

Enzyme-labeled  fluorescence  (ELF).  The  presence  of  AP  activity  was  determined  us¬ 
ing  the  ELF  97  Endogenous  Phosphatase  Detection  Kit  (Molecular  Probes-Invitrogen) 
following  the  improved  protocol  described  in  Dyhrman  et  al.  (2002)  for  two  colonies  of 
each  morphology.  The  slides  were  stored  damp  and  in  the  dark  at  4°C.  ELF  slides  were 
viewed  with  epifluorescent  microscopy  using  a  4,  6-diamidino-2-phenyl-indole  (DAPI) 
long-pass  filter  set  (360  nm).  The  entire  area  of  the  slide  was  examined.  Digital  pictures 
of  Trichodesmium  from  MP09  and  the  Atlantic  cruises  were  taken  at  200X  using  a  Zeiss 
Axioplan2  microscope  equipped  with  a  Zeiss  Axiocam  digital  camera  and  Openlab  soft¬ 
ware  (Impro vision).  Digital  pictures  from  KM0701  were  taken  at  400X  using  a  Zeiss  Ax- 
iostar  microscope  equipped  with  a  MRC5  Axiocam  and  Axiovision  software  (Zeiss). 
Samples  were  rated  ‘-  -’  for  ‘no  labeling,’  ‘+  -’  for  ‘some  labeling’  if  many  filaments 
showed  moderate  labeling  but  were  not  wholly  covered  or  if  there  was  mixture  of  labeled 
and  unlabeled  filaments,  and  ‘+  +’  for  ‘high  labeling’  if  many  filaments  showed  dense  la¬ 
beling  and/or  were  wholly  covered.  For  the  North  Atlantic  samples,  multiple  slides  were 
made  for  each  colony  type  so  one  rating  was  given  based  on  all  the  slides.  Filament  color 
was  noted  as  ‘normal’  if  the  orange  color  of  the  phycoerythrin  autofluorescence  was  vis¬ 
ible  or  ‘dim’  if  the  filaments  appeared  abnormally  dark,  green  or  blue.  Digital  pictures 
of  pennate  diatoms  associated  with  Trichodesmium  colonies  were  taken  at  400X  with  the 
DAPI  long-pass  filter  and  a  chlorophyll  filter  (660  nm)  using  AxioVision  4.2  software. 

Nutrient  Analysis.  For  the  North  and  South  Pacific  (KM0701)  and  the  North  Atlantic 
(EN355  and  EN361),  surface  water  was  collected  using  either  niskin  bottles  from  a  rosette 
or  Go-Flo  bottles,  filtered  through  0.4  /zm  acid-washed  polycarbonate  filters,  and  frozen. 
Samples  for  KM0701  were  sent  to  the  College  of  Oceanic  and  Atmospheric  Sciences, 
Oregon  State  University  and  DIP  was  analyzed  using  a  Technicon  AutoAnalyzer  II  by 
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J.  Jennings  with  a  detection  level  of  6  nmol  L  1 .  North  Atlantic  samples  were  sent  to  the 
Ocean  Data  Facility  of  the  Scripps  Institution  of  Oceanography  (http://sts.ucsd.edu/sts/chem) 
for  analysis  by  a  Skalar  San  Plus  autoanalyzer  with  a  detection  level  of  3  nmol  L-1.  DIP 
concentrations  for  the  North  Pacific  during  MP09  were  provided  by  K.  Bjorkman  (pers. 
comm.).  Two-way  analysis  of  variance  (ANOVA)  Tukey-Kramer  tests  were  used  to  de¬ 
termine  significant  differences  in  DIP  concentrations  among  ELF  ratings  and  oceans  for 
stations  with  both  DIP  concentrations  and  ELF  samples  (a  =  0.05).  DIP  concentrations 
below  detection  level  (BDL)  were  considered  to  be  0  in  the  calculations. 

Nitrogenase  activity.  Nitrogenase  activity  was  assayed  during  KM0701  and  the  Atlantic 
cruises  by  the  acetylene  reduction  assay  described  by  Capone  and  Montoya  (2001).  Five 
to  twenty  colonies  were  placed  in  20-30  mL  of  0.2  //m-filtered  sea  water  in  75  mL  square 
polycarbonate  culturing  bottles  with  silicone  septa.  Two  to  three  replicate  bottles  were 
used  for  each  incubation.  Acetylene  was  produced  on  board  by  mixing  50  mL  double 
distilled  water  with  15  g  calcium  carbide  and  stored  in  a  bladder.  Six  mL  of  acetylene 
were  injected  into  the  bottles.  The  bottles  were  incubated  in  a  Percival  incubator  at  am¬ 
bient  surface  seawater  temperatures  under  soft  white  fluorescent  bulbs  at  50  //mol  quanta 
m”2  s”1  light.  Zero  times  were  targeted  at  11:00  h,  12:00  h,  and  13:00  h,  local  time.  At 
times  0,  1  h,  and  2  h  after  injection,  duplicate  0.2-0. 3  mL  samples  were  removed  from  the 
headspace  of  the  bottles  and  injected  into  a  Shimadzu  GC-8A  gas  chromatograph  and  in¬ 
tegrated  by  a  Shimadzu  CR8A  Chromatopac  to  measure  the  ethylene  peaks.  Standards  of 
9.1  fj L  L  1  ethylene  were  used  to  calibrate  the  peak  heights.  Ethylene  formed  was  cal¬ 
culated  according  to  Capone  and  Montoya  (2001)  with  Bunsen  coefficients  calculated 
from  surface  salinity  and  temperature  according  to  Breitbarth  et  al.  (2004).  Nitrogenase 
activity  was  calculated  as  the  average  rate  of  ethylene  production  per  colony  (nmol  C2H4 
h  1  colony”1).  Acetylene  reduction  rates  for  MP09  were  provided  by  D.  Capone  (pers. 
comm.).  One-way  ANOVA  Tukey-Kramer  tests  were  used  to  determine  statistically  dif¬ 
ferent  rates  between  colony  morphologies  and  among  the  cruise  means  (a  =  0.05). 
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Results 


Physical  conditions.  The  North  Atlantic  upper  ocean  was  stratified  with  late  spring  SST 
ranging  from  23°C  to  25°C  and  MLD  ranging  from  4  to  37  m  during  EN355  and  autumn 
SST  ranging  from  25°C  to  28°C  and  MLD  ranging  from  70  m  in  the  northern  part  of  the 
transect  to  26  m  in  the  south  during  EN361.  During  MP09,  the  summer  North  Pacific  up¬ 
per  ocean  was  highly  stratified  with  SST  ranging  from  26.5°C  to  27.5°C  and  MLD  rang¬ 
ing  from  25  m  to  58  m.  KM0701  cut  across  the  west  Pacific  warm  pool,  with  early  win¬ 
ter  deep  mixing  up  to  1 1 8  m  in  the  north,  early  summer  stratification  as  shallow  as  1 1  m 
in  the  south,  and  a  maximum  SST  of  over  30°C  just  south  of  the  Equator  (Fig.  5-2).  As 
the  MLD  shoaled  in  the  southern  end  of  the  west  Pacific  transect,  DIP  concentrations  de¬ 
creased  (Fig.  5-2). 
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Distance  along  transect  (km) 


Figure  5-2:  Hydrographic  conditions  of  KM0701.  The  grayscale  pseudocolor  plot  is  po¬ 
tential  temperature  (°C)  as  a  function  of  depth  (m)  vs.  distance  along  the  transect  (km) 
starting  from  Hawaii.  The  connected  black  circles  are  mixed  layer  depth  (m)  while  the 
white  circles  indicate  DIP  concentrations  (nmol  L  1 ).  The  vertical  dashed  black  line 
represents  the  Equator  and  the  vertical  solid  white  line  denotes  a  turning  point  in  the 
transect. 
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ELF.  North  Atlantic  samples  consistently  showed  more  ELF  labeling  than  Pacific  sam¬ 
ples.  A  spatial  representation  of  the  maximum  ELF  labeling  at  each  station  is  shown 
along  the  cruise  tracks  in  Fig.  5-1.  During  the  timeframe  of  this  study,  North  and  South 
Pacific  Trichodesmium  field  samples  did  not  show  high  ELF  labeling.  All  samples  were 
rated  ‘+  -’  or  ‘-  -’  ,  indicating  that  there  was  low  to  no  AP  activity  (Table  5.1;  Figs.  5-1 
A,  B;  5-3  A-H).  Labeling  of  South  Pacific  Trichodesmium  colonies  was  highest  on  the 
shelf  near  Australia  at  ‘+  -  ’  (Fig.  5-1  A).  In  contrast.  North  Atlantic  ELF  samples  ranged 
from  ‘-  -’  to  ‘  +  +’  (Table  5.1;  Figs.  5-1  C,  5-3  I-P).  Labeling  of  the  two  colony  morpholo¬ 
gies  was  not  always  at  the  same  level,  as  seen  at  Sta.  16,  18,  and  20  of  MP09  (Table  5.1, 
Fig.  5-3  A,  B)  and  Sta.  14,  26,  and  29  of  KM0701  (Table  5.1,  Fig.  5-3  G,  H).  Labeling 
of  colonies  at  the  same  station  on  different  days  could  vary,  as  seen  at  MP09  Sta.  17  and 
22  (Table  5.1,  Fig.  5-3  C,  D)  and  KM0701  Sta.  14  (Table  5.1).  Two  colonies  on  the  same 
slide  could  also  differ.  For  example,  Sta.  7  and  8  from  EN355  had  some  highly  labeled 
filaments  among  unlabeled  filaments,  giving  the  slides  a  rating  of  ‘+  -’  (Fig.  5-3  J-L).  Al¬ 
though  Trichodesmium  colonies  collected  were  identified  by  morphology,  two  to  three 
different  filament  sizes  were  observed  in  over  half  of  the  samples  (Fig.  5-3  L).  Color  from 
the  autofluorescence  of  phycoerythrin  in  the  filaments  varied,  ranging  from  the  usual  or¬ 
ange  color  (Table  5.1;  Fig.  5-3  A,  B,  E-H)  to  a  dim  green  or  blue  (Table  5.1;  Fig.  5-3  C, 

D,  J). 

In  spite  of  efforts  to  wash  the  colonies  for  the  ELF  assay,  a  variety  of  microbes 
were  found  in  association  with  Trichodesmium  colonies  including  diatoms,  heterotrophic 
bacteria,  and  the  filamentous  N2-fixing  cyanobacteria  Plectonema  sp.  and  Calothrix  sp. 
(Fig.  5-4  A-E).  As  in  Dyhrman  et  al.  (2002),  some  of  these  organisms  were  found  to  be 
ELF-labeled  even  when  Trichodesmium  was  not,  indicating  AP  activity  in  these  other 
taxa:  diatoms  (North  Pacific  MP09-21,  North  Atlantic  EN361-1;  Fig.  5-4  A-B),  het¬ 
erotrophic  bacteria  (over  half  of  the  stations  on  all  four  cruises;  Fig.  5-4  C),  and  Plec¬ 
tonema  sp.  (North  Pacific  MP09-19,  North  Atlantic  EN361-1,  2,  3;  Fig.  5-4  D).  Calothrix 
sp.  was  found  in  association  with  Trichodesmium  near  Hawaii  during  MP09  and  was  not 
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Table  5.1:  Data  summary  for  stations  with  ELF  data.  For  stations  with  multiple  colony 
morphologies,  data  are  separated  by  7’ .  NA  =  not  available. 


Date 

Station 

Latitute 

Longitude 

Colony 

ELF 

Color 

01 Jun  01 

EN355-1 

37°  6'N 

71°  42’W 

puff/tuft 

dim/dim 

02  Jun  01 

EN355-2 

35°  52'N 

71°  43’W 

puff 

+  - 

dim 

03  Jun  01 

EN355-3 

33°  0'N 

71°  40’W 

puff 

+  + 

dim 

04  Jun  01 

EN355-4 

31°  33’N 

71°  39’W 

puff/tuft 

+  +/+- 

dim/dim 

05  Jun  01 

EN355-5 

30°  29'N 

70°  0.0’W 

puff/bowtie 

NA/NA 

NA/NA 

07  Jun  01 

EN355-6 

32°  2.4'N 

70°  19.51'W 

NA 

+  + 

dim-normal 

09  Jun  01 

EN355-7 

33°  16.67’N 

68°  59.32'W 

NA 

+  - 

dim-normal 

1 1  Jun  01 

EN355-8 

37°  27.32’N 

69°  45.28’W 

NA 

+  - 

dim-normal 

04  Oct  01 

EN361-1 

36°  31.53’N 

72°  26.78'W 

puff/tuft 

+  -/+- 

dim-normal/normal 

05  Oct  01 

EN361-2 

33°  44.52’N 

70°  59.79'W 

puff/tuft 

dim-norma/dim-normall 

05  Oct  01 

EN361-3 

32°  40.89’N 

70°  9.3'W 

puff/tuft 

+  -/+- 

dim-normal/dim-normal 

06  Oct  01 

EN361-4 

30°  20.35’N 

68°  22.25'W 

puff/tuft/bowtie 

+  -/+  -/+  - 

dim-normal/ dim-normal/ dim-normal 

06  Oct  01 

EN361-5 

30°  6.42’N 

67°  32.88'W 

puff/tuft 

+  -/+- 

normal/dim-normal 

07  Oct  01 

EN361-6 

29°  10.21’N 

64°  43.528’W 

puff/tuft 

+  -/+- 

dim-normal/ dim-normal 

06  Aug  03 

MP09-13 

20°  29.59’N 

157°  41.32’W 

puff/tuft 

normal/normal 

07  Aug  03 

MP09-14 

18°  30.00’N 

157°  00.00’W 

puff 

normal 

10  Aug  03 

MP09-16 

19°  32.49’N 

156°  55.74’W 

puff/tuft 

normal/normal 

1 1  Aug  03 

MP09-16 

19°  30.55’N 

156°  55.74’W 

puff/tuft 

+ 

normal/normal 

12  Aug  03 

MP09-17 

19°  31.23’N 

158°  58.75’W 

puff/tuft 

normal/normal 

13  Aug  03 

MP09-17 

19°  30.10’N 

158°  59.94’W 

puff 

+  - 

normal 

14  Aug  03 

MP09-18 

20°  31.84’N 

160°  59.81’W 

puff/tuft 

dim/dim 

15  Aug  03 

MP09-19 

21°  00.79’N 

158°  59.70’W 

puff/tuft 

dim/normal 

16  Aug  03 

MP09-19 

21°  00.18’N 

158°  59.69’W 

puff/tuft 

dim/normal 

17  Aug  03 

MP09-20 

20°  16.09’N 

156°  31.72’W 

puff/tuft 

normal/normal 

18  Aug  03 

MP09-21 

20°  15.61’N 

159°  11.21’W 

puff/tuft 

+  -/+- 

normal/dim 

19  Aug  03 

MP09-22 

19°  06.67’N 

162°  13.05’W 

puff/tuft 

+  -/+- 

normal/dim 

20  Aug  03 

MP09-22 

19°  06.29’N 

162°  13.60’W 

puff/tuft 

--/-- 

dim/dim 

21  Aug  03 

MP09-23 

20°  30.68’N 

161°  30.27’W 

puff/tuft 

normal/dim 

07  Jan  07 

KM0701-5 

12°  26.11’N 

167°  43.75’W 

puff/tuft 

normal/dim-normal 

10  Jan  07 

KM0701-8 

5°  38.86'N 

174°  32.18’W 

tuft 

dim-normal 

1 1  Jan  07 

KM0701-9 

3°  14.40'N 

176°  53.04’W 

puff 

normal 

12  Jan  07 

KM0701-10 

0°  22.02'N 

179°  38.60’W 

tuft 

normal 

15  Jan  07 

KM0701-12 

4°  43.34'S 

174°  43.85’E 

puff 

normal 

16  Jan  07 

KM0701-13 

7°  3.64’S 

172°  18.98’E 

tuft 

normal 

17  Jan  07 

KM0701-14 

9°  15.03’S 

169°  59.96’E 

puff/tuft/bowtie 

.  ./-  -/-  - 

dim-normal/normal/normal 

19  Jan  07 

KM0701-14 

9°  15.01'S 

170°  0.01’E 

puff/tuft 

normal/dim-normal 

20  Jan  07 

KM0701-15 

12°  34.55’S 

169°  51.55’E 

tuft 

normal 

21  Jan  07 

KM0701-16 

15°  53.63’S 

169°  43.02’E 

puff 

normal 

22  Jan  07 

KM0701-17 

19°  13.95’S 

169°  34.63’E 

puff/tuft 

+  -/+- 

dim-normal/normal 

23  Jan  07 

KM0701-19 

21°  37.42’S 

168°  39.49’E 

puff/tuft 

normal/normal 

26  Jan  07 

KM0701-20 

25°  40.30’S 

165°  25.03’E 

NA/puff/tuft 

-  -/NA/NA 

normal/NA/NA 

27  Jan  07 

KM0701-21 

29°  2.43'S 

164°  20.26’E 

puff/tuft 

--/+- 

normal/normal 

03  Feb  07 

KM0701-26 

32°  25.29’S 

159°  5.33’E 

puff/tuft 

--/+- 

normal/dim-normal 

07  Feb  07 

KM0701-28 

30°  15.51’S 

157°  18.19’E 

puff/tuft 

+  -/+- 

normal/dim-normal 

08  Feb  07 

KM0701-29 

29°  45.56’S 

156°  37.43’E 

puff/tuft 

--/+- 

normal/normal 

09  Feb  07 

KM0701-30 

28°  45.68’S 

155°  22.16’E 

puff/tuft 

+  -/NA 

normal/NA 
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Figure  5-3:  Selected  digital  micrographs  of  ELF-assayed  Trichodesmium  from  the  four 
cruises:  (A-D)  MP09  during  August  2003,  in  the  North  Pacific;  (E-H)  KM0701  during 
January  and  February  2007,  across  the  west  Pacific  warm  pool;  and  (I-L)  EN355  during 
June  2001,  and  (M-P)  EN361  during  October  2001,  in  the  western  North  Atlantic.  White 
scale  bars  correspond  to  10  /j m.  (A)  Sta.  MP09-16,  10  Aug  03,  tuft,  (B)  Sta.  MP09- 
16,  10  Aug  03,  puff,  ‘+  (C)  Sta.  MP09-22,  20  Aug  03,  puff,  (D)  Sta.  MP09-22, 

19  Aug  03,  tuft,  ‘+  (E)  Sta.  KM0701-14,  17  Jan  07,  puff,  -  ’;  (F)  Sta.  KM0701-17, 

22  Jan  07,  puff,  ‘+  (G)Sta.  KM0701-26,  03  Feb  07,  puff,  (H)  Sta.  KM0701-26, 

03  Feb  07,  tuft,  ‘+  (I)  Sta.  6-6,  07  Jun  01,  ‘+  +  ’;  (J)  Sta.  6-7,  09  Jun  01,  ‘+  -’;  (K)  Sta. 

6-7,  09  Jun  01,  ‘+-’;(F)  Sta.  6-8,  11  Jun  01,  ‘+-’;(M)  Sta.  10-1,  04  Oct  01,  puff, 

(N)  Sta.  10-2,  05  Oct  01,  tuft,  ‘+  +’;  (O)  Sta.  10-3,  05  Oct  01,  puff,  ‘+  -’;  (P)  Sta.  10-5,  06 
Oct  01,  puff,  ‘+ 
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Figure  5-4:  Organisms  associated  with 
Trichodesmium  colonies.  White  arrows 
indicate  the  organisms  and  white  scale 
bars  correspond  to  10  /rm.  (A)  ELF- 
labeled  pennate  diatoms  (Sta.  MP09- 
17);  (B)  the  same  pennate  diatoms  in 
(A)  viewed  under  a  chlorophyll  filter; 
(C)  ELF-labeled  heterotrophic  bacteria 
(Sta.  MP09-16,  11  Aug  03,  puff);  (D) 
ELF-labeled  cyanobacterium  Plec- 
tonema  sp.  (Sta.  MP09-19,  15  Aug  03, 
puff);  (E)  INh-fixing  cyanobacterium 
Calothrix  sp.  (Sta.  MP09-23,  21  Aug 
03,  tuft). 


ELF-labeled  (Sta.  13,  17,  18,  19,  21,  23;  Fig.  5-4E). 


DIP  and  ELF.  ELF  labeling  varied  with  ocean  basin  and  corresponded  negatively  with 
DIP  concentrations.  North  Atlantic  stations  had  significantly  lower  DIP  concentrations 
than  the  North  and  South  Pacific  stations  in  a  2-way  ANOVA  ( p  =  0.0001)  (Table  5.2). 
DIP  concentrations  ranged  from  below  detection  level  (BDL)  to  25  nmol  L  1  in  the  North 
Atlantic,  30  to  320  nmol  L-1  in  the  North  Pacific,  and  50  to  220  nmol  L”1  in  the  South 
Pacific.  Stacked  histograms  of  DIP  concentrations  were  plotted  for  each  ELF  label  rat¬ 
ing  (Fig.  5-5).  Trichodesmium  colonies  with  little  or  no  ELF  labeling  were  found  in  sta¬ 
tions  from  all  cruises,  spanning  the  full  range  of  DIP  concentrations  (Fig.  5-5  A).  Sam¬ 
ples  with  moderate  levels  of  labeling  spanned  a  moderate  range  of  DIP:  BDL  -  169  nmol 
L_1  (Fig.  5-5  B).  High  levels  of  labeling  were  found  only  in  samples  from  North  Atlantic 
stations  with  DIP  concentrations  less  than  or  equal  to  40  nmol  L_1  (Fig.  5-5  C).  In  a  2- 
way  ANOVA,  DIP  concentrations  for  stations  with  an  ELF  rating  of  were  greater 
than  DIP  concentrations  for  stations  with  an  ELF  rating  of  ‘+  -’  or  ’+  +’  ,  significant  at  the 
94%  confidence  level  ( p  =  0.0554)  (Table  5.2). 
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Table  5.2:  2-way  ANOVA  on  DIP  concentrations  among  ELF  ratings  and  ocean  basins 


Factor 

F 

P 

Groupings 

ELF  rating 

3.04 

0.0554 

‘++’  and  ‘+-’; 

Ocean  basin 

11.55 

0.0001 

North  Atlantic;  North  and  South  Pacific 

DIP  (nmol  L  ') 


Figure  5-5:  Histogram  of  binned  surface  DIP  concentrations  (nmol  L-1)  separated  by 
ELF  label  ratings:  (A)  ‘-  (B)  ‘+  (C)  ‘+  +.’  Colors  represent  ocean  basin:  North 

Atlantic  (black).  North  Pacific  (gray),  and  South  Pacific  (white).  Bin  width  is  25  nmol 
L-1. 


Nitrogenase  activity.  Nitrogenase  activity  varied  significantly  between  the  basins  (Table 
5.3,  Fig.  5-6).  In  the  P-stressed  western  North  Atlantic,  mean  nitrogenase  activities  mea¬ 
sured  during  EN355  and  EN361  were  significantly  higher  than  rates  in  both  the  North  and 
South  Pacific  ( p  =  5.35  x  10-8):  0.80  ±  0.22  and  0.79  ±  0.14  nmol  C2H4  h-1  colony-1 
compared  to  0.26  ±  0.03  nmol  C2H4  h-1  colony-1  (D.  Capone  pers.  comm.)  and  0.09  ± 
0.02  nmol  C2H4  h-1  colony-1,  respectively  (Fig.  5-6).  For  stations  with  nitrogenase  ac- 
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tivity  measurements  for  both  puff  and  tuft  colonies,  rates  were  generally  greater  for  tuft 
colonies  than  puff  colonies  although  the  rates  were  not  significantly  different  between 
colony  morphs  according  to  one-way  ANOVA  Tukey-Kramer  tests  (p  >  0.1355).  There 
were  two  exceptions  to  this  pattern.  At  Sta.  KM0701-14,  puff  nitrogenase  activity  was 
greater  than  that  of  tufts  but  not  significantly  ( p  =  0.2379),  and  at  Sta.  KM0701-30  near 
Australia,  the  rate  for  tuft  colonies  was  significantly  greater  than  that  of  puff  colonies 

(p  =  0.0022). 
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Station 


Figure  5-6:  Bar  plot  of  nitrogenase  activity  (nmol  C2H4  produced  colony-1!!"1)  for  each 
station.  Error  bars  show  standard  error  of  the  mean.  Data  is  separated  by  puffs  (black), 
tufts  (gray),  or  mixed  puffs  and  tufts  (white).  The  North  Pacific  is  represented  by  the  aver¬ 
age  rates  for  MP09  Stas.  14-22  (D.  Capone,  pers.  comm.).  Ocean  basins  are  separated  by 
vertical  dotted  lines. 

Discussion 

Trichodesmium  is  a  source  of  new  N  in  tropical  and  subtropical  oceans  around  the  world, 
but  the  factors  limiting  N2  fixation  in  Trichodesmium  are  not  completely  understood.  The 
South  Pacific  in  particular  is  an  undersampled  region  (Campbell  et  al.,  2005).  Our  study 
investigated  P  stress  and  N2  fixation  in  Trichodesmium  in  three  major  ocean  regions:  the 
western  North  Atlantic,  the  North  Pacific,  and  the  western  South  Pacific.  Our  results  show 
patterns  of  both  higher  P  stress  and  higher  N2  fixation  rates  in  the  western  North  Atlantic 
than  in  the  Pacific. 

During  the  2001  western  North  Atlantic  cruises,  colonies  of  Trichodesmium  showed 
more  frequent  and  more  complete  ELF  labeling  than  the  Pacific  samples  in  this  study. 

This  difference  in  ELF  labeling  is  a  qualitative  indication  of  higher  AP  activity  in  the 
North  Atlantic.  The  samples  near  Australia  also  showed  more  labeling  than  other  regions 
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of  the  South  Pacific.  In  this  region  the  MLD  was  shallower  and  the  DIP  concentrations 
were  lower.  Trichodesmium  thrives  in  warm,  stratified  waters,  and  blooms  are  found  in 
oligotrophic  waters  where  the  sea  surface  temperature  (SST)  >  25°C  (Karl  et  al.,  2002). 
However,  supply  of  P  to  surface  waters  can  be  cut  off  by  warming  and  stratification  of 
surface  layers.  ELF  is  not  quantitative,  so  finding  a  direct  parametric  correlation  with  DIP 
concentrations  is  not  possible.  However,  by  examining  the  range  of  DIP  concentrations 
associated  with  the  different  levels  of  ELF  label,  it  is  apparent  that  higher  labeling  is  as¬ 
sociated  with  lower  DIP  concentrations.  This  is  in  contrast  to  Webb  et  al.  (2007),  where 
they  found  no  significant  relationship  between  ELF  labeling  and  DIP  during  a  transect  in 
the  western  Equatorial  to  western  South  Atlantic.  In  our  study,  we  had  more  data  points  of 
ELF  and  a  broader  range  of  DIP  than  (Webb  et  al.,  2007)  (5-320  nmol  L”1  compared  to 
10-120  nmol  L^1),  and  this  breadth  of  data  may  have  enabled  us  to  find  a  significant  re¬ 
lationship  between  DIP  and  ELF  labeling.  Also  in  contrast  to  (Webb  et  al.,  2007),  we  did 
not  see  consistent  differences  in  labeling  between  puff  and  tuft  colonies.  In  a  study  of  the 
North  and  South  Atlantic  subtropical  gyres,  surface  water  column  AP  activity  was  higher 
and  DIP  was  lower  in  the  North  than  the  South  Atlantic  (Mather  et  al.,  2008).  AP  activity 
in  the  NPSG  measured  during  cruise  MP09  was  4.4  times  lower  than  AP  activity  in  the 
North  Atlantic  (Sohm  et  al.,  2008)  which  is  consistent  with  our  ELF-based  observations. 
The  difference  in  AP  activity  observed  by  ELF  labeling  and  hydrolytic  rates  among  these 
regions  suggests  that  the  Atlantic  populations  were  generally  more  P-stressed  than  those 
in  the  Pacific. 

In  addition  to  AP  activity,  indicators  used  to  determine  P  stress  or  potential  lim- 

o _ 

itation  include  DIP  concentrations,  maximum  uptake  rate  of  PO4  (V,MY),  cellular  N:P 
ratios,  and  turnover  times  of  DIP.  Results  from  these  metrics  used  in  the  tropical  and  sub¬ 
tropical  Atlantic  and  Pacific  Oceans  are  consistent  with  our  observations  of  AP  activity. 
The  critical  concentration  of  DIP  for  Trichodesmium  growth  has  been  determined  to  be  9 
nmol  L  1  (Moutin  et  al.,  2005).  While  there  is  some  overlap,  Atlantic  DIP  concentrations 
are  typically  lower  than  this  critical  value  while  Pacific  concentrations  are  typically  higher 
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(Table  5.3)  (Wu  et  al.,  2000).  The  \max  of  PO4  uptake  by  Trichodesmium  colonies  mea¬ 
sured  during  a  cruise  in  the  North  Atlantic  was  four  times  more  than  that  of  the  North 
Pacific  during  cruise  MP09  (Sohm  et  al.,  2008).  Concentrations  of  DIP  in  the  Sargasso 
Sea  region  of  the  Atlantic  are  1-2  orders  of  magnitude  lower  than  those  in  the  subtrop¬ 
ical  North  Pacific  (Wu  et  al.,  2000).  Molar  ratios  of  cellular  N:P  may  indicate  nutrient 
limitation  and/or  storage.  The  reference  Redfield  ratio  of  N:P  is  16.  The  P  limitation  of 
growth  in  Trichodesmium  begins  at  N:P  of  40-50  (Kustka  et  al.,  2003),  and  in  P-restricted 
cultures,  N:P  can  reach  >  90  (White  et  al.,  2006a).  Trichodesmium  N:P  ratios  averaged 
40.1  ±  2.53  during  the  North  Pacific  MP09  cruise  while  they  ranged  from  60  in  the  west 
to  30  in  the  east  in  the  subtropical  North  Atlantic  in  2000,  indicating  more  P  limitation  in 
the  western  North  Atlantic  (Krauk  et  al.,  2006).  Turnover  time  is  calculated  as  the  inven¬ 
tory  of  a  nutrient  divided  by  its  uptake  rate  and  can  be  considered  as  the  period  it  would 
take  to  deplete  a  nutrient  if  input  was  halted  (Ammerman  et  al.,  2003).  Short  turnover 
times  of  DIP  (Xdip)  result  from  low  inventory  and/or  high  uptake  rates  and  can  indicate 
low  DIP  availability.  At  the  Climax  station  in  the  North  Pacific,  Tqip  was  48  h  in  the  sum¬ 
mer  (Bjorkman  et  al.,  2000).  In  contrast,  Sargasso  Sea  x dip  was  measued  to  be  around  9 
h  (Cotner  et  al.,  1997).  Turnover  times  of  DIP  in  the  South  Pacific  near  New  Caledonia 
varied  widely  and  decreased  from  468  h  in  the  winter  to  4  h  in  the  early  fall,  indicating 
P  depletion  after  the  stratification  of  the  water  column  in  the  summer  (Van  Den  Broeck 
et  al.,  2004).  These  P  stress  indicators  are  consistent  with  higher  potential  for  P  limita¬ 
tion  in  the  subtropical  western  North  Atlantic  than  the  subtropical  North  and  South  Pacific 
with  a  region  of  summer  P  limitation  in  the  South  Pacific  near  New  Caledonia. 

Colonies  used  in  the  ELF  assay  were  chosen  in  an  attempt  to  represent  the  com¬ 
munity,  but  a  number  of  issues  arise  because  of  the  small  sample  size.  Field  populations 
of  Trichodesmium  have  a  range  of  characteristics  including  nutrient  history,  species,  and 
degree  of  senescence  which  might  not  be  reflected  in  two  colonies  per  morphology  per 
station.  ELF  labeling  varies  within  a  sample,  a  colony,  or  even  along  a  filament.  While 
colonies  of  the  same  morphology  were  targeted,  filaments  observed  on  the  ELF  slides 
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were  from  a  variety  of  Trichodesmium  species  based  on  cell  width  and  length.  Samples 
also  displayed  a  range  of  phycobiliprotein  autofluorescence.  The  dim  filaments  observed 
may  be  due  to  senescent  colonies,  bleaching,  or  fading.  In  some  cases,  labeling  on  a  slide 
was  higher  on  dim  filaments  than  normal  filaments,  which  may  indicate  elevated  AP  ac¬ 
tivity  in  senescent  colonies  or  in  heterotrophic  bacteria  associated  with  these  colonies  in 
some  cases.  While  the  colonies  picked  for  the  ELF  assay  may  not  capture  the  diversity  of 
field  populations,  they  do  give  useful  snapshots  of  the  state  of  the  community. 

A  variety  of  organisms  including  bacteria,  fungi,  diatoms,  dinoflagellates,  chrys- 
ophytes,  ciliates,  amoebae,  hydroids,  and  crustacean  larvae  can  be  enriched  by  2-5  orders 
of  magnitude  in  association  with  Trichodesmium  colonies  (Sheridan  et  al.,  2002).  Hy¬ 
drolysis  rates  measured  on  Trichodesmium  colonies  may  reflect  the  activity  of  the  whole 
community  rather  than  Trichodesmium- specific  activity  (Dyhrman  et  al.,  2002;  Webb 
et  al.,  2007)).  The  activity  of  AP  from  associated  organisms  such  as  pennate  diatoms, 
heterotrophic  bacteria,  and  Plectonema  may  result  in  an  indirect  source  of  DIP  that  can  be 
utilized  by  Trichodesmium.  The  presence  of  AP  in  these  microorganisms  may  also  result 
in  the  overestimation  of  Trichodesmium  AP  activity  measured  through  bulk  assays.  With 
ELF,  AP  activity  can  be  visually  localized  to  specific  cells. 

In  addition  to  AP  measurements,  N2  fixation  rates  attributed  to  Trichodesmium 
may  be  influenced  by  the  diazotrophs  tentatively  identified  as  Plectonema  sp.  and  Calothrix 
sp.  found  in  close  association  with  Trichodesmium  colonies.  Both  of  these  cyanobacte¬ 
ria  occupy  a  wide  variety  of  habitats,  including  marine.  The  long,  thin  filaments  of  non- 
heterocystous  Plectonema  are  found  in  Trichodesmium  colonies  around  the  world  (Sid- 
diqui  et  al.,  1992).  Calothrix  is  identified  by  its  terminal  heterocyst  and  tapering  trichome 
(Waterbury,  2006),  and  is  found  in  the  upper  intertidal  zone  (Whitton  and  Potts,  1982)  and 
as  an  epiphyte  of  the  pelagic  diatom  Chaetoceros  (Foster  and  Zehr,  2006).  While  Plec¬ 
tonema  is  known  to  fix  N2  only  in  lowered  oxygen  tension  (Rippka  et  al.,  1979),  both 
Plectonema  and  Calothrix  may  be  contributing  to  N2  fixation  rates  measured  from  Tri¬ 
chodesmium  colonies. 
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During  the  four  cruises  analyzed  in  this  study,  Trichodesmium  colony-based  N2 
fixation  rates  were  significantly  higher  in  the  subtropical  North  Atlantic  than  in  the  tropi¬ 
cal  and  sub-tropical  North  and  South  Pacific.  North  Atlantic  values  averaged  0.80  ±  0.12 
nmol  C2H4  h'1  colony'1  while  Pacific  values  averaged  0.13  ±  0.02  nmol  C2H4  h'1 
colony'1.  As  in  Webb  et  al.  (2007),  tuft  colonies  generally  had  higher  per  colony  rates 
than  puffs,  however  the  difference  between  colony  morphologies  were  not  significant  in 
our  study  with  one  exception.  For  comparison,  literature  values  of  N2  fixation  rates  mea¬ 
sured  by  acetylene  reduction  averaged  87  ±  14  mmol  N  m'2yr'!  in  the  western  North 
Atlantic  from  1994-2003  (1.6  x  1012  mol  N  yr'1  over  17.8  x  106  km2;  (Capone  et  al., 
2005)),  while  the  North  Pacific  averaged  31  ±  18  mmol  N  m'2  yr”1  from  1990-1992 
(Karl  et  al.,  1997).  Biogeochemical  estimates  of  the  annual  input  of  new  N  by  N2  fixation 
were  2.0  x  1012  mol  N  yr” 1 ,  or  72  mmol  N  m  2  yr”1  in  the  subtropical  North  Atlantic 
(Gruber  and  Sarmiento,  1997)  and  4±  1  x  1012  mol  N  yr”1,  or  39  ±  9  mmol  N  m”2  yr”1 
in  the  tropical  and  subtropical  Pacific  (Deutsch  et  al.,  2001).  In  a  global  biogeochemical 
elemental  cycling  ocean  model,  N2  fixation  rates  averaged  18  mmol  N  m  2  yr  1  in  the 
tropical  and  subtropical  Atlanic  and  11  mmol  N  m  2  yr'1  in  the  tropical  and  subtropical 
Pacific  (Moore  and  Doney,  2007).  Our  nitrogenase  activity  measurements  are  in  congru¬ 
ence  with  the  pattern  found  in  these  estimates;  N2  fixation  rates  are  higher  in  the  tropical 
and  subtropical  North  Atlantic  than  in  the  tropical  and  subtropical  Pacific. 

The  difference  in  N2  fixation  rates  between  the  Atlantic  and  Pacific  Oceans  may 
be  due  to  disparate  nutrient  limitation  regimes.  N2  fixation  in  the  western  North  Atlantic 
is  more  P-limited  than  Fe-limited  (Wu  et  al.,  2000;  Sanudo-Wilhelmy  et  al.,  2001;  Moore 
et  al.,  2004).  When  Trichodesmium  cells  have  met  their  Fe  cell  quota,  N2  fixation  cor¬ 
relates  with  cellular  P  quotas  (Sanudo-Wilhelmy  et  al.,  2001).  In  the  western  Equato¬ 
rial  and  South  Atlantic,  N2  fixation  in  shipboard  incubations  was  higher  when  P  stress 
was  relieved  by  the  addition  of  DIP  than  in  control  or  Fe  addition  incubations  (Webb 
et  al.,  2007).  Trichodesmium  AP  activity  is  higher  in  the  North  Atlantic  than  in  the  Pacific 
Ocean  (Table  5.3)  (Mulholland  et  al.,  2002;  Sohm  and  Capone,  2006;  Sohm  et  al.,  2008), 
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and  this  difference  is  reflected  by  the  visualization  of  AP  activity  through  ELF  (Table  5.3, 
Figs.  5-1,  5-3). 

In  contrast,  N2  fixation  in  the  North  Pacific  may  be  limited  by  factors  other  than 
P.  The  North  Pacific  receives  less  Fe-rich  dust  input  than  the  North  Atlantic  (Fung  et  al., 
2000),  and  Fe  is  necessary  for  Trichodesmium  growth  and  the  enzyme  nitrogenase  used 
in  N2  fixation.  In  the  South  Pacific,  regions  near  New  Caledonia  and  Australia  receive  Fe 
from  continental  margins  and  Australian  dust  (Moore  and  Braucher,  2008).  Fe  stimulates 
the  growth  of  diazotrophs,  drawing  down  P  (Karl,  2002;  Mather  et  al.,  2008),  so  Fe-rich 
areas  such  as  the  western  North  Atlantic  and  western  South  Pacific  near  Australia  are 
more  likely  to  be  P-limited.  In  a  coupled  biogeochemistry  ecosystem  circulation  model, 
diazotrophs  were  limited  by  P  in  the  Atlantic  and  by  Fe  in  the  Pacific,  and  increasing 
Fe  inputs  in  the  model  increased  N2  fixation  in  the  tropics  and  subtropics  (Moore  et  al., 
2004).  These  data  suggest  that  N2  fixation  rates  of  Trichodesmium  in  the  western  North 
Atlantic  and  in  the  South  Pacific  near  Australia  are  limited  by  P  while  the  Pacific  subtrop¬ 
ical  gyres  may  be  limited  by  Fe  or  co-limited  by  Fe  and  P. 

While  P  limitation  may  be  alleviated  by  the  use  of  DOP  by  Atlantic  Trichodesmium, 
Fe  limitation  in  Pacific  populations  cannot  be  alleviated  by  the  use  of  an  alternative  pool 
(Webb  et  al.,  2007;  Mather  et  al.,  2008).  In  a  comparison  of  the  North  and  South  Atlantic, 
DIP  released  from  AP  activity  in  the  P-stressed  North  Atlantic  supported  up  to  30%  of 
primary  production  and  carbon  fixation  rates  there  were  comparable  to  rates  in  the  P- 
replete  South  Atlantic  (Mather  et  al.,  2008).  Fe  may  be  the  proximal  limiting  nutrient  ac¬ 
cording  to  Leibig’s  law  while  P  may  be  the  ultimate  limiting  nutrient,  controlling  produc¬ 
tivity  over  long  timescales  (Tyrrell,  1999).  Additional  assessments  of  Trichodesmium  P 
and  Fe  physiology  in  open  ocean  gyres  are  necessary  to  elucidate  nutrient  controls  on  N2 
fixation,  but  our  study  supports  previous  observations  that  the  North  Atlantic  is  P-limited 
while  the  Pacific  is  likely  to  be  Fe-limited.  This  study  underscores  the  contrasting  impor¬ 
tance  of  P  (and  likely  Fe)  to  the  activity  of  the  key  diazotroph  Trichodesmium  in  different 
ocean  basins. 
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CHAPTER  6 


Concluding  remarks 


Thesis  summary 

My  thesis  investigated  the  potential  for  niche  differentiation  in  the  filamentous  diazotroph 
Trichodesmium.  In  Chapt.  2, 1  found  that  there  were  two  clades  of  Trichodesmium  among 
our  cultured  strains.  Clade  I  was  comprised  of  Trichodesmium  erythraeum  and  Trichodesmium 
contortum,  and  clade  II  consisted  of  Trichodesmium  thiebautii,  Trichodesmium  tenue,  Tri¬ 
chodesmium  hildebrandtii,  and  Trichodesmium  pelagicum.  There  was  a  surprising  amount 
of  diversity  in  cell  morphology  within  each  clade.  The  makeup  of  accessory  pigments 
(phycobiliproteins)  within  each  clade  was  consistent,  indicating  that  Trichodesmium  niches 
may  be  partially  defined  by  light  quality.  Using  sequences  of  hetR  from  Chapt.  2, 1  devel¬ 
oped  a  real-time  quantitative  polymerase  chain  reaction  (qPCR)  protocol  that  was  capa¬ 
ble  of  distinguishing  between  the  two  clades  in  mixture  (Chapt.  3).  The  qPCR  protocol 
was  applied  to  depth  profiles  across  transects  of  the  Equatorial  Atlantic  Ocean,  the  west 
Pacific  warm  pool,  and  the  South  Pacific  Ocean  (Chapt.  4).  I  found  that  concentrations 
of  clade  II  were  generally  higher  than  those  of  clade  I,  and  both  clades  reached  higher 
concentration  in  the  Atlantic  than  in  the  Pacific.  Concentrations  of  the  two  clades  were 
correlated  in  the  Pacific  Ocean,  indicating  that  similar  factors  were  controlling  their  dis¬ 
tribution.  However,  the  two  clades  were  not  correlated  in  the  Altantic  Ocean,  indicating 
that  different  factors  were  limiting  the  distribution  of  each  clade.  Clade  I  had  deeper  and 
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cooler  ranges  than  clade  II,  suggesting  that  light  and/or  temperature  may  define  the  niches 
of  different  types  of  Trichodesmium.  Concentrations  of  Trichodesmium  did  not  corre¬ 
late  with  dissolved  inorganic  phosphorus  (DIP)  or  iron  (Fe)  concentrations,  but  ratios  of 
the  clades  (clade  ILclade  I)  were  consistently  high  in  stations  with  low  dissolved  DIP, 
notably  the  Atlantic.  Differential  competition  for  resources  may  allow  the  two  clades  of 
Trichodesmium  to  coexist,  but  more  studies  are  required  to  understand  how  light  and  nu¬ 
trients  limit  the  growth  of  different  species  of  Trichodesmium.  In  Chapt.  5, 1  investigated 
phosphorus  (P)  stress  and  No  fixation  in  populations  in  the  western  North  Atlantic,  North 
Pacific,  and  the  west  Pacific  warm  pool  (Hynes  et  al.,  2009).  Atlantic  populations  had 
more  P  stress  and  higher  N2  fixation  rates  than  Pacific  populations,  indicating  that  factors 
other  than  P  were  limiting  Trichodesmium  N2  fixation. 

Evolution  of  Trichodesmium 

There  is  diversity  within  the  genus  Trichodesmium ,  but  where  does  this  diversity  come 
from?  There  are  two  major  clades  of  Trichodesmium ,  even  if  one  includes  the  field  sam¬ 
ples  of  T.  contortum,  T.  tenue,  and  Trichodesmium  aureum  not  represented  in  culture  (Jan- 
son  et  al.,  1999;  Lundgren  et  al.,  2005).  These  two  clades  are  circumtropical,  existing 
in  tropical  and  subtropical  waters  around  the  globe.  This  implies  that  there  was  a  diver¬ 
gence  early  in  the  evolution  of  Trichodesmium  before  the  spread  of  this  genus  around  the 
globe.  Within  each  clade,  species  are  similar  genetically  but  have  vastly  different  phe¬ 
notypes  (Chapt.  2).  Consider  clade  I;  T.  erythraeum  has  mostly  square  cells  about  7  n m 
wide  while  T.  contortum  has  disc-shaped  cells  about  25  qm  wide  and  3  /.mi  long.  Given 
the  high  genetic  similarity  of  these  two  species  based  on  16S,  ITS,  and  hetR,  their  dis¬ 
parate  phenotypes  may  be  due  to  differences  in  gene  regulation  or  post-translational  con¬ 
trols. 

Nitrogen  fixation  is  probably  the  most  important  ecological  role  for  Trichodesmium. 
Based  on  ni/H  sequences,  Trichodesmium  diverged  from  other  diazotrophs  early  in  evo¬ 
lution  (Capone  et  al.,  1997).  While  Trichodesmium  diverged  from  other  groups,  niJH 
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sequences  among  species  of  Trichodesmium  are  highly  similar  (Lundgren  et  al.,  2001), 
conserving  the  structure  of  the  N2 -fixing  enzyme  nitrogenase.  Regulation  of  N2  fixa¬ 
tion  in  Trichodesmium  may  be  carried  out  in  part  by  hetR  (El-Shehawy  et  al.,  2003).  Pu¬ 
rifying  selection  on  hetR  was  stronger  for  T.  erythraeum  than  for  other  lineages  of  Tri¬ 
chodesmium  (Mes  and  Stal,  2005).  While  the  structure  of  nitrogenase  in  Trichodesmium 
is  conserved,  its  regulation  might  vary  between  clades. 

The  copepod  Macrosetella  gracilis  is  one  of  the  few  known  grazers  of  Trichodesmium 
(O’Neil  and  Roman,  1994;  O’Neil,  1998),  but  grazing  pressure  is  considered  to  be  low 
(Carpenter  et  al.,  1987;  Eberl  and  Carpenter,  2007).  T.  thiebautii  is  known  to  produce  a 
neurotoxin,  which  has  not  been  reported  in  T.  erythraeum,  and  copepods  graze  prefer¬ 
entially  on  T.  erythraeum  (O’Neil  and  Roman,  1994).  Differential  grazing  pressure  may 
have  led  to  the  production  of  neurotoxins  in  clade  II,  but  not  clade  I. 

In  Chapt.  4, 1  found  that  clade  II  was  numerically  more  dominant  in  both  the  Pa¬ 
cific  and  Atlantic  Oceans,  but  clade  I  had  a  deeper  range  than  clade  II  in  field  popula¬ 
tions.  The  dominance  of  clade  II  and  the  wider  physical  range  of  clade  I  are  reminiscent 
of  the  classic  intertidal  study  of  barnacles  by  Connell  (1961).  Cthamalus  stellatus  was 
out-competed  by  Semibalanus  balanoides  in  lower  reaches  of  the  intertidal  zone,  but  C. 
stellatus  was  better  able  to  withstand  dessication  in  the  upper  intertidal  zone  than  S.  bal¬ 
anoides.  Two  characteristics  that  could  affect  the  ability  of  Trichodesmium  spp.  to  thrive 
at  depth  include  pigment  composition  and  gas  vesicle  strength.  In  Chapt  2, 1  found  that 
clade  I  strains  had  high  phycoerythrobilin  (PEB),  which  absorbs  light  with  wavelengths  of 
about  545  nm  (green)  and  clade  II  strains  had  high  phycourobilin  (PUB),  which  absorbs 
at  495  nm  (blue-green)  (Chapt.  2).  Blue  wavelengths  penetrate  deeper  in  the  water  col¬ 
umn  than  green  wavelengths,  so  it  would  seem  that  PUB,  and  hence  clade  II,  would  be 
better  adapted  for  life  at  depth.  Unexpectedly,  we  found  clade  I  at  depth.  T.  erythraeum 
cultures  have  been  shown  to  alter  pigment  ratios  in  response  to  light  quality  and  quantity 
(Bell  and  Fu,  2005).  Deep  clade  I  populations  might  have  higher  PUB:PEB  ratios  than 
what  we  measured  in  cullture.  Subramaniam  et  al.  (1999)  has  suggested  that  high  PUB 
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in  T.  thiebautii  field  populations  at  midday  might  be  a  mechanism  to  reduce  photooxida- 
tive  stress.  The  relationship  between  pigmentation  of  Trichodesmium  and  light  utilization 
is  a  subject  of  further  study.  Gas  vesicle  strength  is  also  variable  among  species  of  Tri¬ 
chodesmium.  Gas  vesicles  in  Trichodesmium  spp.  are  much  stronger  than  those  in  other 
cyanobacteria,  a  characteristic  which  may  be  the  result  of  natural  selection  for  the  ability 
to  withstand  pressure  and  maintain  buoyancy  at  depth  (Walsby,  1992).  Critical  pressures 
for  gas  vesicles  in  T.  thiebautii  (3.7  MPa)  and  T.  contortum  (3.4  MPa)  were  higher  than 
that  of  T.  erythraeum  (1.2  MPa)  (Walsby,  1992).  This  variation  in  gas  vesicle  strength 
does  not  fall  along  phylogenetic  lines,  so  either  there  is  variability  within  clade  I,  or  due 
to  the  overlapping  physical  characteristics  of  T.  erythraeum  and  T.  thiebautii ,  they  were 
misidentified.  It  is  difficult  to  resolve  species  identification  in  the  absence  of  genetic  data. 
The  interpretation  of  pigment  and  gas  vesicle  data  is  unclear,  but  inter-clade  competition 
may  have  resulted  in  the  selection  for  characteristics  which  would  enable  clade  I  species 
to  survive  at  greater  depths  than  clade  II. 

Future  studies 

Understanding  the  factors  that  influence  the  distributions  of  Trichodesmium  spp.  enables 
us  to  predict  the  contributions  they  make  to  the  global  marine  nitrogen  (N)  and  carbon 
(C)  cycles.  Current  ecosystem  models  incorporating  N2  fixation  are  often  based  on  Tri¬ 
chodesmium,  particularly  T.  erythraeum  (clade  I).  However,  the  dominant  phylotype  in 
the  open  ocean  is  clade  II.  If  these  two  types  of  Trichodesmium  are  significantly  different 
in  their  niche  spaces,  the  models  may  not  accurately  represent  the  ecologically  impor¬ 
tant  species  of  Trichodesmium.  Future  studies  should  aim  to  define  the  differences  among 
species  of  Trichodesmium  and  identify  the  mechanisms  which  enable  them  to  coexist. 
Physiological  experiments  should  be  extended  to  include  T.  contortum  and  representatives 
of  clade  II.  The  responses  of  the  two  clades  to  nutrient  limitation  is  under  investigation: 

P.  D.  Chappell  has  been  investigating  Fe  (Chappell,  2009),  and  E.  D.  Orchard  has  been 
investigating  P.  Previous  experiments  that  focused  on  T.  erythraeum  should  be  extended 


146 


to  include  more  species:  chromatic  adaptation  in  response  to  light  quality  (Bell  and  Fu, 
2005),  growth  and  INF  fixation  rates  in  response  to  temperature  (Breitbarth  et  al.,  2007), 
and  the  possible  effects  of  increasing  atmospheric  CO2  (Hutchins  et  al.,  2007).  Ecosystem 
models  should  incorporate  the  diversity  of  Trichodesmium  and  test  whether  the  differ¬ 
ences  between  the  two  clades  are  ecologically  important. 

To  be  able  to  define  the  niche  space  for  each  species  of  Trichodesmium  and  under¬ 
stand  the  mechanisms  which  enable  them  to  coexist,  culturing  efforts  must  be  extended 
to  encompass  the  full  diversity  of  Trichodesmium  spp.  The  1-80  melting  type  of  clade  I 
found  in  Chapt.  4  had  a  deeper  range  than  the  rest  of  clade  I.  It  also  resembled  sequences 
of  Trichodesmium  field  samples  found  at  75  m  (Janson  et  al.,  1999).  We  should  target 
deep  Trichodesmium  for  culturing.  Higher  diversity  within  the  culture  collection  leads  to 
the  possibility  of  developing  a  species-specific  method  for  enumerating  Trichodesmium. 

My  thesis  has  been  a  step  toward  investigating  Trichodesmium  at  a  higher  res¬ 
olution,  both  spatially  and  taxonomic  ally.  While  we  are  far  from  resolving  Hutchinson’s 
“plankton  paradox,”  we  are  beginning  to  understand  that  there  is  important  diversity  within 
the  genus  Trichodesmium. 
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